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The Schumann-Runge O, Emission Bands 
in the Region 3100 A.-2500 a. 


By M. W. FEAST 
Imperial College, London * 


MS. received 19th October 1949 


ABSTRACT. Rotational and vibrational analyses of the emission Schumann—Runge Oy, 
bands in the region 3100 a. to 2500 a. are presented. Tables of wavenumbers are given for 
fine (il, UD), () 1D), Gh aD, (Oso), GIO), GE Oy COE Sy (aay Cs Dy lorravely, able 
origins of the bands, the rotational differences, and constants for the levels v’=7 to 11 are 
given. ‘Three bands, found in a high voltage arc in oxygen between platinum electrodes, 
are attributed to PtO. 


Sie ENE ROD CA LON: 


N a previous paper (Feast 1949) an investigation of the Schumann—Runge O, 
] bands as emitted by a high voltage arc in oxygen was reported. ‘The pressure 
dependence of the bands was discussed and an explanation for this was 
suggested. It was concluded that there is no predissociation in the upper state 
of the system. During this work it was found possible to extend the system 
throughout the quartz ultra-violet region of the spectrum and to identify 17 new 
bands (see Feast 1948 and also the plate given in Feast 1949). ‘The present paper 
gives details of the rotational and vibrational analyses of those of the new bands 
which lie in the region 3100. to 2500 a. (except the (2, 10) and the (0,9) bands). 
Above 3100 a. only one new band (the (2, 22) band) has been found and no measure- 
ments of this band have yet been made. Although four new bands have been 
found at wavelengths shorter than 2500 a. no exact measurement of them has yet 
been attempted. Mr. W. R. S. Garton and the author have recently obtained 
absorption spectra from oxygen at high temperatures in this region of the spectrum 
and it has therefore been decided to defer any further measurements of the 
emission spectra until the absorption work has been completed. 

Since the Schumann—Runge bands have been identified in the solar spectrum 
(Babcock 1945) and since rocket experiments are enabling this spectrum to be 
extended further and further into the ultra-violet it appears important that ~ 
measurements of the Schumann—Runge bands in the region covered by the present 
investigation should be available. 


* Now with the National Research Council, Ottawa, Canada. 
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§2. ANALYSIS 

The analysis of the Schumann—Runge bands by Lochte-Holtgreven and 
Dieke (1929) extended over the region 3100a. to 42004. Tentative analyses of 
some of the bands at wavelengths shorter than 3100a. have been published 
(Millon and Herman 1944, Lal 1948); both these analyses appear however to. 
be in error (see Feast 1948, 1949). No other work appears to have been published 
on the analysis of the emission spectrum. Lochte-Holtgreven and Dieke gave 
data on the vibrational and rotational constants for v’ =0,1,2 and v” =12 to 21. 
Analysis of the absorption bands by Curry and Herzberg (1934) and Knauss and 
Ballard (1935) has provided data on the constants for the levels v’ =0 to 15 and 
vw =0 and 1. 

Previous work has shown that the bands are emitted by the transition 
3y,, >82, and thus each band has six branches (three R and three P). Further- 
more, since the nuclear spin of oxygen is zero all the even K value lines are missing. 
On the dispersion used here the spin triplets are not resolved and hence each 
band appears simply as consisting of a single P and R branch. ‘The spectrograms 
used for the present work were taken on a Hilger E.1 spectrograph which has a 
dispersion of about 56cm~1/mm. at 31004. and about 42cm™!/mm. at 2475 a. 
The region 3100. to 2475 a. was measured accurately against iron arc lines, the 
wavelengths of the latter being taken from Burns’ Iron Arc Tables (1913). 

Table 1 gives the wave-numbers and quantum numbers K for the bands 
analysed; they are: (1, 12)—part only of the band since too many blended lines 
occur in this region of the spectrum for a complete analysis to be made—(0, 11), 
(1, 11), (0, 10), (1, 10), (1, 9), (2, 9), (1, 8), (2,8), (2,7). There are too many blended 
lines in the region of the (2, 10) and the (0,9) bands for a detailed analysis to be 
made. ‘The intensities of the lines are given on a visual scale of 0 to 10 and some 
notes (on blends etc.) are included. ‘The measurements should be accurate to: 
about 0:2cm™1. 

Since the triplet splitting is not resolved the band may be treated as a 
‘X12 transition for the present analysis. The formulae for the rotational 
term differences in the upper and lower states are then (using standard spectro- 
scopic nomenclature) as follows: 

R(K) — P(K) = F'(K +1) — F’'(K—-1)=A,F'(K) =(4K +2)[B’ +2D'(K?+K +1)] 
and 

R(K—1)—P(K +1) =F"(K +1)— F"(K-1) 

=A, F"(K) =(4K +2)[B” +2D"(K?+K +1)]. 

The last parts of these equations involving the constants B and D are only valid 
if higher terms than DK*(K +1)? in the power series for F” and F’ are negligible 
as they probably are in the present case. In Table 2 the differences A,F’(K) for 
the levels v’=7 to 11 are given. The values are the means of the differences 
obtained from the various bands with a common lower vibrational level. The 
agreement of common differences from different bands showed that the measure- 
ments were within the expected accuracy (0-2. cm). 

The values of B” were obtained graphically by plotting A,F’(K)/(4K +2) 
against K and extrapolating to K =0 (see above equation). The values of B” are 
given in Table 3. They are plotted on the graph, on which the values of Lochte- 
Holtgreven and Dieke (1929) for v”=12 to 21 and Curry and Herzberg (1934) 
for v" =0 and 1 are also shown. It will be seen that they fit quite accurately to 


O, Emission Bands 55r 
the straight line and so may be expressed in the form B,” =B,’ —a(v+4). The 
slope of the present line gives « =0-0157 cm-! whereas Curry and Herzberg (1934): 
gave «=0-0158cm™ the difference being within the error of the experiments. 
The value of B,” obtained from the experiments is 1-447cm~!._ The constant D” 


varies little with v" and is usually taken as equal to D,” where D,” = —4(B,”)3/w,2. 


The value D,” = 4-88 x 10-® cm~ has been given previously by Curry and Herzberg 


(1934). Curry and Herzberg give the following formula to represent the band 
origins (vp) of the system, 
¥9 = 49357-55 +(698-86v' — 10-708v'2 — 0-19403) 
— (1568-330" — 11-9930"? +.0-0517v"3 — 0-001430"4). 


1-40 >< Previous Emission 
© Present Emission 
@ Absorption 


BS (cm*') 


Values of Bu” found by various observers. 


The values of v, for a band may be calculated from the wave-number v of any 
line of the band by the formula 

vy =v) +(B’ +B") M +(B’ —B’ +D! — D’)M? +2(D’ +D") M3 +(D' — D") M4, 
where M=K +1 for the R branch and M=—K forthe P branch. In the present 
case the values for the origins were calculated in this way and it was found that if 
M is less than about 15 the terms in D’ and D” become unimportant, and hence 

vg =v —(B’ + B")M—(B' — B") M?. 

The result of this calculation is shown in Table 4 together with the values: 
calculated according to the Curry—Herzberg formula. The latter values seem 
to be consistently slightly higher than the measured values. However this is. 
probably within the accuracy of the formula since the term in v™ has a value of 
about 10 to 20cmz! and further terms would be necessary to give the highest 
accuracy. 
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32200°9 
32095°8 
S1OOIES, 


33601 -2 
835982:5 
33460-0 
33383-0 
$3308e 
33220°9 
33132-3 
S304 2a 
32947°8 
328483 
32746°8 
32640-2 
32530-7 
32415-7 
322979 
32181°8 


5 
6b 
3 
2 
8b 
4 


VP WODUAUADAWOUIMAAWOIN A, 
To a 


O, Emission Bands 


Table 1 (cont.) 


(0, 10) band 


P branch (v) I Rbranch(v) I K 
34824°8 3 OF, 
34813-7 1 34834:6 6 29 
34797-6 = 5 34822:°9 6b 31 
C477 5° ae 34809-6 = 6 33 
34755:7 6 34793 °3 5 335) 
34728:9 6 34772:9 6 37 
346986 86 34749:0 = 7 39 
34664:6 37 34721-2 7 41 
34626:6 6 34689°5 7 43 
SYN, 37/ 34654:-4 7 45 
34540:3. 0 7 34616°1 7, 47 
34490-4 8 SEI a 

(1, 10) band 

P branch (v) I R branch (v) JI K 
35500:7. 4 SSO) — 7b) 39 
35484-7 4 35508°8 7b 41 
35465:3. 5 35495:2 7b 43 
35441-6 = 5 354780 397 45 
354146 5 35457:9 6 47 
353836 6 35433°5 7b 49 
35348:5 6 35404°3 7 51 
353099 6 SYS 7/85 7 53 
B5267/-08= 6 Seson) 7 55 
35220°6 6 35295:-4 6 oi 
SS WOK Fis SDSS 6 59 
SSG Wes 3520550 Neb 61 
5058225 7; Soul Sule} 63 
34996:8 8 35096°1 a, 65 
34930°8 8 $5036.78 67 
34861°5 8 34973:2 8dv 69 
BA 8) en, 34906-1 8 

(1, 9) band 

P branch(v) ZL Robranch(v) I K 
368616 386 31 
36850:0 O 36866:7 5 33 
SOOO — S 36857°5 5 35 
3681375" 5 36843 °3 6 37 
36789:4 5 36826:0 6 39 
36761:1 7b 36804°3 6 41 
86729-2556 367788 6 43 
36693:2 5 36749-1 6 45 
KGHSSIG 7 OVI © 47 
36609°5 47 36678:6 8b 49 
36561:9 7 36636°8 8b 51 
36509°5 = 6 36590:°9 7b 53 
36454:3. 7b 36541:2 6 55 
36394-1 6 36488°1 6 


P branch (7) 
34438-0 
34381-3 
34321-4 
34257°5 
34190°-1 
34118:2 
34043 -7 
33964-7 


P branch (v) 


34710°9 
34629°4 
Sos el 
34456°3 
34364°3 
34267-2 


34064-2 
8305025 
33844°9 
SOTAIeT 
33609°8 
33486°7 
33361 °3 
33230°3 


P branch (v) 
36330°5 
36259-6 
36191-2 
36115°8 
36035°7 
359522 
35864°8 
BS 77/103) 
35678°7 
35579°6 
35475°7 


If 
a 
6 
7 
i, 
6 
7b 
6 


7 do 


IR branch (7) 


8 


ND 


R branch (v) 


34526°6 
34476°3 
Sasi 
34364°3 
34303°5 
34238 -3 
34169°-7 
34097°3 
34020°8 
33940°6 
33856°5 


34836°3 
34759°6 
34681°1 
34598°5 
345011 
34422-1 
343286 
34228°5 
34128-2 
34020°8 
8391222 
33796°3 
33677-5 
33560°3 
33434°8 
3530923 


R branch (7) 


36430°3 
36368°5 
363034 
SOAS) 
36159°6 
36083 -0 
36000°4 
35916:2 
35826°7 
35732-4 
35635°5 
BI5S3e9 
35428°7 


553. 


554 


K Pbranch(v) I R branch (v) 


837/536 7) 
37498 -1 
37477°5 
374530 
37424-2 
373914: 
37354°6 
BUSES! 
372684 
37219-4 
37165°9 
37108°7 


P branch (vy) 
38204°3 
38183°7 
38159°5 
38130-4 
38097°5 
38060-7 
38019°-3 
37974°3 
37925°1 
37872:0 
37814-4 
37752°5 


P branch (v) 
38899°8 
38885°-9 
38869°1 
38847°-9 
38822:-7 
38793 °3 
38759°8 
38717°-4 
38679°6 
38633 :2 
38582°6 
38527:°9 
38468 -6 
38405°9 
383386 


DNDNDNDADADAHDAAY $F fF W 


i 
4 
4 
6 dv 


Unb NAPUNUN 


NNNNNDNNDDDMNUMADBAYN WY 


S7o21e0 
SIO OT 
37489 -0 
37466°3 
37440-0 
37409 °3 
SSA 
37335°8 
SPT) 
372453 
37194-6 


R branch (v) 


38221-0 
38191-1 
Sol753 
38148 -2 
38116°5 
38081-6 
38042-7 
37999 -6 
37952-6 
37902-0 
37846:-2 


R branch (v) 
38909-9 
38903°5 
38892:-7 
38877:9 
38858°8 
388356 
38807:-9 
38776:3 
38740°5 
38700-4 
38656-0 
386074 
38554°5 
38497-3 
38436-2 
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Table 1 (cont.) 


(2, 9) band 


NN 
a 


DDANUMADOUA un 
alts, 
— 


(1, 8) band 


(2, 8) band 
I 


oO. 
Ww 
NI 


DDANDDADDADAUDADAUDYS 4, 
aS 
NI 


P branch (v) I Rbranch(v) J 


37047°6 
36982°1 
36912-9 
36839°0 
367611 
366786 
36590°9 
365041 
36409°8 
36312-0 
36209°7 


P. branch (7) 


37687°3 
37617°-4 
37543°8 
37466°3 
37383°7 
SIL 
37208 -6 
S7AISEO 
37012-4 
36912°9 


P branch (v) 


38266°7 
38191-1 
Sol TZ 
38027°3 
37939-0 
37846:2 
37749 -6 
37648 -4 
37543°8 
37434-0 
37320°6 
37202°8 
37080°-4 
36954°7 


I R branch (v) 


5 
6b 


3:7139°2 
37080°4 
37016-0 
36949-0 
36877°4 
36800-9 
36720°9 
36636°8 
36549°3 
36454°3 
36360°6 
36262-7 


R branch (v) 


37786°7 
377234 
37655°7 
37583°1 
37507°7 
37427-8 
31343°9 
37255°4 
37163 -4 
37067 -7 
36966-7 


383708 
38300-6 
38227-0 
38148-2 
38066-5 
37979-0 
37888-9 
377939 
37694-7 
37590-6 
37482-4 
37370°6 
37255-4 
37136-0 
370124 


6 
7b 
Wf 

4 do 


MWwWARUONDAD wn, 
Qu. 
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Table 1 (cont.) 
(2, 7) band 

K Pbranch(v) I Rbranch(v) J K Pbranch(v) J R branch (v) J 
3 40293:-4 O 21 40021:2 5 40088:0 6 
5 40280:6 3 40298-0 5b 23 39969°3) 5 40042:7 5 
i 40262:9 4 40286:7 3 25 39913:0 6 39992°5 5 
9 40241°5 5 40271-4 5b Gl 39852:5 5 39938°1 5 

11 40215:9 4 40251-6 5 29 39787:0 35 398792 5 

13 40185:8 5 40227:7 5 31 39718°1 5 39815°3 5 

15 40150°8 5 40199°5 5 33 39644:7 5 39748 -4 5d 

iy 40112:0 5 40166-6 5 35 39566:5 6d 39676°3 5) 

19 40068:8 5 40129-6 5 37 39600:1 5 

No measures beyond P (35) and R (37) 

Abbreviations in Table 1 :— 

b—line blended. dr—tline diffuse to red. 
b*—blended with a Pt line. dd—line very diffuse. 
d—line diffuse. do—line probably double. 
dv—line diffuse to violet. +—line on diffuse patch. 
Table 2. A,F’(K) 

TG =O = 9 re) 1) ed Ke Ue Ue On Ua eae oa lel 
4 24-0 DD 34 181:9 179°:7 177-2 174:7 172-9 
6 35:1. 34:4 33:1 34:3 33-3 36 189-4 187°8 185:5 182-3 
8 45:2 44:8 44-1 440 42:8 38 199:6 197-6 194:9 193-0 

10 5525 D5e2 D4 Om S88 53-2 40 209°5 207:4 206-0? 202-9 

i 65:8 65:5 64:9 63:9 63°3 42 Mier, ilges  pileto’y ~ Dilpcil 

14 76:9 75:8 75:0 74:3 73-3 44 229-9 227:0 224:8 221-4 

16 87:5 87:5" 85:5 84-7 83°8 46 240°5 237:4 234:2 231-4 

18 97°38 97:0 95:8 94:5 93-9 48 250°3 247-1 243-9 241-0 

20 108-4 107°3 106-1 104:°5 103-6 50 260°7 256:7 249:7? 

De AS 7 ails OcO eld5- Ol dl 327 52 270-0 264:4 260-1 

aoe Ome led 25-0 124-0 54 279°6 272:0 271-0 

26 140-0 138°5 136°6 135-1 133-9 56 290:2 283°3 

28 151-1 149:1 147-1 145-1 143-5 58 300:7 293+1 

30 161-1 158°8 157-4 155-0 153-4 60 302°4 

32 170-6 169°4 167:5 165:0 163-4 

Table 3. New Values of B,,” 
v’ 11 10 9 8 7 
Be 1:269 1-282 1:298 1:316 1-329 
Table 4. Wave-numbers of Band Origins 

Band Formula Experimental Band Formula Experimental 

Db, Uf 40305:9 40304:6 iL 10) 35524:1 35522°7 
38911-2 38910°1 0, 10 34836-2 34835°7 


Poy ies 

1,8 38245°8 38245-2 ils ils 34197:2 341973 
Le 36873°7 36872°8 Omit 33509:2 335089 
Dee) 3/539°0 SSCS) Pez 328929 32892-2 
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§3. THE VIBRATIONAL INTENSITY DISTRIBUTION IN THE 
SCHUMANN-RUNGE SYSTEM 

The bands found in the present investigation and in that of Lochte-Holtgreven 
and Dieke indicate a wide Franck—Condon parabola in the v’, v” array as would be 
expected since the emission of the system involves a large change in equilibrium 
inter-nuclear distance. Miss M. E. Pillow has recently been engaged on a 
theoretical treatment of the vibrational intensity distribution in this system. 
Though her results are not yet quite complete they are sufficiently extensive for 
it to be stated that there is good agreement with the present experiments. 


S4ee BANDS Ach RB Ulin ED ah Omet® 


Whilst working with a high voltage arc in oxygen between Pt-Rh electrodes. 
(see Feast 1949 for a-discussion of suitable electrodes) a few bands were found in 
the visible region. ‘These bands also occurred when using Pt electrodes but not 
with Cu electrodes. In the latter case the CuO bands appeared. ‘The bands are 
therefore attributed to PtO. The wavelengths of their heads are 56644.; 
59064. and 6180a. They are degraded to the red. There appears to be no 
previous record of bands due to this molecule. 
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MS. received 19th October 1949 


ABSTRACT. The spectrum of an electrodeless discharge in pure oxygen has been 
studied. New bands of the O,* Second Negative system have been found and the intensity 
distribution in this system is discussed; subsidiary parabolae as well as the main Franck— 
Condon parabola are obtained as is predicted by the theory. Emission bands reported by 
Johnson as due to O 3 are due to O,*+ and O11. Some abnormal intensities of O11 lines 
excited in the electrodeless discharge are noted. The electrostatic and electromagnetic 
types of excitation are briefly compared. 


See No ROD, Carl ON 

HE O,* Second Negative band system was excited by an electrodeless. 

discharge in pure oxygen. In addition to bands previously reported, 

others were found which have been shown to belong to this system. When 
the new bands are taken into account the system is found to show a vibrational 
intensity distribution with secondary and tertiary parabolae as well as the normal. 
Franck—Condon parabola. ‘These subsidiary parabolae are to be expected for 
the transition since it involves a large change of equilibrium inter-nuclear 
distance. 

Interest in these bands was originally stimulated by Johnson’s observation 
(1924) of weak bands, which he attributed to O3, on his spectrograms of the 
O,* Second Negative system. ‘These bands are not identical with those observed 
by Janin (1938) in emission from an ozonizer discharge, and also attributed 
to O;. It was thus felt desirable to re-examine the oxygen spectrum in order 
to reach a more definite conclusion with regard to the O, emission spectrum. 
The present work with an electrodeless discharge has revealed no bands that 
can be attributed to O;. A study of Johnson’s results has shown that his ozone 
emission spectrum may satisfactorily be accounted for as Or lines and O,* bands. 
It has not been possible to decide whether or not Janin’s bands are a true O3. 
spectrum. Finally the relative intensities of the Or lines excited in the high- 
frequency discharge are discussed, and a complex spectrum obtained from 
another form of high-frequency discharge is mentioned. 


§2. EXPERIMENTAL 


The spectrum emitted by a glow discharge in oxygen often consists mainly 
of bands and lines due to impurities in the gas (NO, N,, CO,+, CO). To excite 
band spectra from oxygen itself it is necessary to remove all sources of impurity 
from the vacuum system employed. In the present investigation a fused quartz 
discharge tube 15 cm. long and 3 cm. in diameter with a quartz window fused 
on one end was used. This tube was evacuated by a mercury diffusion pump: 
with liquid air trap and a ‘Hyvac’ pump. Pure oxygen was obtained by heating 


* Now with the National Research Council, Ottawa, Canada. 
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potassium permanganate in a side tube. After the apparatus had been baked 
.out and swept out once or twice with pure oxygen it was found that the spectrum 
was free from all impurities except faint Mn raves ultimes and one or two faint 
‘SiO bands. 

To excite the spectrum a Metropolitan—Vickers induction furnace was used. 
‘The output of this generator was applied to two external electrodes in the form 
of tin-foil bands wrapped round the outside of the discharge tube. The region 
2000 a. to 7000 a. was investigated with a medium quartz (Hilger E.2) spectrograph 
using Ilford ‘Astra III’ and ‘Q’ plates. A prism grating instrument was 
employed to examine the region 6000a. to 10,0004. using Kodak sls awe Oe 8a 
and ‘1.Q.’ plates. These plates were hypersensitized in ammonia before use. 


§3. APPEARANCE OF THE SPECTRA 

At a pressure of about 0 05 mm. Hg a yellow-green glow is obtained, and this 
.consists of bands due to the O,+ First and Second Negative systems together 
with lines due to Ou and weaker ones due to Or and Orr. This discharge fills 
the whole of the tube except for the region inside the electrodes, where the 
luminosity is confined to the axis of the tube. At a slightly higher pressure 
the discharge is non-luminous except for a mauve glow inside the electrodes 
and strongest at the walls. ‘The glow shows the same main features in its 
-spectrum as that shown by the lower pressure discharge, but the lines are stronger 
relative to the bands. At a still higher pressure just prior to extinction of the 
discharge there is a faint glow along the axis of the tube between the electrodes. 
This last type of glow gives a continuous spectrum. ‘The continuum was 
photographed with a small quartz (Hilger E.3) spectrograph because of its weakness 
_.and was found to extend from the blue to the ultra-violet region with gradually 
decreasing intensity. 


§4. NEW O,+ SECOND NEGATIVE BANDS 

The Second Negative system of O,* consists of a large number of red-degraded 
-double-headed bands extending from about 2000a. to 60004. The bands have 
‘been studied by several workers (Johnson 1924, Ellsworth and Hopfield 1927, 
Stevens 1931, Mulliken and Stevens 1933, Lal 1948). This work has established 
that the bands are due to the transition a*II1,—x?II, with a doublet separation 
‘of about 200 cm. It has been tentatively suggested (Lal 1948) that the 
Franck—Condon parabola of this system may be much narrower than has 
previously been assumed. However, this suggestion is not supported by the new 
‘data given here, which are consistent with the large change in equilibrium 
inter-nuclear distance involved in the transition (7, =1-41a., 7,” =1-14a.). 

The present spectra show, in addition to the bands previously reported for 
this system, a considerable number of new bands. These new bands are red- 
degraded and show a 200 cm™ doublet separation. Hence an attempt was made 
to fit them into the Second Negative system. It was found that all the new 
bands could be satisfactorily accounted for in this way. Head measurements 
are presented in Table 1. 

In Table 1 A is the wavelength of the head measured from E.2 spectrograms; 
over most of the region investigated this should be correct to 0:5a. The wave 
numbers obtained from these wavelengths are given in the third column. The 
vibrational quantum numbers assigned to the bands are given in the last column 
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In this column i and ii denote the shorter and longer wavelength heads of each 
‘band. The analysis was made on the basis of wave-number calculations from 
‘the constants of the system. The calculated values given in the column are v 
the formula used for these was 


vy = 38596 
38796 


cal» 


} +{898-9(0" +4) —13-7(o" +4)%} 
— {1876-4(v" +4) —16-53(v” +4)7}. 


Table 1 

m\ il v Veal vv" A I v Veal vv 
2562°5 2 39018" 39012" 115 4 3257°6 3. 30689 30691 8, 8i 
2588-9 4 38615 38615 10,41 SSI DE7 aS O17 Seem OLSZ By wi 
.2602:9 4 38407 38415 10,4 ii 3334-2 a 29984 29982 Se7 i 
2677°8 1 S/SSs0 TOSS 8, 41 ISIS) Se XSL 29812 7,811 
.2691°5 1 37143 37138 8, 4i1 3465-9 1 28844 28844 10,101 
2860-0 3. 34955 34947 io 3542-0 4 28225 28219 9,101 
2877-3 3. 34745 34747 Uo 3567°7 4 28021 28019 9,10 11 
2936-9 4 34040 34040 6, 5 il 3653-0 Pe SLY SG) LOA 8, 10 i 
3005-0 2 33268 33269 On 3700°5 a 27016 27019 3,81 
3022°:8 2 33072 33069 7, O11 3928-9 a 25445 25445 SLO; 
3070-0 5 32564 32562 6,61 3959-1 a 25251 25245 55 Oar 
3088-7 De 52307 SL S02 6, 6 il 4760 3° WOO — AUOOE Sel Sar 
3096 ? Dy SSI 32303 Saal 4803 3 20815 20806 yy Hg a8 
3113-6 2 SONOS Seis} 8, 711 5102 30 1959519592 5,141 
3160-7 a 31629 31624 iia 5143 ? 3 19439 19392 5, 14 ii 
B220°9 a 30991 SOIST wl, Oi 


i, 11 denote shorter and longer wavelength heads. 

+, blend with 3,91 and 3,911 (see text). 

a, difficult to judge intensity due to overlapping structure. 

?, wavelength may be inaccurate due to overlapping structure. 


This formula uses the constants given by Jevons (1932) and Birge (1929), 
which differ from those used by Ellsworth and Hopfield (1927). The new 
measurements do not justify any change in this formula, as deviations from it 
appear to be random and hence must be attributed to experimental error. A 
consideration of band intensities suggests that the bands marked f are actually 
blends of the bands given in the last column with 3,9i and 3,911 (25,440 cm" 
and 25,240 cm-!). The intensities recorded in the second column are visual 
estimates. An attempt has been made to keep them on the scale used by Johnson 
(1924). On this scale the strongest band of the system has an intensity 9. 


§5. VIBRATIONAL INTENSITY DISTRIBUTION 

The Figure shows the distribution of vibrational intensities in the system. 

t is seen that all the bands previously reported fall on a wide Franck-Condon 

parabola. However, when the new bands are taken into account, secondary and 

tertiary parabolae are found. Observation of these extra loci is in accordance 

with wave-mechanical theory, which predicts them for transitions involving 

a large change in equilibrium inter-nuclear distance (see for instance Gaydon 

and Pearse 1939). It will be noted that there are bands falling on the subsidiary 

loci which have not been observed. Calculations show that all these missing 
bands would be obscured by main parabola bands. 
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Miss M. E. Pillow has kindly supplied details in advance of publication 
of her theoretical work on the intensity distribution in this system. These 
calculations were made on the assumption of an equal distribution of molecules: 
amongst the vibrational levels of the upper state (i.e. an infinite temperature 
distribution). Comparison with the results of the present experiments shows 
that good agreement is to be expected if the distribution of molecules in the 
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Distribution of vibrational intensities in the Second Negative system of O,*. 


upper state is changed to weight the higher wv’ levels even more heavily that in 
the infinite temperature distribution. Such a distribution is in fact likely if 
the system is excited by electron collisions with molecules in the ground state 
of either O, or O,*. A consideration of the various potential curves shows that 
if the Franck—Condon principle is obeyed the higher v’ levels of the O,+ upper 
state will be preferentially populated. ‘lhe agreement of theory and experiment 
may therefore be stated to be satisfactory. 


$6. THE REPORTED ©, EMISSION SPECTRUM 


On Johnson’s spectrograms of the O,+ Second Negative band system (1924) 
appear a number of bands which were not fitted into any system. From a 
comparison of the wavelengths of these bands with the wavelengths of the ozone 
absorption bands Johnson concluded that these formed an O, emission spectrum. 
On the other hand, Janin (1938) records a large number of bands emitted by an 
ozonizer discharge in oxygen. These bands are attributed to O., but they are not 
identical with those observed by Johnson. No structure was found on the present 
spectrograms which could be interpreted as due to O3. A study of Johnson’s 
measurements and his published photographs has shown that his ozone bands 
may be accounted for as O,+ bands and Om lines. The identification of 
Johnson’s bands is given in Table 2. It should be noted that the O1 lines quoted 
were all found to occur together with the bands on the present spectrograms. 

In two cases one O,+ band is used to identify two ‘O,’ bands. This is 
justified by the appearance of Johnson’s spectrograms, two positions on the band 
having been measured. Some confusion exists about the wavelengths of the 
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two ‘O;’ bands marked ¢. The wavelengths 4333-3. and 4316-54. appear 
in Johnson’s list, and if they are indeed correct they may be identified as the 
Ou lines in Table 2. However, on the published spectrogram a band is marked 
4316a. which actually lies near 4240a. Another band is marked 42164. and 
appears quite strong, though it is not to be found in the table of wavelengths. 
A comparison of spectrograms indicates that these last-mentioned bands are 
probably 4238-0 1,9ii and 4218-3 1,91 O,+. Lines are present above 4300. in 
the photograph which are probably the O1 lines given in Table 2, but these are 


Table 2 

O,; band if Identification O, band I Identification 
4464-7 1 4465-45(4) O11 3337-8 6 : 
4333-3 2 4336-86(6) Ort 3331-5 a S08 eS aot On 
4316-5 5 4317-16(8) Ort 3313-1 5 3314-3 5,71 O,t 
4276-9 2 4275-52(4) O11 3284-6 2 3287-59(9) O11 
3568-9 2 . : 3278-1 2 3277-69(7) O11 
3565-2 i See Oh ae Oe" 3606-1 5 ? 
3465-3 0 3465-9 10, 10i Os* 3097-5 4 3093-1 8,7i Ost 
3378-3 1 3377-2(7) O11 3090-4 8 3090-8 6,6 ii O,t 

inseeatexts 


not marked ‘O,’. ‘The present interpretation of Johnson’s bands still leaves 
one of them unclassified. This band does not occur in the present spectra, and 
may possibly be due to an impurity. 

It is more difficult to be certain of the origin of the bands reported by Janin 
(1938). Comparison of the wavelengths shows that many may be due to O,+ 
and O1r. However, some strong O,+ bands and O11 lines do not appear in Janin’s 
list, and thus the correctness of this identification is in doubt. 


§7. NOTES ON INTENSITIES IN THE SPECTRUM O1 


Emeléus (1934) has pointed out that the following multiplets which are 
present in the positive column of a condensed discharge are missing in a negative 
glow of low current density: 


4956 a.4941 a.(3P)3p2P° —(3P)3d2D 
4406 a.4359 a.(3P)3p 2D° — (2P)3d2D 
2419 a.—2407 a.(3P)3p 2D° — (#D)3d2D 


An examination of the high-frequency discharge spectrograms shows that these 
multiplets are missing there; it has revealed that a number of other multiplets 
are also missing or weakened. ‘The multiplets are listed in Table 3. since 
there appears to be no fundamental difference between a line being weakened 
or completely missing from a particular spectogram, the two types are not 
distinguished in the table. Lines with a common initial level are grouped 
together, beginning with those lines whose initial levels are in common 
with the above three multiplets. The identifications were made using the 
data of Mihul (1928) and Moore (1933). For each initial level given in Table 3 
all known multiplets were studied and they are all listed in the table, showing 
that there are no inconsistencies in the data. The multiplets marked + are in 
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any case rather weak, and to find them missing on moderately exposed spectro-- 
grams is not, of itself, surprising. They are, however, listed for completeness. 
since each has a common initial level with a multiplet which definitely shows. 
the effect. 


Table 3 
ACA.) d(a.) 
+4871-4861 (:D)3p?P°-(:D)3d?D 4448-4443 ({D)3p?F°-(?¢D)3d?F 
44332-4327  (1D)3p2D°-(4D)3d2D 4703-4698  (1D)3p2D°-(:D)3d2F 
2575-2572 (?P)3p2P°-(-D)3d?2D 2530-2518 (?P)3p2D°—-(??D)3d?F 
4705-4753 (?P)3p2D°—(@P)3d?F 4169-4121 @P)3p7P°=CP)sd-& 
4112-3876 (?P)3p4*P°—(?P)3d2F 4925-4891 (3P)3p4S°—(?P)3d4P 
3857-3833  (3P)3p4D°-(*P)3d2F 3097-3864 — (@P)3p4D°-(2P)3d4P 
4699-4751 (?P)3p?D°—(P)3d4D +3488 (?P)3p?S°—-(@P)3d4P 
4865-4856 (?P)3p4S°—(?P)3d4D 3390-3377 (?P)3p7S°—(?P)3d?P 
4121-4097 (?P)3p4P°—-(?P)3d4D 5207-5160 (?P)3p7P°— (2P)3d2P 
3883-3848 (?P)3p4*D°-(?P)3d4D +3985 (?P)3p *P°—(?P)3d?2P 
43475 (2P)3p2S°-(*P)3d4D 


+, weak multiplets (see text). 
In the multiplet (?P)3s?P—(P)3p2D° the line 4452-4 a. is abnormally weak. 


There is no very obvious cause for these irregularities. Emeléus (1934) 
discussed various possible explanations for the behaviour of the multiplets he 
examined, but reached no very satisfactory conclusion. However, as he pointed. 
out, it is of importance to recognize these irregularities even if no explanation 
can immediately be found for them, especially since O1r lines are met with in 
astrophysical work (stellar spectra and possibly the aurorae). 


S35 DHE SOLENOID MEDH OD FOREX Clr VATTON 

The work described in the preceding sections was carried out using external! 
electrodes as described in §2. Some experiments were also carried out using 
a solenoid wound round the discharge tube as the method of applying the high-. 
frequency impulse to the gas in the discharge tube. At about 0-05 mm. Hg of 
oxygen a very intense, whitish-blue discharge is obtained. The spectrum of 
this discharge contains lines of Or and O11 but, in contrast to the previous method 
of excitation, the Or lines predominate. O,* First and Second Negative bands. 
also occur. ‘The Second Negative bands have been observed in detail and show 
a much higher rotational temperature than in the electrostatic method of 
excitation; this leads to an almost unrecognizable appearance of the spectrum 
at the short wavelength end. There also appears to be a lowering of the: 
vibrational temperature. 

It is hoped to make a more detailed comparison of the two modes of excitation 
in the near future. 
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The Spectra Emitted by the High Voltage Arc in Nitrogen,. 
Hydrogen, Nitrogen-Hydrogen Mixtures and Ammonia 


By M, Wa FEASL* 
Imperial College, London S.W.7 


MS. received 19th October 1949 


ABSTRACT. The spectra emitted by the high voltage arc in nitrogen, hydrogen,. 
hydrogen—nitrogen mixtures and ammonia have been studied. 

‘The appearance of the N, Second Positive bands at various nitrogen gas pressures is. 
discussed and the divergence of rotational and vibrational temperatures is noted. Traces 
of oxygen greatly reduce the strength of the N, bands emitted by the arc in nitrogen. 
In N.—-H, mixtures the NH 3360 a. band is found to be emitted strongly only when hydrogen 
is present at a very low partial pressure. An attempt is made to explain the various 
observations. 


Soe EN ER ©) Ui Cae ON 

N a previous paper (Feast 1949) a high voltage, low current, arc source was 

described and its use in the excitation of the O. Schumann—Runge bands 

was discussed. ‘The investigation of this source has now been extended 
to include several other gases. ‘The present paper reports spectroscopic 
observations on the excitation of nitrogen, hydrogen, mixtures of these two: 
gases and ammonia. A study, such as this, of the variations produced in band 
spectra by changing the experimental conditions has become of interest lately, 
due to the possibility of applying the results obtained to astrophysical and upper 
atmosphere problems. 


§2. THE HIGH VOLTAGE ARC IN NITROGEN AT ALMOSPHE RIC 
PRESSURE 
Preliminary observations of the spectrum of the high voltage arc were made 
with a Hilger small quartz spectrograph which covered the region 2000 a. to 
65004. In the infra-red region a prism-grating instrument was used from 
6000. a. to 9000. Features of particular interest were later examined in detail 
in the first order of a 10 ft. Eagle Mounting concave grating (5-5 a/mm.), 
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When cylinder nitrogen at atmospheric pressure was used in the are without 
purification the spectrum in the ultra-violet region was found to consist of strong 
N, Second Positive bands and weak Nj First Negative bands together with 
NOy, OH, CN and sometimes NH bands. If, however, a piece of metallic 
-sodium was placed in the arc chamber, the NO and OH bands were removed. 
The CN bands (the violet system) were very persistent unless the vacuum system 
was particularly clean, the smallest trace of carbon compounds bringing the 
bands up strongly. It was noticed that the Nj, bands were considerably 
strengthened by the removal of oxygen from the gas. This is of interest in 
connection with the theory of the formation of NO, in discharges. Brewer and 
Westhaven (1930), by chemical methods, and Wansbrough-Jones (1930a), by 
controlled electron experiments, conclude that N;, or Nz (excited), ions are 
necessary for the production of NO,. It cannot be said that this point is definitely 
settled yet—see for instance Willey (1937), and compare with the analogous 
controversy in the case of O, (Wansbrough-Jones 1930 b, Brewer and Westhaven 
1930b, Deegan and Emeléus 1948)—but the present observation links the dis- 
appearance of the Nj excited states with the presence of oxygen in the gas, and 
thus tends to support the Nj theory. It is relevant to remark here that though 
there is copious production of NO, from a high voltage arc in N,—O, mixtures 
(i.e. air), there is little or no production of O, from the high voltage arc in oxygen. 
This is probably due to the fact that any O, that is formed will be almost certainly 

dissociated at the high temperature of the arc. 

In the infra-red region of the spectrum 10-ft. grating spectrograms revealed 
nitrogen First Positive bands. The appearance of these bands differed 
considerably from that found in a low pressure discharge. ‘This may be explained 
as due to the much higher rotational temperature in the arc. 


$35) DAE NITROGEN SECONDS? OST TIVE SEAINIDSs LNG EE a renG Et 
VOLTAGE ARC AT VARIOUS PRESSURES 

The appearance of the nitrogen Second Positive bands in the high voltage arc 
at various nitrogen gas pressures has been studied in detail on 10-ft. grating 
spectrograms. 

The general appearance of the Second Positive bands in the arc at atmospheric 
pressure may be seen from the Plate at g. A spectrogram of a positive column 
discharge in nitrogen at low pressure is included for comparison. The change 
in structure of the bands may best be illustrated by discussing the (0,0) sequence. 
It will at once be seen that the rotational temperature of the bands is considerably 
higher in the arc than in the discharge. A similar result was obtained by 
Nakamura (1927) who employed a low voltage tungsten arc in nitrogen at 30 to 
50 cm. Hg pressure. He used this mode of excitation to measure high 
rotational lines (up to J =112) in the (0,0) band. Microphotometer measure- 
ments of the present spectra give rotational temperatures of about 4,000°k. for 
the are and 1,000°k. for the discharge. 

Further study of the arc and discharge spectra reveals that whilst the 
rotational temperature of the Second Positive bands has increased from the 
discharge to the arc, the vibrational temperature has decreased. Thus in the 
(0,0) progression the bands (3,3), (4,4) are clearly visible in the discharge and 
, entirely absent in the arc spectra. There is, of course, emission from v’ levels +0 
‘since sequences to the short wavelength side of the (0,0) sequence do occur, but 
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N, 2nd Pos. 


NOy e uaa 
a — 76 
b 2 74 
c 10 66 
d 38 38. 
e USS LS 
if 76 — 
asbo 3000 3500 4000 4500 A. 
(1,0) (0,0) (0,1) (0,2) 


is) 


(4,4) (3,3) (0,0) 


a-f. Spectra of a high voltage arc in H.-N, mixtures. Equal exposure times for all spectrograms. 
E.3 spectrograph. Spectrogram f is too weal to show the H, continuum. 
§. Spectrum of a high voltage arc in N, at 76 cm. Hg pressure with a spectrum of a normal low pressure 
discharge in Ny for comparison. 10 ft. Eagle Mounting. spectrograph. 
h-k. Spectra of a high voltage arc in N, at various pressures. Equal exposure times for all spectrograms. 
10 ft. Eagle Mounting spectrograph. 
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the present work-shows that they are weakened relative to the discharge spectra. 
Further evidence of this divergence of rotational and vibrational temperatures 
has been obtained by photographing the arc in nitrogen at various gas pressures. 
The Plate (h-k) shows a set of equal exposures at pressures of 760 mm., 270 mm., 
150 mm. and 75 mm. Hg. _ It will be noted that the strength of the bands increases. 
with decreasing pressure, at first very rapidly. This behaviour should be 
contrasted with the rapid decrease in intensity of the Schumann—Runge O, bands. 
under comparable conditions (Feast 1949). The spectrograms show the increase 
in vibrational temperature and the decrease in rotational temperature as the 
pressure is reduced. 

It may be noted here that the First Positive bands also show this divergence 
of rotational and vibrational temperatures. The bands of this system which 
lie in the visible region correspond to high v’ values and are all quite weak compared 
with the infra-red emissions. 

The explanation of the difference in behaviour of the rotational and 
vibrational structure of the bands is probably to be found in the high collision 
rate at atmospheric pressure. The excited molecules with v’>0 will suffer 
sufficient collisions before radiating (at these high pressures) for a large number 
to be reduced to v’=0. This assumes of course that vibrational de-activation 
is more easily affected under these conditions than electronic de-activation; 
this is probably correct. Since the arc may approach thermal equilibrium at 
the rotational temperature, rotational de-activation is not to be expected to be: 
important. 

Predissociation by rotation has been observed in all levels of the upper state 
of the Second Positive system except v’=0. Coster, Van Dijk and Lameris. 
(1935) have attempted to fix the position of the intersection of the potential 
curve of the dissociating state with that of the dissociated state. From the 
levels in which predissociation has been observed they deduce a point, B, on the 
potential curve of the dissociated state below which this intersection must lie. 
Measurements of the predissociation can do no more than this since the dissociating 
state probably has a maximum in its potential curve at large inter-nuclear distances. 
Coster et al. have examined a discharge in nitrogen and measured lines of the 
(0,0) band to J =90. Since no predissociation occurs up to this level they are 
able to deduce another point, A, above which the intersection of the two curves: 
must lie. They then tentatively place this intersection mid-way between A 
and B on the dissociated curve. This is of course somewhat unsatisfactory, 
especially as the J =90 point is completely arbitrary. They have apparently 
overlooked Nakamura’s work in which J=112 was attained. No further 
experiments have yet been made with the present source, but it might be possible 
to fix the point of intersection of the two curves exactly by photographing the 
(0,0) band from the arc at very large dispersion to obtain a limit for predissocia- 
tion. The very high dispersion would be necessary since the predissociation is 
probably in the region where the (0,0) and (1, 0) bands overlap. 


§4. THE HIGH VOLTAGE ARC IN HYDROGEN AT ATMOSPHERIC 
PRESSURE 


In hydrogen at atmospheric pressure the high voltage arc between Pt-Rh 
electrodes is barely visible, being faintly bluish in colour. Spectrograms of 
this discharge show the Balmer lines, the hydrogen ultra-violet continuum and 
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the H, many-lined spectrum. It was often found that absorption bands appeared 
near 2100. in the continuum. These bands became much stronger when a 
little nitrogen was added to the arc, and they have been identified as due to NH. 
There is apparently sufficient nitrogen in the hydrogen gas cylinders to give 
these bands. 


$5. THE HIGH VOLTAGE ARC IN HYDROGEN-NITROGEN MIXTURES 


The experiments in hydrogen-nitrogen mixtures were all carried out with 
a total pressure of one atmosphere. A series of spectra was obtained which 
covered the whole range from 100°, hydrogen to 100% nitrogen. ‘Typical 
spectrograms for the ultra-violet and blue regions are shown in the Plate (af); 
these were all taken with the same exposure times. 

The 100% cylinder N, spectrum shows NOy, N, Second Positive, OH 
and CN bands (probably also N,+ bands). Introduction of hydrogen (2 cm. Hg 
pressure) removes the NOy bands completely. This is presumably due to 
the hydrogen reducing any NO present. The other alteration in the spectrum 
is the appearance in great strength of the NH3360 band (the 6 band). The 
most interesting observation which has resulted from these experiments is the 
peculiar behaviour of the NH£ band as the percentage composition of the 
gas mixture is varied. The Plate shows that as the percentage of hydrogen is 
increased the band reaches a maximum and then falls in intensity. It is missing, 
or is at least very weak, in a mixture containing 50°% H, and 50° N., which might 
have been expected to show the band most strongly. No accurate determination 
has been made of the pressure of hydrogen at which the band reaches its maximum 
intensity. A small quantity of hydrogen (a few mm. Hg) produces the band very 
strongly, and the intensity does not appear to vary appreciably until it starts 
falling at between 1 and 3 cm.Hg of H,. Observations in the visible region 
with a mixture of 50°, H, and 50° N, reveal the presence of the ammonia « band 
(probably due to NH,). 

With smaller proportions of hydrogen the « band is weaker and is overlapped 
by nitrogen First Positive bands. 

The only other experiments of which the author is aware which are relevant 
to the present observations on the NH band are those of Gaviola and Wood 
(1928) (see also Willey 1934). Gaviola and Wood succeeded in obtaining the 
NH band in a sensitized fluorescence experiment of the Cario—Franck type. 
‘The fluorescence was excited in the normal manner by mixing the gas (H, and N,) 
with mercury vapour and illuminating it with the mercury line at 2537. Atoms 
of nitrogen and hydrogen were formed by collisions with excited Hg atoms 
and these combined to give NH and the Bband. It was found that the NH band 
fell in intensity as the hydrogen pressure was raised. This behaviour was 
attributed to the reaction NH +H,—NH,, which would reduce the intensity 
of the band. Such an argument seems plausible, although quantitative 
confirmation does not appear to be possible. The same reasoning can be applied 
equally well to the present observations, and the reaction NH +H— NH, may be 
added to express the increase in strength of the « band with the decrease in strength 
of the 6 band, though of course the NH, might also be formed in other ways, such 
as N +H, NH). 
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§6. THE HIGH VOLTAGE ARC IN AMMONIA 


In studying the emission spectrum of ammonia in the high voltage arc at 
atmospheric pressure it was necessary to use a constant flow of the gas as the arc 
caused a considerable amount of dissociation into nitrogen and hydrogen. 
Apart from weak N, Second Positive bands and the NH B band (also quite weak) 
the spectrum of the arc consisted of the so-called « band of ammonia. It has 
been found possible to photograph the « band in the region 6000 a.-9000 a. on 
the 10-ft. grating, and it is hoped to communicate measurements of these spectro- 
grams and further details of the work on ammonia at a later date. 
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Rotational Analysis of the (1,0) Band of the 
N, First Positive System 
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ABSTRACT. The (1, 0) band of the Ny First Positive system has been photographed at 
a dispersion of about 1:2 cm~1/mm. ‘The wave numbers and quantum numbers of the 
lines, and the rotational term differences in the upper and lower states are given. The 
spin tripling in the *S state and the A doubling in the *II state are evaluated. 


SiG LONMEROUND ECAP LOIN 

URING the course of other work it was found necessary to have at hand a 
Io rotational analysis of the (1,0) band of the First Positive system of 

nitrogen. This band lies near 8900 a. in the photographic infra-red and 
no measurement or analysis of it appear to have been published. Since the 
photographic infra-red region is becoming of increasing importance both astro- 
physically and in laboratory spectroscopy, and since the nitrogen molecule is very 
important for both these fields of study, it seems advisable to make the present 
analysis available. 

The First Positive system of nitrogen is strongly excited by an uncondensed 
discharge in nitrogen and, much less strongly relative to the Second Positive 
system, in air. ‘The system consists of several well marked sequences of bands 
extending from about 50004. to beyond 10,0004. ‘The (5,2) band at 6705 a. 
and the (6,3) band at 6623 a. have been rotationally analysed by Naudé (1932). 
His work showed that the transition is 9I[>* and that the system has its initial 
level in common with the final level of the Second Positive system. The 27 
branches of each band (of which 9 are main branches) are arranged in three 
sub-bands. In the present paper wave numbers of the band lines are given 
together with rotational term differences in the upper and lower states. In the 
3X state the spin tripling has been evaluated for the level v’’=0 and in the II 
state the A doubling has been evaluated for the three sub-levels with v’=1. 

The (1, 0) band was photographed in the first order of a 20-ft. Eagle Mounting 
concave grating with 30,000 lines/inch. The dispersion was 1a/mm. which 
corresponds to about 1-2cm™!/mm. in this region of the spectrum. The source 
used was a simple end-on discharge tube run off a transformer at 3,990 Vv. and 
0-lamp. ‘The pressure of the nitrogen (which passed continuously through the 
tube) was adjusted to give maximum brightness of the bands. Kodak 
I.R.E.R. plates, hypersensitized in ammonia, were used. 


§2. ANALYSIS OF THE (1,0) BAND 
The lines of the (1,0) band were measured against iron arc lines (Meggers 
1935). Due to the small number of these lines and to their considerable breadth 
the wavelength measurements are somewhat less accurate than the instrument is 
capable of giving. Wave numbers of the lines are given here to two decimal 
places, and for unblended lines they are probably good to +0:01cm-1, 
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In Table 1 are listed the wave numbers of the lines and their estimated 
intensities (visual scale 0-10). No lines of the following satellite branches were 
ROWS MER, oP eo Ria Ova 

It will be seen from Table 1 that the lines of low K value are often missing. 
In the (1, 0) band this is very much more noticeable than in the (2, 1) band which 
has also been studied and it may be due to self-absorption for the metastable 
lower level with v'’=0. The lines in any one of the branches show alternating 
intensities in the ratio of about 1:2 as predicted for a molecule such as N, with 
a nuclear spin of unity. 

(1) Lower State, =X 


The spin fine structure separations in the ®& state are, given by: 
Afey(K) ="Rya(J) — QV +1) ="Qp(J) — PJ +1) 
=R(J) —*Qn(J +1) =Q.) —°Pa(J +1) 
=*Rgo(J) —°Qai(J +1) =" Q3(J) — ®Paa(J +1) 
Afes(K) =" Quo) —*Ras(J — 1) =°Pra(J) —°Qua(J — 1) 
= Q.(J) — ®RST Dara) =? Qos(J =) 
="Qso(J) —Ra(J — 1) =°Pso(F) — Qa — 1). 

In Table 2 the values of these differences that can be obtained from Table 1 
are given. ‘The triplet splittings are close to those found by Naudé for the levels 
v’=2and3. ‘This should be so since the splitting varies slowly with vibrational 
quantum number. ‘The values are however slightly higher than those found 
by Naudeé. 

In Table 3 the rotational term differences A,F’’(K) are listed. These 


differences can be obtained from several combinations of branches. In Table 3 
they have been computed from the main branches as follows: 


eal gh y= Real) = Pee )), 
KeaJ ME Ve RI 1) Ps 1); 
Ree eer Ce Rh=1)— Pa 1). 


(ii) Upper State #11 
The rotational term differences in the upper state 
ALF’ (J) =R(J)—P(J) t= 1, 2 or 3, 
are not given. They agree with those obtained by Hulthén and Johansson (1924) 
from analysis of the N, Second Positive bands as should be the case. 

To obtain data on the A-type doubling in the *II state it is necessary first to 
obtain the rotational term differences A, F’,,(J) and A, F’,,(J). These are found 
by means of the following relationships: 

AyFroal J) =R(J)-Q), 
Table 4 gives these differences. 
The mean of the A doublings (8,,(J +4)) for the levels J=K and J=K +1 
is found from the relationships: 
| AF .dJ) —AiFacdlJ +1) =Avnd +1) -Av) 
=28,q(J +1). 
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Table 5 gives the A doublings for the three substates of the 31 state. The sign 
of the values for the II, level is reversed in Table 5 since in the theoretical treatment 
(Hebb 1936) components which show the same symmetry for a given value of J 
are grouped together rather than those components with the same designation 
c or d. 

The A doublings in the v’ =1 level show the same general trends as found by 
Naudé in the v’ =5 and 6 levels. That is: 

(1) 5, large at low J values and decreasing as J increases, 

(2) 5, and 8, small at low J values and increasing as J increases, 

(3) 8) >5,>53. 

The actual values are, as would be expected, slightly different from Naudeé’s. 


Table 1 
De if Q, I R, iT FO; I K 
Tit 7-47 1 
11217-20 0:5 11223-96 0-5 2 
11217-47 3 11227-20 2°5 3 
11217-78 0 11230-69 2-5 eld 25deot 4 4 
11218-23 6 11234-18 3 11257-52 1:5 5 
11218-90 1 11237-84 3 11263:91? 1°5 6 
11219-56 3 11241-71 5 11270-80 5 7 
11220-49 f'5 11245-78 3 11221-71 0:5 8 
11221-58 B52 1250-4 eee aS 11285-16 5-5 © 1222-78 3°5 9 
11222-78 3-5 11254-57 Ae5 alt 29) coe bumee 11224-09 1 10 
11224-39 4 11259-20 5 11300-40 5°5  11225-61 3°5 11 
11226-13 2 11264-11 5 11308-37 4 11227-35 1 12 
11228-10 4:5  11269-29 6 11316-60b 6:5 11229-33 4 13 
11230-33 2-5 11274-73 Deo abloo aa 5-5 11231-56 0:5 14 
11232-81 4-5  11280-39 7 11333-89 6 11234-03 1 15 
11235-56 2 11286-33 5°55 11342-92 5 11236-79 0 16. 
11238:56? 4 11292-53 af 11352-25 a 11239-81 1:5 17 
11241-89 1:5 11299-03 5-5 ett GGl-o7 5 11243-19 1:5 18 
11245-49 5 11305-83 7 11371-77 Zi 19 
11249-44 2 11312-88 6 11381-95 5 20 
11253-62 4 11320-26 6:5 11392-44 6 21 
11258-09 D 11327-91 5 4 22; 
11262-81 4 11335-87 i 11414-31 6 23 
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Table 1 (cont.) 


RQa I Pos I @Ps, f Regs I ss 
11372277 11199-81? 0-5 11312-68 0 11312-34 2:5 18 | 
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Spin Triplings in the 3% state 
Afss(K)  Afe(K) 
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Rotational Term Differences in 3 state 


KO ALF, "(K) Ask 2"(K) 
Ni MORIA) 95-09 
£7 101-03 100-94 
18 106-76 106-69 
19 112-43 112-46 
20 = 118-15 118-19 
21 123-86 (23°91 
220 129°63 
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24 141-02 141-02 
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Table 4. Rotational Differences A, F’j.q(K) and A,F"jq-(K) in #1 state 


311, 311, STI, 

K AFialK) AF ’iead(K) ALF ’2acK) AF ’2ca(K) AF ’gac(K) ALF ’3ca( FX) 
2 4-66 
g 8-68 
4 12-10 13-95 
5 15-74 21-05 14-14 L749 
6 18-94 26:07 19-30 20:94 
7 22:16 29-08 22-66 25-23 20°88 24-28 
8 BIDS 25°85 28°41 24°29 27-66. 
9 28°53 SOS DDS: 3165 27°65 31-03 
10 Sile7) 38-20 S55 34:74 31-30 34-11 
11 34°81 41-20 SOO: SUL37 34:33 37°72. 
12 37:98 44-26 39-08 aaO2 SEO: 41-04 

13 41-19 47-31 APS 44-21 40-99 44-32 

14 44-40 50-39 45-50 47-39 44-18 47-63. 

15 47°58 53:50 48-80 50:57 47-61 50:93: 
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17 Soy S12 55:21 56°95 53-96 57:46 
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21 66:64 72:18 68-28 66°81 70°45, 

Jap 69°82 TIES O YORI 69-99 73°69 

23 73-06 78-44 Weevil 76:02 F332 

24 76°30 81:57 77°89 

25 79-42 81-09 

26 84-27 

27 87:47 

Table 5. A Doublings in II state, d,4(J +1) 

K sTI, 311, STI, K 311, sIT, eh 
4 — Oe Or 15 ilo3}7/ 0:74 0-11 
5 16 Ses 0-72 0-12 
6 1-96 0:03 i 29 0-80 0:15 
y 1-90 0-31 =—0-01 18 125 0-83 0-14 
8 0-41 0-01 19 1-25 0:88 0:20 
9 1:63 0-52 —0-14 20 ikea? 0:91 0-21 

10 1-70 0-55 etl 21 1:18 0:23 

11 1-61 0-61 0-01 22 0:95 0-19 

1? 1-54 0-67 0-03 23 1:07 0:94 

13 1-86 0:65 0:07 24 1-08 

14 1-41 0-71 0-01 
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An Absorption Comparison of the 8-Particle Spectra of 787AcC’ 
(allowed), 7;3;RaE (second forbidden) and 3:5 yr.-?°T1 
(third forbidden) 
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ABSTRACT. Experiments are described in which a direct comparison is made by an 
absorption method of the f-particle spectra of actinium C’, radium E and 3°5 yr. thallium 
representing allowed, second forbidden and third forbidden transitions respectively. 

New values of 623+4mg/cm*® and 300+3 mg/cm? are proposed for the absorption 
limits in aluminium of the 8 particles emitted by actinium C” and thallium and it is shown 
that, of the three elements, radium E emits most and actinium C” least low energy electrons. 

An attempt has been made to derive energy distributions from the absorption curves, but 
the results for radium E depart markedly from those obtained with magnetic spectrographs. 
The general findings, however, agree with those of other workers who have attempted the 
same procedure. 

The formula used in the derivation of the energy distribution includes an approximation 
which necessarily leads to an apparent deficiency of low energy particles, but the accuracy of 
the data is not sufficient to permit the use of a better approximation. Strictly for purposes of 
comparison, however, the method has obvious applications. 


§1. INTRODUCTION 
HERE is still very little experimental information on the differences in. 
shape amongst the f-particle spectra of bodies undergoing simple transitions 
of which the degree of forbiddenness is known. Most of the investigations 
so far have been carried out by means of magnetic spectrographs (cf. Martin and 
Richardson 1948), but there is no record of an absorption investigation of the 
problem, probably because of the difficulty of deriving an energy distribution 
from an absorption curve. It is desirable, however, that the absorption method 
should be extended to this work since it allows investigation of radio-elements 
from which sources cannot yet be made sufficiently strong for spectrographic 
investigation and it simplifies the apparatus required to a considerable degree. 
The first serious attempt to solve the problem of deriving an energy distri-- 
bution from an absorption curve is that described by Das Gupta and Chaudhury 
(1948). These workers compare their experimental results for radium E with 
the theoretical predictions of the Fermi and so-called Konopinski-Uhlenbeck 
theories and conclude that radium E f-particles follow closely a K.U. form of 
spectrum. From the results given in their paper, however, it would appear that 
their absorption measurements are not sufficiently accurate to warrant the analysis 
to which they are subjected. In addition since the Fermi distribution applies 
strictly only to allowed transitions it is not to be expected that radium E (second 
forbidden) should exemplify it—and the K.U. theory has effectively been discarded. 
in recent developments. 
A similar method has been described by Hughes, Eggler and Huddleston 
(1949). These authors compare the energy distributions obtained from the: 
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absorption curves with those obtained by spectrographic methods. ‘This is a 
more realistic approach to the problem in its present state than that adopted by 
Das Gupta and Chaudhury. 

The present paper describes absorption measurements on three B-particle 
emitters, actinium C’, radium E and thallium (3-5 yr.) known to have effectively 
simple B-spectra and to emit little or no y-radiation. ‘These bodies were chosen 
for investigation as providing examples of allowed, second forbidden and third 
forbidden transitions respectively, their comparison being the more profitable 
because their atomic numbers are so nearly the same. As a result of the measure- 
ments new values are proposed for the absorption limits of the B-particles from 
actinium C” and 3-5yr. thallium. Range and energy distributions are derived 
for the B-particles from the three bodies and these distributions are directly 
compared. . 


S27 APPARATUS 


The electronic apparatus was all of A.E.R.E. design and consisted of a power 
unit, type 200, a scaling unit, type 200A, a timing unit, type 1003K, and a probe 
unit, type 1014A. The Geiger-Miiller counter was a G.E.C. type G.M.4, with 
a 6mg/cm? duralumin end window. ‘The counter, sources and absorbers were 


" 
es 


Figure 2. 


Enlarged Section showing 
Clamping of Counter 


Figure 1. 
A. G-M counter. C. Space for absorbers. 
B. Source holder. D. Clamping plate. 


mounted in a Perspex assembly (Figure 1); the source holders (Figure 2) were also 
Perspex and were replaceable in a fixed position and orientation in the assembly 
shown in Figure 1. Perspex was used in an attempt to reduce back-scattering. 
The active material was contained in the recess in the top of the source holder. 

The short half-value period of actinium C” (4-76 min.) made it necessary to 
adopt a method for rapid introduction of fresh sources of this material. A strip 
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of Perspex with a standard source holder suitably built up at each end was rotatable 
in a horizontal plane about a brass rod supported in an open-sided Perspex box. 
The source holders projected through 1 in. square holes in the top of the box by 
the same amount as the single holders in the assembly shown in Figure 1. Over 
one of the holes was fitted the counter assembly with the bottom square plate 
removed and over the other was situated a brass pot into which a brass button 
activated with actinium B could be inserted. The recesses in the two source 
holders were painted with a suspension of colloidal graphite in acetone and were 
connected beneath the Perspex strip by a fine wire which was earthed; the 
actinium B source was connected to a positive 2,000 volt supply. By depressing 
the brass rod against a light spring and rotating it through 180°, an actinium C’ 
source could be changed over and counted, two or three seconds after activation 
was completed, and another source placed in the activator at the same time. 
Moreover—and this is very important—the conditions prevailing in the experi- 
ments with thallium and radium E were reproduced precisely throughout. 


§3. PREPARATION OF SOURCES 


The thallium was in the form of thallium oxide of high specific activity activated 
in the pile at Chalk River, Canada. A small amount of the material was spread 
uniformly in the recess in the top of a source holder and was cemented in with a 
few drops of chloroform which dissolved the top layer of Perspex. Two sources 
were used of approximate strengths 1-6 microcuries and 0-016 microcuries and 
thicknesses estimated at less than 1mg/cm? and much less than 0-1 mg/cm? 
respectively. There was thus practically no self-absorption in the sources. 

For radium E, the most satisfactory method of preparation was found to be 
as follows: a small piece of silver foil was immersed in a solution of 5 microcuries 
of radium D for 48 hours to remove polonium. A piece of copper foil was then 
left in the solution for 24 hours collecting radium E; this foil was dissolved in 
nitric acid and the radium E was precipitated out, with neodymium carrier, as 
the hydroxide. After several re-precipitations to remove all traces of copper, 
the active material was washed and dissolved in hydrochloric acid; a portion of 
the solution was transferred to a source holder, evaporated under an infra-red lamp 
and cemented in with chloroform. ‘The thickness of this source was estimated 
to be less than 1 mg/cm? and its strength was about 0-5 microcuries immediately 
after preparation. 

The primary source for the preparation of actinium C” was in the form of 
actinium fluoride with lanthanum fluoride carrier. Actintum B was collected 
from this on the face of a }in. diameter brass button by the normal recoil method, 
using a potential of 230 volts. This activated button was then used for the 
collection of actinium C’ by recoil, in the source changing apparatus described 
earlier. 

§4. THE ABSORPTION CURVES 


The experimental procedure followed the normal routine of observing the 
average counting rate for a number of absorber thicknesses, using aluminium 
absorbers and applying the usual corrections. seat 

The absorption curve for the 3-5 yr. thallium, shown in Figure 3, curve (a), was 
built up in two sections which were fitted together by taking the two counting 
rates for an absorber thickness in the region where they overlapped, assigning to 
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them the common value unity and expressing all other readings for both sources 
as fractions or multiples of this counting rate. All points had standard errors 
between 1-24% and 0-80%. 

A single radium E source was used for the whole absorption curve for this 
body (Figure 3, curve (6)), the readings for small absorber thicknesses being made 
16 days after the first reading at maximum absorber thickness. The points on 
the straight portion of the curve beyond the f-ray end-point have standard errors 
not greater than 6-3°% and all points for absorber thicknesses less than 331-6 
mg/cm? have standard errors of about 1%; no point has a lower error than 
086%. Error limits are shown only for some points in the straight lower portion 
of the curve. In the figure the actual recorded counting rate is plotted, after 
correction for decay. 


Counts/Second 


J 
500 600 700 800 300 1000 
Al (mg/cm?) 


Figure 3. The absorption curves. 
The background counting rate has been subtracted in all cases. 


Measurements of the half-value period of actinium C” indicated that it was 
being collected radioactively pure. A ‘control’ count was taken for each freshly 
activated actinium C” source, with no absorber in place, and the counts taken 
with absorbers, using the same source and after applying the usual corrections, 
were expressed as fractions of the ‘control’ count. Each point on the absorption 
curve, (c) of Figure 3, is the mean of at least 5 readings with different sources and 
the same absorber thickness. The intensity of the ‘control’ count, corresponding 
‘to the counter window and air gap thickness of 7-7 mg/cm?, is taken as 10,000 units 
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for the purpose of the figure and all other results are expressed relative to this 
reading. ‘Though the statistical accuracy of the points decreases towards the 
lower end of the curve, errors are in general too small to insert without confusing 
the graph. Three typical points have the following standard errors on the 
vertical scale: 46-2mg/cm?, 0:56%; 161-7mg/cm2, 0-92%; 324-9 mg/cm?, 
2-4°%,; on the part of the curve beyond the f-ray end-point the errors are much 
greater since the total counting rate was very little in excess of background. 
The four points shown with their errors were determined to within 30°, and were 
used exclusively for fixing the straight line portion of the curve. With the material 
available it was not considered worth while to attempt to increase further the 
accuracy of these points since to reduce the standard error to 5°/, would necessitate 
counting for at least 180 hours and the use of more than 1,000 actinium C” sources. 
If a more detailed examination of this part of the curve is required a much stronger 
primary source is desirable. 


§5. ANALYSIS OF THE RESULTS 

The first analytical procedure was to obtain curves for the pure f-activity of 
the three elements. This was done by extrapolation back to zero absorber 
thickness of the straight end portions of the curves of Figure 3. Assuming this 
construction to give the y-ray and bremsstrahlung effect, this was subtracted from 
the total intensity of both 6- and y-rays at convenient intervals on the horizontal 
scale and the resulting f-intensity was replotted on a logarithmic scale, with all 
three curves passing through the same point on the intensity axis, this point being 
taken as unity. The ranges in aluminium of the f-particles from 3-5 yr. thallium 
and actinium C” were then calculated, using the method described by Feather 
(1938); radium E was taken as the reference body on account of the well authenti- 
cated value of 476 mg/cm? for the absorption limit in aluminium of its -particles. 

This analysis yielded 623 + 4 mg/cm? for the absorption limit in aluminium of 
the primary f-particles of actinium C” and 300 + 3 mg/cm? for the corresponding 
limit for the B-particles of 3-5 yr. thallium. The shapes of the ‘calculated range’ 
curves indicated that actinium C” and 3-5yr. thallium both emit relatively fewer 
low energy electrons than does radium E, the deficiency being much greater with 
actinium C” than with 3-5yr. thallium. The former result was to be expected 
since actinium C” undergoes an allowed transition, but the fact that 3-5 yr. thallium 
(third forbidden) definitely emits fewer low energy particles than does radium E 
(second forbidden) is new knowledge. 

The results were also treated according to the method of Sargent (1939). 
This method is based on the conclusion, derived from a limited set of experiments, 
that beyond an absorber thickness equal to about half the effective range R the 
form of the absorption curve (logarithm of f-intensity against fraction of effective 
range) is the same for all simple spectra. The results shown in Figure 4 do not 
bear out this conclusion in detail. In constructing Figure 4 effective ranges found 
by Feather’s method are used and the three curves are made to be concurrent 
(at unit intensity) at a point corresponding to two-thirds of the range in each case. 
Whilst it will be noted that the curves remain indistinguishable for absorber 
thicknesses greater than 2R/3, the curve for actinium C” begins to separate from 
the other two at about 0-65R though the curves for radium E and 3-5 yr. thallium 
do not diverge appreciably until roughly 0-45R (in agreement with the Sargent 
conclusion). As was to be expected the Sargent method of comparison gives the 
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same result as the less direct ‘calculated range curve’ method of Feather regarding 
the relative proportions of low energy electrons; quite obviously from Figure 4 
actinium C” emits least and radium E most low energy electrons amongst the 
three bodies investigated. 

Very little work has been described in the literature on more detailed methods 
of analysis of absorption curves. It should, however, be possible to derive range 
or energy distribution curves from them. This has been attempted by Das Gupta 
and Chaudhury (1948) and by Hughes, Eggler and Huddleston (1949) ; inde- 
pendently a similar method was applied to the results above described. In this 
connection a slight modification in procedure may be claimed as an improvement 
in the general method of analysis. 


Intensity 
S 


0 0-2 0:4 0-6 0-8 
Range Fraction 


Figure 4. Sargent plot. 


Let J(R) represent the distribution with respect to range of primary f-particles 
emitted by a source ; then the fraction of f-particles having ranges between R 
and R+dR is f(R)dR. ‘The fraction I(x) of these particles recorded through an 
absorber of thickness x may be written 


Iay=d +a(F) +65) +5) cheesy 


a, b, c.... being constants (Rutherford, Chadwick and Ellis 1930). Experiments 
with homogeneous f-particles show that to a reasonable first approximation we 
may take I(x)=(1—x/R). Then if, for a source having a continuous spectrum 


of primary B-particles as specified above, A(«) is the fraction recorded through a 
thickness x 


A(s)={ (2 = ) {Duk (1) 


whence (x) =xA"(a), 
dashes denoting differentiation with respect to «. 
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Saleen (2) clearly provides a method of deducing the range distribution 
Ww e the intensity A(x) is plotted on a linear scale (as in the methods of Das Gupta 
a hide aes mn of eee et al.). It is customary, however, to plot A(x) 
ogaritimic scale, as has been done in the present work. In this c:; i 
L(x) =1n A(x) it can be shown that twa s 


F(x) =«A(x)[L"(~)+{L'(w)] (3) 
Also if 4(£) is the distribution function with respect to energy, since 
$(E£) dE =f(x) dx, 


it 1s only necessary to know the appropriate range—-energy relation in order to 


be able to derive the distribution function relative to energy from the absorption: 
curve. 


x10> 
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Figure 5. Graphs of L’(x) against absorber thickness. 


The numerical values only of L’(x) are plotted; L’(x) must be prefixed by a negative sign for 
insertion in the formula. 


Comparing equations (2) and (3) it will be noted that the former involves the 
second differential of the experimentally determined A(x) whilst the latter 
involves both the first and the second differentials of L(x). Formally, this. 
difference represents the improvement in procedure already mentioned: the 
determination of the second differential of an experimental curve is always attended 
with very considerable uncertainty, unless the accuracy of the original points is 
exceptionally high and unless very many points have been taken. An analysis 
based upon (3) would appear to be preferable to one based upon (2) since F(x) 
determined from (3) is dependent more upon the value of L’(x) than upon that of 
L"(«). In applying (3) to the experimental absorption curves of the present 
investigation it was found that, except at the low energy ends of the curves, L(x), 
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was generally only of the order of 10°% to 20% of {L’(x)}*. Thus it seems that 
there is good reason for preferring such an analysis to the type employed by 
Das Gupta and Chaudhury, and by Hughes e¢ al., which is based on (2). 

The graphs of the first differentials of the B-particle absorption curves are 
shown in Figure 5. For convenience both Figure 3 and Figure 5 have been 
plotted using logarithms to the base 10; conversion to the base e was made later 
in the calculation of F(x). 

The calculated distributions with respect to range for the three 8 spectra are 
shown in Figure 6. The scale of abscissae gives the fraction of the maximum 
effective range of the B-particles from each body, the area under each curve being 
maintained at unity. 
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Figure 6. Range distribution curves. 


_ The distribution with respect to energy was calculated with the help of data 
given by Marshall and Ward (1937) and Curie (1935) from which a range-energy 
curve was constructed. ‘The three energy distributions so obtained are shown 
in Figure 7. In Figure 8 the calculated energy distribution for radium E is 
compared with the currently accepted distributions obtained using a magnetic 
spectrograph (Flammersfeld 1939, Neary 1940). It was not possible to compare 
similarly the distributions for 3-5 yr. thallium and actinium C”, there being no 
reference to magnetic analysis of these spectra in the literature. 
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§6. DISCUSSION 

Here the general results of the absorption experiments will be referred to 
briefly before the calculated distribution curves are discussed in detail. Using 
the range-energy relation of Feather (193 8) the maximum energy of the B-particles 
of actinium C” is obtained from the absorption limit of 623 +4mg/cm? as 
1-442 +0-008 Mev. This value is to be compared with 1-47 mev. deduced by 
Sargent (by the alternative method already mentioned) from an estimated 
absorption limit of 0-68gm/cm?. The measure of agreement between these 
values differently determined is not unsatisfactory. For 3-5yr. thallium the 
absorption limit of 300 +3 mg/cm? deduced in the present work is sufficiently 
close to the less accurate 0-31 gm/cm? originally reported by Fajans and Voigt 


(1941) to support the conclusion that the maximum energy of the f-particles of 
this body is 0-85 + 0-01 mev. 


x10? 
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Figure 7. Energy distribution curves. Figure 8. Comparison of the derived 


energy distribution for RaE with 
magnetic spectrograph determina- 
tions. 


The other point of some interest in the results concerns the magnitude of the 
(y-ray + bremsstrahlung) effect with the three bodies. From the extrapolations 
in Figure 3 it will be seen that the values of the ‘y/f ratio’ at zero absorber thickness 
are (3-6 + 0-6) x 10-4, (0:53 + 0-01) x 10-4 and (1:38 + 0-02) <Osetor. actinium C ; 
radium E and 3-5 yr. thallium, respectively. Assuming that the value for radium 
E represents a pure bremsstrahlung effect, and taking count of the difference in 
average f-particle energy as between the three spectra, 1t may be concluded that 
the true y-ray effect with 3-5 yr. thallium is certainly no greater than that with 
actinium C”. In the latter case it is known that a y-ray of 870 kev. energy is 
emitted about once in 200 disintegrations (Surugue 1941). In any case the 
assumption that the present work has been concerned with essentially simple 


f-spectra is fully substantiated. 
39-2 
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Turning now to the calculated distribution curves, it will be noted that pate 
they give approximately the correct end-point energies for the three materials, 
the maximum of the radium E curve certainly falls at too high an energy. It iS 
reasonable to suppose that this is also the case for 3-5 yr. thallium and actinium C’”. 
Such an effect is to be expected since the initial absorption of homogeneous 
B-particles is slower than that given by the linear relation I(x)=1—x/R. 
Experimental curves obtained by Marshall and Ward (193 7) are of the form shown 
in Figure 9 while the linear relation used is given by the line AB. The line AC 


Intensity 


mg/cm? 


Figure 9. Absorption of homogeneous f-particles (Marshall and Ward). 


might be thought to give a better approximation to the whole curve, and to use it 
reduces slightly the energy corresponding to the maximum of the distribution 
but it also reduces greatly the end-point energy, so that no overall improvement is. 
derived. Clearly, the true state of affairs can only be approached by considering 
more terms in the expression 


(x)=1+0(5) +6(3) so ites 


for the absorption of homogeneous f-particles; this, however, introduces. 
differential coefficients of higher order than the second and these cannot be obtained 
with any reliability from the absorption curves. Hughes et al. (1949) have 
obtained very similar results with ®Cu, ?”8Au, 18°W and radium E comparing the 
energy distributions derived from absorption curves with magnetic spectrograph 
data. ‘They found that the end-point energy was correct in all cases, but the 
maxima of the calculated curves occurred about 100 key. higher than was expected. 
To overcome this difficulty, they proposed the empirical method of shifting the 
calculated distribution by 100 kev. towards lower energies, for an end-point 
energy of 600kev., so that the maximum would fall at the correct energy. An: 
obvious objection to this procedure is that the end-point energy is then low by the 
amount of the shift. The distribution obtained (Figure 7) for radium E agrees. 
fairly well with these findings since the maximum of the curve is about 200 kev. 
too high; this is double the value reported by Hughes et al., but the end-point 
energy is also approximately double that for which they propose a shift of 100 kev. 
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Excited Electronic Levels in Conjugated Molecules— 
V: A Valence Bond Estimation of Energy Levels in Aromatic 
Hydrocarbon Molecules * 


BraG. Rk BAL DOCK 
King’s College, London 


Communicated by C. A. Coulson; MS. received 28th November 1949 


ABSTRACT. This note is concerned with the calculation of the energy of the ground 
state and the lowest excited state for hydrocarbons of the condensed ring type. The differ- 
ence between these two energies is the energy of the first electronic transition. ‘The normal 
valence bond method is used, taking only the Kekulé structures into consideration. ‘The 
results are compared with the experimental absorption wavelengths, and are in good 
agreement, giving a value of the exchange integral « between —1-8 and —2 electron volts. 


SiN ROC hON 
T is well known that the simple molecular orbital treatment yields a good 
picture of ground state phenomena in condensed ring moli-cules, but in calcula- 
tions on the excited states considerable difficulties arise. At the very least it 
becomes necessary to use antisymmetrized wave functions (Jacobs 1949), and at 
present this method is too cumbersome to be applied to large molecules. A 
similar complexity occurs if the valence bond, or resonance, method is used, 
unless nearly all the excited structures are neglected. This approximation has 
therefore been made in all our present calculations, where Kekulé structures 
alone have been taken into consideration. We used the standard technique, in 
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Table 1 
Compound E,~2 Bene Eels 
a a fo 
Hexacene 6-04 5:14 0-90 
Heptacene 6:60 5-84 0-76 
Benznaphthacene 6:13 4-69 1-44 
Benzpentacene 6°88 5-66 jen 
Dibenzanthracene 6:60 4-82 1-78 
Dibenznaphthacene 7-40 5-84 lesz 
Dibenzpentacene 8-10 6-78 1:32 
Picene 6:70 4-91 1E79) 
Benzpicene egy) 6:05 1-74 
Benzpyrene 5-68 4-03 1-65 
3489—Dibenzpyrene 6°57 Siedul 1-46 
Anthanthrene a) 4-27 1-48 
Coronene 6:64 S07 Wosy7/) 


* Cook, Schoental and Scott (private communication). 


+ Patterson (1942). 
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* Calculated by Forster (1938). 
t+ Calculated by Seel (1947). 


which the terms of the secular determinant are calculated by means of Pauling’s. 
method of superposition diagrams (Pauling 1933). Although this approximation 
is known to give false values of the ground state energy (Jonsson 1942, Pullman 
and Pullman 1946), it has been shown by Forster (1938) and Seel (1947) to give 
in general a good estimate of the energy of the first transition. 


§2. RESULTS 


Table 1 shows the results of the calculations, expressed in terms of the exchange 
integral «. Q is the coulombic integral, Z, the ground state energy, E, the energy 
of the lowest excited state, and FE, — E, the energy of transition. The last column 
shows the wave number of the longest wavelength absorption maximum, quoted 
from Cook, Schoental and Scott (1950) unless otherwise indicated. 

In the Figure we have plotted the experimental energy of transition against 
the calculated values of (E,—E,)/x. A few previous results for other molecules 


are also included. 
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S2h IDUSCO SWORN! 

The conclusions from this work are as follows: 

(a) The fact that in the Figure all the points lie approximately on a straight 
line shows that theory and experiment are in essential agreement. This agree- 
ment is unexpected in view of the crudity of the approximation adopted in the 
calculations. ; 

(b) The results shown in the diagram suggest a value of « in the region of 
—1-9ey., which compares well with previous estimates listed in ‘Table 2. 


Tables 
Author Method —«a (ev.) 
Pauling and Wheland (1933) Thermochemical data 155 
Sklar (1937) Fluorescence of benzene 1-92 
Lloyd and Penney (1939) Vibration spectra 185 


(c) The excited structures are known to contribute to ground state pheno- 
mena, such as the resonance energy, increasingly with the size of the molecule. 
In our calculations these structures have been entirely neglected. Therefore 
deviations from the proportionality of the calculated and observed energies of 
transitions might be expected to increase with the size of the molecule. In fact 
they do not increase with the size, as can be seen from the diagram, since a single 
value of « correctly describes a series of molecules varying in size. It is the shape, 
or ‘angularity’ (see Seel 1947), which is an important factor, since the points 
corresponding to molecules with different shapes lie on distinct straight lines. 
. Thus the points representing the polyacenes lie approximately on the straight line 
A, the molecules 7, 8, 9, 10 give the line B, and the molecules 12, 13, 14 give a 
third line C, as is shown in the Figure. These three lines (which all pass through 
the origin) are represented by three distinct values of «, —1-98, —1-86 and 
—1-78 ev. respectively. ‘The more angular the molecule, the smaller the numerical 
value of «. (Highly condensed molecules, such as pyrene and coronene, are 
excluded from this discussion.) 

(d) The numerals on the right of the figure indicate the minimum number g of 
non-benzenoid rings in the Kekulé formula for the molecule, i.e. the number of 
non-benzenoid rings when the structural formula is written according to the 
Fries rule. It has been remarked by Cook, Schoental and Scott (1950) that there 
is a correlation between q and the longest wavelength absorption maximum ». 
The figure shows that q certainly increases as the experimental value of v decreases, 
and that the points representing the molecules may be grouped between hori- 
zontal lines as shown. ‘There are two exceptions, picene (15) and benzpicene 
(16), which lie in the region q=1, although in their structural formulae q=0. 
The calculated values of v also show this effect but not so clearly as the experimental 
results, for vertical lines cannot be drawn in the figure to separate the points into 
groups with anything like the distinctness of the horizontal grouping. 

(e) The symmetry of the polyacene molecules is such that all the Kekulé 
structures represent states which are symmetrical about the axis parallel to the 
length of the molecule. Thus the method we have used assumes that the transi- 


tions will be polarized along the length of the molecule. This is in direct 
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‘contradiction with the simplest molecular orbital interpretation. But in view 
of the present uncertainty in the exact interpretation of the long-wave ultra-violet 
bands, no clear deduction can be drawn (Jacobs 1949). However, if the 
transitions we have discussed are not polarized along the length of the molecule, 
the agreement with experiment becomes all the more difficult to explain. 


§4. POLYACENES 


The polyacene series can be treated generally by the method we have used. 
If there are 7 rings in a polyacene molecule, there will be 47 +2 nuclei, n +1 Kekulé 
structures, 27 +1 7-bonds in each structure, and 5” +1 non-vanishing exchange 
integrals between neighbouring atoms, all supposed equal in magnitude. 

The energy Ey of each separate structure is given by 


wes +3(2n +1)—4(5n +1) =0, 


EEE Ta Per 
a 2, 


i.e. 


—E 


lo 


Let 2 


+5 el =x} 


where F is the energy of the electrons. ‘Then the resonance energy is E, — E=x«. 
The secular determinant, of order +1, is conveniently written in terms of 
= 2K] 9. 
The following equation for y is obtained: 


y yt yt2 yt+3 CGO Otechtin Ove y+n 

y+1 4 Gr Oem He) 4%(y +-n—1) 

yal syst 1) ay oO) Aaa aaa oe 44(y +n —2) 

y+3 42(y +2) 44(y +1) ABy ver eeeer ees 45(y +n —3) =): 
yen AXytn—1) My tn—2) Ay tn—3) 0++ Py 


By subtracting each column from its successor, and performing the same 
operation on the rows, we eliminate the variable y from the off-diagonal terms. 
Next, multiplying certain rows and columns by suitable factors and making the 
substitution z =15/(15y—17), we get 


172z4+15 PB Z 4 fopalie oberta eta qen ete 
4 z+ b Pa ig 
g x 42z+1) 42 Az corse 
Backs B 422 ACEH Doe wee EPs Orsay fo) 
2 » 42x 44 46x41) ---- 
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By similar operations we obtain 


17u—16 4u—17 


0: 0 
15 60 ° e 
4u—17 44 —2 1 
iz Of © 2) 
60 60 4 ® 
0 t u 1 Wie corner O 50) 
0 0 1 u i Naeycais KG 
0 0 0 1 u-:::: 0 
0 0 (0-0 u 


where u=(15z2 +17)/4. 
Putting u=2 cos 0, the equation reduces to 


256sin(n+2)0—256sin(n+1)0 +96 sinnd — 16sin(n—1)6+sin(n—2)0 _ 
sin 6 


A,= 


| for n>3, 

h o— 85x — 96 212 lieost—s 
yes een oe Seo 1 
A,, is a polynominal in cos@ of degree n +1, with n +1 real zeros since it is the 
determinant of a Hermitian matrix. 

Now A, =4,,B,,, where 
A, =[16 cos (4n +1)6 —8 cos $n + cos (4n — 1)6] sec $0 

and B,,=[16 sin ($n +1)0 —8 sin $n8 + sin ($n — 1)6] cosec $0. 


It can easily be verified that 4, and B,, arise from the secular determinants 
obtained by factorizing the original determinant into two parts involving the 
coefhcients of eigenfunctions respectively symmetrical and antisymmetrical about 
the central plane normal to the length of the molecule. 

From the form of A, and B,, we see that if mis even A, has $n +1 zeros and 
B,, has $n zeros, and for n odd A, and B,, both have $(n +1) zeros. In each case 
all the zeros are real values of 6, so that the m + 1 zeros of A,, are accounted for and 
lie in the range —1<cos@<1. 

Furthermore, for m even or odd, the zeros of A, and B,, are interlaced through- 
out, and each zero @, is restricted within a range 


9 


rm/2—sin- 


ete 
is 2 gees Aces 
n+ 1 


in+1 

Hence the spacing of the energy levels, which all lie in the range 
—2:4<x<2-4, 

becomes infinitely close (being uniformly dense in terms of @) at every point as 

n—> ©. 


The ground state resonance energy thus tends steadily to 2-4% as n> o. 
The ground state energy corresponds to the lowest root 6,, and the lowest 
excited state to the next root ,. 6, is always a root of A,,=0 and @, a root of 
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B,,=0, so that the first transition is always from a symmetrical to an antisym- 


metrical state. It is therefore an allowed transition polarized along the length 
of the molecule. 


a Qn 
For large n, ae apd er eC 
377? 
2 2~ 
so O,°— 6, ~ (n+ 
} 12s Rebs G8 Sl Dif 25 
N C= ~ me 
ow 5 17—8cosd ‘5 ( ra”) eee 
10 1072 
that Hy — Xy = (05% — 01”) = a5 : 
iia Vy MSG (0,2—6, Soren as N—> 00 


Thus we should expect the longest wavelength absorption maximum to. 
increase proportionately to the square of the number of rings, for long strips. 
A similar variation would be expected for all other bands in this series, in agree- 
ment with Coulson’s molecular orbital theory (Coulson 1948). 
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Relation between Bond Structure and the Longest Ultra-violet 
Absorption Band of Polycyclic Aromatic Hydrocarbons 


By J. W. COOK; R. SCHOENTAL AND Es JXe COR 


Chemistry Department, University of Glasgow 
Communicated by C. A. Coulson; MS. received 28th November 1949 


ABSTRACT. A comparison has been made between the values expressed in wave 
‘numbers of the longest ultra-violet absorption bands of a series of polycyclic aromatic 
hydrocarbons and the minimum number of quinonoid rings contained in their molecules. 
It is found that the hydrocarbons fall into series of groups related to this number of quinonoid 
rings. 


N their survey of fluorescence spectra of polycyclic aromatic hydrocarbons 
if Schoental and Scott (1949) noticed a close relation between the position 

of the maximum of the first (shortest) fluorescence band and the minimum 
number of quinonoid rings in the molecules, so as to conform to the Fries rule 
that as many rings as possible should have Kekulé structures. The number of 
-compounds which have been available for the fluorescence studies were, however, 
too small to allow a satisfactory generalization. These authors have found also 


V(cm-') x 1077 


G42 


Figure 1. 


‘that for this series of compounds the shortest fluorescence band coincides 
approximately in wavelength with the longest absorption band. If this is assumed 
to be generally true, then it is possible to examine more extensively the relation 
between the number of quinonoid rings and the wavelength of the maximum 
of the longest ultra-violet absorption band from the numerous ultra-violet 
absorption data recorded in the literature. 

In the Figures and Table are collected the values of the longest ultra-violet or 
~visible absorption bands, expressed in wave numbers, in relation to the molecular 
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formulae of the respective compounds. (Values referring to hydrocarbons 
with highly condensed ring-systems are grouped in Figure 2.) It will be noted 
that values for isomers lie on straight vertical lines. 

The horizontal lines divide the figures into zones in which the compounds 
possess an equal number of quinonoid rings Q regardless of their structure and 
total number of rings. There is some overlap between the zone of compounds 
which have no quinonoid rings (i.e. structures which can be built wholly from 
Kekulé rings) and that corresponding to one quinonoid ring. With the increase 
of the number of quinonoid rings the divisions are more striking. 
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Figure 2. 


There is, however, a class of compounds, namely those having the perylene- 
structure in their molecules, where the spectra consistently fall into zones. 
corresponding to one quinonoid ring more than the minimum number derived 
from the application of the Fries rule. Clar (1932, reference 2) has already drawn 
attention to the fact that perylene in its chemical behaviour does not correspond 
to a simple dinaphthyl-derivative. It may be that other excited structures come 
into play for this class of condensed compounds. Assuming for perylene two. 
quinonoid rings, which would agree with the position of its spectrum, its benzo- 
and naphthologues exhibit the expected shifts in their spectra.* 

The regularities observed by Clar (1941, reference 14) in the ultra-violet 
absorption spectra of series of aromatic hydrocarbons, built up by linear or 
angular condensation, find to some extent their expression in the figures. The 
values of the maxima of the longest ultra-violet absorption bands of each series. 
lie on approximately straight (diagonal) lines. The values corresponding to: 
the linear hydrocarbons lie in the lowest parts of their zones, but that of anthracene: 
is exceptional, in that it is situated in the upper part of its zone. 


* Data for the longest ultra-violet absorption maximum of several polycyclic hydrocarbons 
recorded in the literature have not been included in the table and figures. These need checking 
with the values of their fluorescence spectra; the shortest fluorescence band, being usually the most 
intense, is not likely to be missed, as may be the case with the longest absorption band, of a low- 
extinction coefficient. 
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This simple relationship between the number of quinonoid rings and the 
maximum of the longest ultra-violet absorption (or the shortest fluorescence) 
band of polycyclic aromatic hydrocarbons seems to reflect some fundamental 
connection between the chemical reactivity and the spectra of the respective 
substances, and is worthy of the attention of theoretical chemists. 
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Data refer to the shortest fluorescence band maxima, 


O=number of quinonoid rings. 
* 
+ 
+ 


This structure is assigned by Clar (1948, Chem. Ber., 81, 52) 
to a different hydrocarbon (terrylene), 
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The Scattering of High Energy Charged Particles 
by Thin Foils of Matter 


Bye Sai, BUTLER 
Department of Mathematical Physics, The University, Birmingham 


Communicated by R. E. Peierls; MS. received 3rd November 1949 


ABSTRACT. The distribution of high energy charged particles scattered by a thin foil 
has been investigated and compared with the formula for single scattering. It is shown that 
at sufficiently large angles the distribution approaches that for single scattering and can be 
expressed in the form of a rapidly convergent series of which the first term is identical with 
the single scattering formula. The remaining terms of the series form an accurate expression 
for the deviation from true single scattering. The simplified case of projected scattering has 
also been considered. 


§1. INTRODUCTION 


HEN high energy charged particles are scattered by a thin foil of matter, 
it is well known (see e.g. Williams 1939) that at sufficiently large angles 
of scattering, the angular distribution of the scattered particles 
approaches that for scattering by a single atom, multiplied by the surface density 
of atoms in the foil. Due, however, to the fact that a particle which is scattered 
once through a large angle always receives, in addition, a large number of 
deflections through small angles, the actual distribution is never exactly equal 
to the single scattering distribution. ‘The percentage difference between them 
is less, at any angle, the smaller the thickness of the foil, and it is often of 
importance in scattering experiments to know how thick a foil can be used (and 
thus how great a yield of scattered particles in any one direction can be obtained) 
without the deviations from single scattering at angles of interest being too great. 
It is the purpose here to deal with this question by deriving a formula for the 
angular distribution at angles for which the deviation from single scattering is 
less than about 100%. The corresponding formula for the projected scattering 
on a fixed plane is also obtained (cf. Snyder and Scott 1949). 
In this treatment small angle approximations are employed, and energy 
losses by the particles are neglected. 


§2. GENERAL FORMULAE 


Let a particle be incident normally to a scattering foil, and let f(n,s) dQ be 
the probability that at a distance s through the foil, the particle is moving in the 
direction of the unit vector n within the element of solid angle dQ. Let n make 
an angle 6 with the normal to the foil, and have an azimuthal angle ¢. We measure 
s along the normal. Then f(n, 5) satisfies the transport equation 


2 (n, )=N | 2(0) (ey 9 Ce erie es (1) 


cos 0’ cos @ 


where @ is the angle between n and n’, and is the differential cross section 
per unit solid angle for scattering by a single atom. N is the number of atoms 


per cm? in the foil. 


40-2 
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Divide p into g+h, such that g takes account of most probable deflections 
through small angles, and h gives the probability of large deflections.* The 
exact nature in which this splitting is carried out will be unimportant, and for 
convenience we assume a sudden transition from g to h at some angle 6, (Figure 1). 
Equation (1) now becomes 


He s)=N { o(0) {ay - fo dQ’ +N] h(@) {4 (tise fa dQ’. 


cos 0’ cos cos 0 cos 8 
eeeree: (2) 
For small angles © (those for which g(@)>0), write n’ =n+e and expand 
f(n’,s) in powers of e. Thus 
2 2 
cei +S + shams) 


yur ClO rake aon 
f(n’,s)= f +E ae tIay tH axe * axay 787 aye 


where (x,y) are cartesian coordinates of the projection of n on a plane parallel 
to the foil, and (&,7) the corresponding coordinates for € such that ©? =£7 +7. 


4, 4, 
Figure 1. The division of the cross section p into small angle scattering (¢) and 
large angle scattering (/). 


We now substitute this expansion in the first integral of (2) and, in this integral, 
put also cos @’=cos@. It is evident, on account of the spherical symmetry of g, 
that terms containing odd powers of € and 7 vanish and, moreover, that 


| €29(0) a0" = | 29(0) da’ =4 | 02¢(@) aor 


and | &4g(@) dX’ = | n#g(0) a0’ =3 | @17¢(0) do’ =8 | @*¢(@) da”, 
‘Therefore (2) becomes 
ee / 40) 4 sor fe 103) are (4) 


ds 2cos@ cos 0’ cos @ 8 cos @ 


r 


oy Ne 1 ANE 
inwhich = k= z (2g(@) dQ’ = =| O2e(@) dQ, 


b= 2 | @%(0) a0, 


V? =0?/Ox? +07/dy? and V4=04/dx4 + 204/Ax2dy? +a4/dy4, 


* ‘The author is indebted to Professor Peierls for suggesting this procedure, 
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Although, as will be seen, the only terms of the expansion (3) of importance 
in (4) are those up to second order, we have considered also the fourth order 
terms in order to have some guide in our choice of transition angle 0,. We shall 
choose 6, such that the effect of these terms is small. 

If we make the small angle approximation cos @=1 and sin =0, the vector n 
becomes identical with the vector @ having plane polar coordinates (0,6). In 
this approximation (4) can be written as the well known diffusion equation 


ee Vf =ul0, S)k a ae Nee (5) 
where u(8,s)= N| h(® — 6’) {/0", s) —f(8, 9} dQ! + = Vf. 


The general solution, subject to the condition that it be normalized to unity such 
that for all s 


| f(0,s)dQ=1 
is 
, exp(—6?/2ks) _/*® exp { —(@—8’)?/2R(s—s’) her ee 
GOO) oo S| Tame u(0’, 5’) dQ’ ds’, 
ee (6) 
"The gaussian exp {—(" +y")/2ks} = exp (— 4"/2ks) te et (7) 


2aks 2aks 


is the solution of the equation 
EL a es 
=~ —-V%f=0 
Cs? ee 
and the integral in (6) is a particular solution of (5). 
Both terms of the function u(8,s) can be made small by suitable choice of 6). 
For f(8,s) we thus substitute in u the pure gaussian (7). Then (6) becomes 


exp ( — 6?/2ks) ON len ue 64 
10S) Oaks (1 t at\ is t Bt 


LN | : | | exp {= (0 —6)/2R(s=5')} jg, _ gry XP (KF ARS) ey acy as! 
; 2k(s—s’) 2nks 


where z= Ns | hdQ, and is the probability that a particle, on traversal of the foil, 


undergoes a deflection through an angle greater than 4,. 

The physical interpretation of this solution is clear ; the first term represents the 
multiple scattering due to many collisions through angles smaller than 4; and 
the second term represents the probability of a particle emerging in a direction 
@ after traversing a foil of thickness s and having received one major collision 
(i.e. through >6,) somewhere along its path. : 

Consider the single scattering integral. In general the maximum values 
of the integrand occur for 6” small and for (6 — 9’) small. We expand, therefore, 
h(0’ — 6”) in the neighbourhood of 8: 

Oh (x — x — x)" Oh fe 


" oh Y Mf SS 
h(@’ — 6) =(0) +(x’ —x—x")=° + (9-9 I"), + 9 
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On carrying out the integration, all terms containing odd powers of (x’—.), 
(y’—y), x” and y” vanish, and we get 


ks Res? Rr's(2n — 1)(2n—3)..++1 oonp \ 
No) + pie me Nee emer yre 


Since h is spherically symmetric, this can also be written 


ks (d*h_ 1 dh Peay ete BN NEGRI. Neate 
Ns (0) + 5 (a +595) é = (ae : 5 a5) ae a) i: =) 

It is seen that the first term Nsh(@) gives the well known single-scattering 
probability. The next terms represent the modifications to this probability 
due to the many small additional deflections which a particle receives. 

By substitution for p we shall show that, for suitable adjustment of 6,, the 
multiple scattering portion of (8) differs very little from the gaussian (7). At 
angles for which this gaussian is negligible, f is thus given by (9) which here 
converges rapidly. 

A check on the accuracy of our results can be made by proceeding to the next 
approximation in which (8) is substituted in (6). In this case the probability 
of a particle undergoing two major collisions appears. Although this is not 
carried out here, it can easily be shown that for such angles 6, as we shall choose, 
the results are changed by a negligible amount. 


§3. SUBSTITUTION FOR THE SINGLE SCATTERING LAW 
The single scattering law p(@) has the form 
7(6) = (aa) ( ” =) 1 a KS) sin? 40+ .... ; 
2Mv? e {sin? 46 + (40,)?}? 
where Z’e is the charge of the particle being scattered, and M its rest mass. In 


the small angle approximation, and introducing the total energy E of the particle 
in units of Mc?, p(@) becomes 


222 e* _ 1 
pO) a eee} (8 + 0,2) coe ae lienerencetys (10) 


The effect of the screening of the orbital electrons is taken into account by the 
factor 0) where 6) =A/27a = h/aMc(E? — 1)"?, a being the usual screening constant 
h?/me?Z"'®. Actually, however, the results are very insensitive to the manner 
in which screening is considered and would be affected very little if, instead of 
the above, we had considered the Rutherford scattering law for an unscreened 
nucleus as merely cut off at the angle 6). 

Assuming 6, >, we obtain from (10) 


ZO 04s) RS 0-65 and kis =0°0,7/8 =04/82, 


LLL, e 2 
h * =n Ns 
where 5 Ns ae a eee (el) 
a 6, ZS 728 E2 
and €=log 1-820, =4-06 +t log oe Fo} Se iter ouonortn (12) 


in which p is the density and W the atomic weight of the scattering material. 
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If we choose (ks)"? =$€1? as unit of angle (and change the unit of solid angle 
accordingly) (8) becomes, on inserting (9) and substituting for / 


f6,9)= 7-| exp(-5) {1-24 gza(1-e+ 5) 


zZ, Cope aa 0S 
talitgtmt gt}: ees (13) 


The series representation of the single scattering integral diverges for Gi 3'5, 
but we shall be considering angles greater than this. 

Unless s is extremely small, when our assumption §,>6,) breaks down, it is 
always possible to make both the z and 1/8z€? of (13) small (say<0-1) by adjusting 
6, such that they are approximately equal. For such values it can easily be verified 
that these terms are of negligible importance in (13) and can be neglected. Some 
typical values, for example, of z and 1/8z€? are those given in the Table, 
corresponding to the scattering of high energy electrons by lead. 


s=0-05 s=0-01 s=0-005 
z 0-2 0-1 0:05 0-2 0-1 0-06 0-2 0-1 0:07 
1/82€? 0-023 0-042 0-070 0:032 0:056 0-083 0-034 0:064 0:084 


€ takes only a limited range of values (say 2-5) and it is found from (13) that 
the gaussian always becomes negligible compared with the ‘single scattering’ 
portion for 644-5. In this range the distribution is thus given by 


{ext Sa? 2768 
een aS lie Medi: SE REG Nn, OB ct 14 
f(9,s) salltata ae ) (14) 


where the first term is the true single scattering. 
The percentage derivation from single scattering is thus given by 


OVA S71 es) 
100(5 +7 +e), CRE aca (15) 


Usually only two terms of (15) are required, but for deviations greater than about 
20°, more terms must be used. 

The unit of angle is 5€”?, where 5 and é are given respectively by (1 1) and (12). 
Although € depends on z (which is adjusted so that z~1 /826") it is, however, 
quite insensitive to the exact value of z taken. For example, in the case of high 
energy electrons being scattered by 0-01 cm. of lead, if we choose z=0:1 (and thus 
1/8zé?=0-056) we find €=4-8; if however, we choose z=0-06 (and thus 
1/8zé2=0-083), €=5-0. In fact for all cases, in the evaluation of €, can be taken 
as something like 0-08. 

If the energy of the electrons in the case just considered be taken as 20 MeV., 
then the unit of angle 5é”2=8-7 x 10-2, and from (15) the deviation from single 
scattering is approximately 30% at 30° and 10% at 45°. , 

Thus, taking S€4? as unit of angle, the scattering distrbution per unit solid 
angle approaches single scattering at angles greater than 40-45. In this range 
the distribution is given by (14), and the percentage deviation from single 


scattering by (15). 
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§4. PROJECTED SCATTERING 


In some cases it may be useful to know the formula corresponding to (14) 
for the projected scattering on a fixed plane. This can easily be derived in the 
same way a¢ (14). 

In thé srhall angle approximation the transport equation becomes 


Ac Jae i x Pyros *){f(%’, 8) —f(%, 8)} a", 


where 


Pooi®)= | P(A dy= | p(x? +y*)8)dy eae (16) 
and x is the angle of scattering on the plane. Dividing p,,.;(x) into 2),0;(#) and 
Aryroj(*) at angles +x, and —.x, and proceeding as before, we find 


exp (—x2/2h's) “exp {—(w—a')?/2k'(s—s')} 


f(x, 8) = Oak ye (1 +small terms) +N] | | (27k’(s —s’))¥? 


’ ny EXP ( =. x"? /2k's’) ” PEL 
x hyroi(* —x ) TRS) dx” dx’ ds 5 
where 


ki =N { 29 ro,(x) dx. 


At angles for which the gaussian is negligible, 
(SAP oes Duet Stud ghee 


f(x, s) = Ns {ye sp eta eu Hie 


Now from (16), 
62 1 
Pproi(*) = INs fx? + 0,2}32 ? 
and we find k’s =67&’, where 


| | 


2-006, 


4/3 7/2 2 
é’ =log EZ ps E \ 


ZO. x peat Rn Ny Se 
=3 62-+4 og { We Bott: 
Here x’ is the probability of a particle, on traversal of the foil, obtaining a 
deflection through >|.x,| and can be adjusted to something like 0-05. 
In units of (k's)? =5€’1/2 (17) becomes 


1 G45 
fe) =a il+ Gt at. eee (18) 
‘The first term is the projected single scattering, and the other terms are corrections. 

In a recent treatment of multiple scattering, Snyder and Scott have given an 
asymptotic formula for the projected scattering on a plane, valid as long as the 
deviation from the projected single scattering is less than about 20°%. Using 
the ordinary unit of angle (i.e. the radian) their formula can be written 


Rey oe 66? heal 4554 | x| a 
hes ove a3 {1 Bien (10g 5 —0 8st) {eae (ioe A +0-616) 


; Lal 
x 0 


~ 3-85 log (7 +0616) +3:33 | = Pe iy nen (19) 
0 
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In Figure 2 the percentage deviations as calculated from (18) and (19) are compared 
for the case of scattering of 20 Mev. electrons by 0-005 cm. of iead. It is seen 
that the Snyder—Scott formula gives slightly larger deviations than does (18). 
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Deviation 


A 


0 4 i L teas un 
0-3 0-4 0-5 
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Figure 2. The percentage deviation of the projected scattering from single scattering for the 
case of 20 Mev. electrons scattered by 0-005 cm. of lead, as calculated (1) from equation (18) 
(2) from the Snyder—Scott formula (19). 


, 


It is of interest to note that we can force (18) into a form almost identical 
with (19) by choosing x, in fact equal to the actual angle x at which a measurement 
is made. In this case €’ =log! «| /2:000), and with the ordinary unit of angle, 
(18) becomes 


se: 65? eles, 4564 Fo en 2 
{(%;5)= za! a (toe —0 698 + SF (los —0 698 ‘i fee 


The similarity between (19) and (20) is evident, and on plotting percentage 
deviations they are found to be almost indistinguishable. However, when x, 
is adjusted in the manner discussed in $3, it is always much lower than angles 
for which the distribution approaches single scattering, and we would expect 
therefore that (18) be more accurate than (20). 
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Penetrating Showers at High Altitude 
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Communicated by P. M.S. Blackett; MS. received 1st September 1949 


ABSTRACT. A study has been made of cloud chamber photographs of penetrating 
cosmic-ray showers, originating in lead placed over the cloud chamber and in a lead plate 
within the cloud chamber, obtained at altitude 3,572 m. in a magnetic field of about 1,000 
gauss. 

The use of well-separated successive lead layers allowed the association of penetrating 
events to be demonstrated, and there is some evidence of a general reduction of energy 
involved in events in the lower lead layers as compared with those in the top layer. Mesons, 
protons and particles of higher charge were identified in the penetrating showers, and a 
mean range for non-coulomb scattering of shower particles is given. Evidence is also 
given of the apparent starting and stopping of charged particles in the lead plate at the centre 
of the cloud chamber. 


Si JONAM RODS CALLIN) 

HE existence of penetrating cosmic-ray showers was first demonstrated by 
Wataghin and his co-workers (1940) and by Jénossy and his co-workers 
(1940), since which time these showers have been the subject of many 
detailed researches. In particular Broadbent and Janossy (1947) distinguished 
two types of penetrating showers: the local penetrating showers which they 
considered to be produced by isolated nucleons entering the apparatus, and 
extensive penetrating showers in which the penetrating events are associated 
with extensive cascade showers in the atmosphere. In the early work, the 
multiplicity achieved in the local penetrating showers was mainly attributed to 
the creation of mesons, and the transition curves for these showers corresponded | 
to a primary range which was in reasonable agreement with the theoretical value 
given by Hamilton, Heitler and Peng (1943). Recently, the more detailed work 
of Rochester and Butler (1948 a,b) has thrown doubt on the simplicity of the 
picture and in particular has shown that many of the fast particles emitted in the 

local showers are not ordinary mesons and may well be protons. 

The present experiment was designed to take advantage of the rapid increase 
in the nucleon intensity with altitude in order to study large numbers of penetrating 
shower events. It was also hoped that if the events attributed by Rochester and 
Butler (1947) to the decay of heavy particles increased with altitude in the same 
ratio as the local penetrating showers, further examples of these processes would 
be observed. Work was accordingly carried out at the International High 
Altitude Research Station at Jungfraujoch (3,572 m. above sea level). 


§2. APPARATUS 

A counter-controlled cloud chamber, 27 cm. diameter and 7:5 cm. deep, and 
containing a lead plate 2:6 cm. thick, was operated in a magnetic field of 900 gauss. 
The chamber was contained in a thermally isolated brass box, the whole apparatus 
being housed in an aluminium cork-lined hut on the roof of the Sphinx meteoro- 
logical observatory at Jungfraujoch. Although the hut was thermostatically 
maintained at a uniform temperature, considerable difficulty was encountered in 
maintaining thermal stability in the cloud chamber. This was due to the large 
temperature fluctuations of the external air brought into the hut to cool the field 
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solenoids. After the first short period of operation, however, in which track 
distortion was severe, a condition was reached in which the oeinanin detectable 
momentum in either half of the chamber was about 2 x 108 ey C 

The counter arrangement used, which followed closely an arrangement used 
by Rochester and Butler (1948 a) is shown in Figure 1. A three-fold counter- 
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Figure 1. 


telescope, built of six 20cm. x 2cm. counters, was surrounded by lead (10cm. 
above, at least 5 cm. at each side) and selected particles in the central plane of the- 
chamber, which was placed below a further lead layer varying in thickness from 
lcm. to about 3cm. Below the chamber, and without further shielding, was 
placed a two-fold counter tray of four 30cm. x3-5cm. counters. The five-fold 
coincidences of the upper shielded telescope and the lower two-fold tray were 
used to operate the chamber. This system of selection was, of course, sensitive 
to ‘knock-on’ electrons produced by fast incident particles in the lead plate inside 
the chamber and, to a certain extent, in the lead immediately above the cloud 
chamber. It is also sensitive to multiple showers of penetrating particles in any 
of the lead layers: in the lower layers an ionizing primary is required, while in 
the top lead layer multiple showers arising from a non-ionizing primary would also. 
be recorded. 
§3. EXPERIMENTAL DATA 


3.1. The Counting Rate 
At sea level (Manchester), the counting rate of the five-fold counter system 


(Figure 1) was 6:3 counts per hour. At Jungfraujoch 4,700 photographs were 
taken in 642 hours, i.e. at a rate of 7-3 photographs per hour. In order to compare: 
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“these rates it is necessary to know the inert time of the cloud chamber, during 
which cleaning expansions were made and thermal equilibrium restored after 
each operation. In our apparatus this inert time was 4 minutes. In Table 1 
the rates quoted above and the simple counting rate recorded at Jungfraujoch 
-are shown, together with the rate of photographing at sea level deduced from the 
value of the inert time. 


Table 1 
Counting rate Photographing rate 
Sea level 6:3 per hour (4-4 per hour) 
Jungfraujoch 14-3 per hour 7:3 per hour 


It will be seen that, although the intensity of penetrating showers at Jungfrau- 
_joch is expected to be many times that at sea level, the actual rate of photographing 
has only increased 1-7 times. 

Suppose the counts to arise only from knock-on electrons, which will be 
produced almost entirely by the ordinary meson component and from penetrating 
‘showers. While the latter increase in intensity about 28 times between sea level 
and Jungfraujoch (George and Jason 1948), the former will scarcely increase at 
all. This small increase of knock-on counts is due to the fact that only those 
knock-on events have a large probability of operating the counter system in which 
the electron is projected at a relatively small angle to its primary, and with radius 
-of curvature at least greater than that of the cloud chamber. Such knock-on 
-events have been shown (Wilson 1938) to arise almost exclusively from mesons 
of large momentum (>3 x 10% ev/c.), and we estimate their increase in number 
to be not greater than 1-4 times. 

If, then, Cy is the counting rate at sea level due to knock-on events, and Cp 
‘that due to penetrating showers, we have 


Cyt+Cp=6:3, (Sea level) 
1-4C, +28Cp=14-3. (Jungfraujoch) 
‘Whence 
Cp/(Cy +Cp) =0-03, 
28Cp/(1-4Cy +28C>p) = 0-40. 


‘Thus, while the rate of photographing increased 1:7 times, it appears likely that 
‘the probability that any operation has been initiated by a penetrating shower 
has increased from 3% to 40°, and the resultant gain is about 13 times. 


3.2. The Experimental Results 


Of 4,700 photographs taken at Jungfraujoch 92 showed ‘ penetrating showers’ 
-according to the definition of Rochester and Butler (two or more particles pene- 
trating the 2-6 cm. lead plate without measurable loss of momentum and leading 
‘to no secondaries), and the main result of the present work is concerned with these 
“events. We shall discuss the angular distribution of the penetrating showers, 
the association of penetrating showers and of comparable nuclear explosions, 


and the evidence from our photographs of the nature of the penetrating particles 
in the showers. 
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3.3. Angular Distribution of Particles in Penetrating Showers 

‘For all penetrating showers, the tracks observed were drawn on a scale 
diagram of the cloud chamber and the various lead absorbers; it was then found 
possible, by extrapolating the tracks backwards, to determine the layer of lead 
in which the shower appeared to have originated. In this way 58 showers were 
referred to the upper lead layer T and 34 to the layer M immediately above the 
cloud chamber. ‘The second group formed in M were particularly well defined,. 
and in them the tracks projected back to a point of formation; in the first group: 
the projection was not so successful, probably on account of scattering of the 
particles in the lower layer M, but this uncertainty was not of a kind which could 
cause these showers to be attributed, erroneously, to a point of origin in the- 
lower plate. 

On inspection, it appeared that the showers from the upper plate T are, on 
average, more collimated than those from the lower plate M, but it is extremely 
difficult to establish this fact unambiguously on account of the great diversity of 
type among the showers. However, several examples existed, for each layer of 
formation, of showers in which two penetrating particles reach the cloud chamber; 
the projected angular separation of these penetrating pairs of particles was. 
compared with the distribution of projected angles which could be imposed by 
the geometry of the cloud chamber, the absorbing plate and the counter array for 
a distribution of formation which is isotropic at least for the limits of angular 
spread set by geometry. The results of the comparison are given in Table 2. 


Table 2. Angular Distribution of Penetrating Showers of 
Two Particles in the Cloud Chamber 


Layer of Mean angular separation (deg.) 
Origin Neos Observed Calculated (isotropic) 
ip it 4°5 10 
M 7 18-4 20 


We thus have some evidence that the showers from the upper layer 'T but not 
those from the layer M are appreciably collimated over angles of the order 10°. 
It is therefore probable that the energy involved in nuclear collisions in the lower 
lead layer is significantly less than that in the top block of lead, and that the nucleon 
spectrum undergoes a transition effect entering lead from the atmosphere. This: 
effect might be expected to arise if the multiplicity in nucleon showers is a function 
of the number of nucleons in the nuclei in which collision takes place, provided 
that some at least of the secondary particles in showers are themselves capable of 
shower production. 

The angular distribution of penetrating particles in some events is clearly 
most complex. In one example of an explosion in the central lead plate (Figure 2),. 
which is probably initiated by an uncharged primary, the penetrating particles. 
emerging fall into three groups: a very narrow group, a, of at least four particles. 
within a total angle of less than 10°, a second group, b, of at least five particles also: 
in an angle of about 10°, the axes of the groups being separated by about 40°, 
and a third group, c, also narrow and containing about four particles with axis. 
about 20° to that of a. In addition one or two other particles make considerable: 
angles with the main groups. In such an event the final distribution of particles. 
can hardly be regarded as the result of a single stage of interaction, and any” 
argument based on one section only of the event is liable to be misleading. 
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3.4. The Association of Penetrating Showers 


Our photographs give evidence that particles in a penetrating shower can 
themselves lead to further penetrating showers, and are thus in agreement with 
recent observations (Tinlot and Gregory 1949, and Chao 1949) which show 
considerable association of showers of the nuclear collision type. 

Penetrating showers originating in layer M were observed to be associated 
either with showers originating in the upper layer 'T (Figure 3), or else with single 
penetrating particles passing through M at a point other than the point of origin 


a 


Figure 2. 


Figure 3. 


of the shower but in the same general direction. The occurrence of such a 

combination as a random association of the single particles with the penetrating 
--event in M is extremely rare (about one, only, in all photographs of penetrating 
events starting in M); hence we regard these also as indicating the production of 
a secondary nucleon shower in M. 

Nuclear explosions in the lead plate C inside the cloud chamber, whether 
produced by ionizing or non-ionizing particles were similarly observed to be 
.associated with single penetrating particles traversing some other part of the plate 
C or with penetrating showers originating in the lead layer Tor M. The occur- 
tence of events in these categories is shown in Table 3. 


Table 3. Number of Showers associated with pre-existing 
Penetrating Particles (N,) 


Summary Np 

Showers in plate M (Total no. 34) 9 

Explosions in plate C (Total no. 32) ; 12 
Showers in M Detail 

Possibly large scattering of particles from T 3 

Shower +1 other penetrating particle 3 

Shower-+shower from T 3 


Explosions in C 


Explosion+1 other penetrating particle 

Explosion+1 other penetrating particle+1 particle stopped 
Explosion-+shower from T 

Explosion -+shower from M 


wud bd 
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3.5. The Nature of Penetrating Particles 


(a) Particles showing increased ionization in traversing plate C. 
Five penetrating shower particles (Table 4) showed increased ionization after 
traversing the lead plate C in the centre of the cloud chamber. 


Table 4. Estimated Ionization of Particles traversing Chamber 
and appearing to show Increased Ionization 


"Top half of chamber L/Imin 2 D 2 Sait 1% 
Lower half of chamber Vi Gees 6 6 6 3 3 


* Although for the last two particles the ionization in the upper half was judged to be 
greater than minimum, little weight can be attached to these estimates. 


The expected behaviour of particles of mass 200m,, 300m,, 1,000m, and 
1,840 m, over a range relevant to the first three particles of Table 4 are shown in 
‘Table 5, from which it can be concluded that the masses of these particles are at 
least of the order 1,000 m,. 


Table 5. Ionization of Particles traversing 2-6cm. Lead 


Mass Ionization as multiple of minimum 
200776 ea oe ms Rees 
300m - oe a ne eee 
1000me se io o ne 
1840 se i a oe 


(b) Particles which are heavily scattered. 

The projection of the angles of scattering in a plane perpendicular to the axis 
of photography of 230 shower particles traversing the lead plate in the cloud 
chamber is summarized in Table 6. 


Table 6. Scattering of Particles in Penetrating Showers of Momentum 
>2 x 10% ev/c. traversing 26cm. Lead in Cloud Chamber 


Angle (degrees) 0-5 5-10 10-20 20-30 >30 
Observed no. of particles 183 18 Beis 6 5 
‘Calculated number 183 JA8s ©. 13 ep 0-04 


The row ‘ Calculated number’ refers to the assumption of coulomb scattering alone 
when all particles scattered more than 5° are of the lowest possible momentum. 


Cc 


The momentum of the particles concerned is in no case less than about 
2x108ev/c. at which momentum the probable angle of multiple coulomb 
scattering in 2-6 cm. lead would be 7°. ‘Thus with the assumption most favourable 
to the presence of coulomb scattering only (that is, that the whole of the group of 
particles scattered in the range of angles 5—10°. are of momentum 2 x 108 ev/c. 
and suffer only coulomb scattering, while the majority of particles scattered 
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0-5° are of much larger momentum), between 5 and 15 particles of those scattered! 
more than 10° must be considered to show additional scattering of non-coulomb: 
type. We thus find that the particles in the penetrating showers are, on average, 
heavily scattered once in traversing 40-100 cm. of lead, a scattering range which may 
be compared with that of 25cm. lead reported by Rochester and Butler (1948 b). 
The longer range in the present work is essentially an upper limit, and has almost 
certainly been increased, as compared with the figure given by Rochester and 
Butler, on account of the much lower maximum detectable momentum obtained 
in our apparatus. It is clear, however, that many of the particles near minimum 
ionization in the penetrating showers are not j-mesons. 


(c) Particles stopped and created in the lead plate C. 

Out of 247 particles with momentum greater than 2 x 10% ev/c. in penetrating 
showers, 17 showing minimum ionization were stopped in the lead plate C. 
Twenty other particles appeared under the lead plate without a visible origin 
above. Reprojecting these photographs through the same stereoscopic lenses used 
in photographing, the position of these tracks in the cloud chamber was located 
in order to establish that with normal scattering they would not pass out of the 
illuminated region. These particles must either be scattered through very 
large angles or be genuinely stopped or created in the lead plate. Projections in 
two planes of photographs showing tracks stopping or starting in the lead plate 
are given in Figures 4 and 5, and show particles stopped or starting which are well 
inside the illuminated part of the cloud chamber. 


B b 


Figure 4. Figure 5. 


We can exclude the possibility that the particles which stop are either /-mesons 
or electrons, the former because an initial momentum of only about 108 6V/Cuis 
required to traverse the plate, and the latter because from an entrant electron 
p>2 x 10% ev/c. several shower particles are to be expected below 2:6cm. lead. 

It is therefore probable that the particles are of a kind interacting strongly with 
nuclei, and that they either undergo an exchange transformation of the type 
proton-neutron in a collision in the plate, or else dissipate the incident energy in 
such a collision. Similarly, the particles created in the plate are likely to arise 
from nuclear explosions initiated by a non-ionizing primary or else be derived from 
a non-lonizing primary by an exchange process. 
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(d) Heavily ionizing particles. 


A total of 94 heavily ionizing particles were observed, 49 in the upper half 
of the cloud chamber and 45 in the lower half, under the plate C. The 
circumstances of occurrence of these particles are summarized in Table 7. 


Table 7. Summary of Occurrence of Heavily Ionizing Particles 
1. Particles in Upper Half of Chamber 


No. associated with penetrating particles 21 
No. associated with soft component 14 
No. identified as multi-charged particles 4 
No. identified as jz-mesons 1 
No. not associated with other events in the chamber 9 

Total 49 


2. Particles in Lower Half of the Chamber 


Single particles associated with explosions 10 
Single particles associated with single penetrating particles 4 
No. of heavy pairs created in the lead plate (9) 18 
No. of particles becoming heavy in traversing the plate 5 
No. of particles occurring in showers of more than two particles 6 
Single particles not associated with other events 2 

Total 45 


The great majority of these particles show ionization estimated to be greater 
than five times minimum and have no measurable curvature: all of these may well 
be protons. However, two particles occur which are certainly light (z- or p-) 
mesons, and five particles which lead to 6-electrons of measurable range will be 
shown to be multiply charged. 

One of the tracks identified as a meson decays to a lightly ionizing particle 
in the upper half of the cloud chamber, and must be identified as a jz-meson. 

The other meson track, which occurs in a explosion from 1cm. depth in the 
plate C initiated by an ionizing primary, is identified by ionization and momentum. 
‘The measured momentum is (2:2 + 0-7) x 10’ ev/c. and the ionization is estimated 
as 5 times minimum. ‘This particle is therefore not identified with certainty, 
but may be either a ~p-meson or 7-meson. 

Several heavily ionizing tracks are of particular interest because the formation 
of $-electrons along the trajectory in the cloud chamber allows an estimate of 
particle charge to be made. 

The probability of production of secondary electrons in elastic collision is 
given (Janossy 1948, p. 94) by 

WE) dE =20rZ2Nme2(Z/pyPdk/E*, . —....... (1) 
where (E) dE is the probability that a particle of charge Ze traversing material 
of atomic number z and N atoms/cm? will give rise to a secondary electron in the 
energy range E to E+dE in one centimetre of path length, and where %,m, c, B 
have the usual significance. ‘The number of 6-electrons per centimetre of path, 


say n, is then given by 
n =2nr,22N me(Z/py | 


-Emax 
dE/E?, pee (2) 


Emin 
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where E,,,,, is the minimum energy of 5-electrons which are selected, and Ejnax, 
the maximum transferable energy, is 2me?6?/(1—?). Equation (2) allows us 
to write Z2/n as a function of 6? for the particular medium traversed (in our chamber, 
oxygen at about 1 atmospheric pressure), and for any chosen value On Lies 
typical curve of Z2/n plotted against 2 is given in Figure 6 which refers to 6 rays 


160 
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Figure 6. 


of minimum energy 15kev. ‘The charge of the fast particle then lies between 
limits at the minimum of this curve and the value of Z corresponding to B=1. 
For all of our examples the low velocity limb of the curve is excluded by the range 
over which the particle has been observed. Results for five particles are given 
in Table 8. 


Table 8. Charge Determinations for Heavily Ionizing Particles 
producing 5-Electrons 


(1) (2) (3) Zmin Zmax 
5 a NE 3 12 

6 6 2/6 2 8 
6 6 2/6 De 8 
4 3 4/4 3 11 
855 3 5/8°5 2 9 


(1) Length of primary track available (cm.); (2) minimum range of 5-electron selected 
({mm.); (3) n=(no. of d-electrons)/(cm. track avalable). 


No account has been taken in this discussion of the statistical weight to be 
given to the very small number of secondary electrons observed for each heavy 
track. However it must be noted that the quantity deduced is actually Z2/n, and 
so for the first, fourth and fifth tracks, at least, there is an extremely low possibility 
that a singly charged particle will give rise to the observed secondaries. 
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§4. DISCUSSION 

The most striking feature of our photographs is the extent to which the results. 
of nuclear interaction, whether they are classified as penetrating showers or as 
explosions in the lead plate inside the cloud chamber, are associated. The 
distinction of successive events is assisted by our experimental arrangement in 
which the lead over the cloud chamber is in two well separated layers: with this. 
arrangement 21 out of 66 events taking place in the lead M immediately over the 
cloud chamber or in the plate C inside the chamber are associated with penetrating 
particles in addition to their own primaries. Similar conclusions have been 
reported by Tinlot and Gregory (1949) and Chao (1949). 

The result of the entry of a high energy nucleon-like particle into a layer of 
matter is thus in some ways similar to an electron—photon cascade, although the 
process is simpler in so far as multiplication takes place only at a very limited 
number of collisions. However, the degree of multiplication at a collision as a 
function of the energy of the incident particles and of the nature of the nucleus 
which is disrupted is still unknown, and observations have not yet allowed us to: 
estimate how many of the particles produced in a penetrating shower are of a kind 
with strong nuclear interaction which may be expected to lead to further stages of 
multiplication. Attention has been drawn to the multiplication process by 
Rossi (1948). . 

When the detailed behaviour of collisions leading to penetrating showers is. 
understood, it will no doubt be possible to formulate transition relations for a 
beam of nucleons (perhaps together with other strongly interacting particles) in 
equilibrium, analogous to those familiar in cascade theory. In this connection 
Rossi shows that the energy spectrum of protons (kinetic energy greater than 
about 100 Mev.) changes only slightly in the greater part of the atmosphere, and 
this result must certainly be regarded rather as an equilibrium condition than as a 
pure absorption of the primary beam. 

We have given some evidence of a change of energy distribution in the nuclear 
beam on traversing successive layers of lead, in so far as the mean angle between 
pairs of penetrating particles from the upper lead plate appears to be smaller 
than the corresponding angle for pairs from the middle plate. ‘This result would 
indicate a transition effect towards a different equilibrium particle spectrum in 
lead. Such a transition might arise either as a result of the change in multi- 
plicity of projected nucleons between air and lead or perhaps as a result of the 
presence, in the condensed material, of unstable particles capable of nuclear 
collisions which in air would decay before making further collisions. 
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ABSTRACT. Models, largely based on the assumptions introduced by Peierls and 
Nabarro in dealing with a single dislocation, are used in calculations on three types of 
intercrystalline boundaries, namely, (I) a boundary due to a difference of atomic spacing, 
(II) a twist boundary, and (III) a symmetrical tilt boundary. With these models the 
resolution of the boundary into a sequence of dislocations is a natural consequence of the 
analysis, which also yields the expressions for the stresses, atomic displacements and energies 
associated with the boundary, as functions of the angle of tilt, etc. By allowing the distance 
between dislocations to tend to infinity, these expressions reduce to the corresponding 
ones for single dislocations. ‘The outstanding feature of the interfacial energy is that it 
increascs initially very rapidly with the angle of tilt, etc. An application of the results to 
the theory of oriented overgrowths, developed by Frank and the present author, is described. 
"The validity of the assumptions and approximations involved and the advantages of the 
treatment as compared with those of other writers are discussed. 


§1. INTRODUCTION 
OR the solution of a number of problems in the physics of solids a knowledge 
of the energy associated with an interface of misfit between two similar 
or dissimilar crystals would be of considerable assistance. The purpose 
of the present paper is to calculate the stresses and energy associated with such 
aboundary. Three simple types of boundary are treated. These are: 

I. A boundary at which the lattice spacing changes: e.g., that formed by 
two crystals of different lattice spacing in contact. 

II. A twist boundary, i.e. the boundary between two halves of a crystal 
rotated with respect to each other about a common axis perpendicular to the 
boundary. 

III. A tilt boundary, which occurs when two portions of the same crystal 


are symmetrically rotated with respect to each other about a common axis lying 
in the plane of the boundary. 


* Now at National Physical Laboratory, Pretoria, South Africa. 


Inter-Crystalline Boundaries On 


Case III has been treated by Burgers (1939, 1940), who points out that the 
surface between two crystals can be treated as, and is geometrically equivalent to,. 
a row of edge dislocations. In cases I and II the surface can similarly be treated 
by regarding it as a system of dislocations. This geometrical equivalence of a 
crystalline boundary and a system of dislocations has been discussed by various 
authors, namely, ‘Taylor (1934), Bragg (1940), W. G. Burgers (1947), Shockley 
and Read (1949) and Frank (1950). 

Burgers gives no detailed discussion of the singularity, i.e. the infinity, which 
occurs in the expression for the stress when the surrounding material is treated 
as a continuum, and he is thus prevented from making an estimate of the energy. 
Shockley and Read obtain an estimate of the energy by assuming that the energy 
associated with an element of volume at any point which is nearer to a dislocation 
than half the distance between dislocations does not vary with this distance, and 
that the strain energy outside these regions is negligible. Further, they do not 
calculate the energy associated with the region around the singularity. 

We avoid these difficulties by an appropriate choice of the model, depending 
on two basic assumptions: (a) that, at the interface, the potential energy of an 
element of area of the one may be taken to be a sinusoidal function of its position 
relative to an element of area of the other; (5) that the material elsewhere may 
be treated as an elastic continuum which obeys Hooke’s law. 

The assumptions are essentially those used by Peierls (1940) and Nabarro 

1947) in dealing with a single dislocation. By using them we introduce certain 
limitations and approximations. We shall represent the surface interaction by 
a single sinusoidal function for each of two rectangular position coordinates in 
the plane and thus restrict ourselves to cases in which the crystal surfaces in 
contact have square or rectangular symmetry. The magnitude of the error 
involved in adopting a single sinusoidal term instead of a Fourier series to 
represent the interaction energy in I and II has been discussed by Mackenzie 
(1950). 

With this model the resolution of the boundary into a sequence of dislocations. 
is a natural consequence of the analysis, and no singularities occur in the expressions. 
for the stresses. The energy density at any point is found to vary with the angle 
of tilt or its equivalent in the other two cases. 


I. Bounpary DUE TO DIFFERENCE OF ATOMIC SPACING 
$2 DEES MODEM 

The first problem is to extend the Peierls~Nabarro model of a single dislocation 
in a perfect crystal to the case where the atomic spacing on one side of the slip 
plane is different from that on the other and a parallel sequence of dislocations. 
exists in the slip plane. 

Such a case is that which arises when two different crystal lattices meet at an 
interface, e.g. a metal and its adherent oxide film or in general a crystal overgrowth 
on a different substance. 

The model used in the present analysis may be described as follows. Let 
A and B (Figure 1) be two semi-infinite crystals meeting in the surface S, which 
is perpendicular to the plane of the paper. The meshes in the atomic planes Sy 
and S,, adjacent to S are rectangular. Let d be the equilibrium distance between 
S, and S,. Cartesian coordinate axes are taken as shown in the figure, the 
y axis being perpendicular to the plane of the paper. At S the x spacings of 
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A and B atoms are different. Let these spacings be a and 6 respectively. Wet 
the lattice be such that Ph=(P+1)a=p, where P is a positive integer. If the 
atoms C and D are in alignment, then E and F, at a distance p, will also be aligned, 
and so on. For simplicity, misfit in only one direction has been introduced. 
However, because the equations are linear to the degree of approximation used, 
solutions may be superposed and the following treatment easily extended to cover 
the case of misfit in both the x and y directions. This will be discussed in greater 
detail in § 6. as 

In order to express the boundary conditions in a simple manner, it 1s convenient 
to introduce a reference lattice L having an x-spacing c. ‘The lattices A and B 
may be regarded as having been generated from L by a homogeneous compression 
of one half, in which the distance (P+1)c becomes (P+1)a=(P+3)c and a 
homogeneous extension of the other half in which the distance Pe becomes 
Pb=(P+$)ec. The molecular forces are assumed to be zero during this 


0 . 
i 
1 


jo are eS 


Figure 1. Jnterface between crystals A and B, having different atomic spacings. 
Dislocations are shown at Q, O, R 


hypothetical deformation and to act again on its completion. The algebraic 
difference of the final displacements with respect to the reference lattice of an 
atomin Spandthenearestatomin S, willberequired. Thisrelative displacement, 
which we denote by U, is seen to increase in every interval p by c. For example, 
if C and D remain in alignment throughout the deformation, the relative 
displacement between E and F is c. 

The atoms in Sp, are subjected to two forces: (i) the interaction with atoms 
in A, particularly those in S,, which tend to bring the atoms in S, into correct 
alignment with the nearest atoms in S,, and (ii) the interaction with other atoms 
in B, which tend to keep them in their natural positions. In equilibrium the 
resultant of these two forces vanishes and the regions of misfit, for example the 
one at the y axis, are localized. The y axis is thus the axis of a dislocation. 
Since the dislocation lines extend throughout the crystal parallel to the y axis, 
the problem will be one of plane strain. The strains are furthermore periodic 
in x with a period equal to the distance p between dislocations. The lattices 
of A and B will be slightly distorted under the influence of the force (i), particularly 
in the neighbourhood of the dislocations. We denote the x and z components 
of the displacement of a point in the distorted lattice with respect to its position 
in the undistorted lattice by u and w respectively. If the radial extent of the 
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strains associated with the dislocations is large compared with the interatomic 
distance, then both w and w will vary slowly from atom to atom. The relative 
displacement of neighbouring atoms in each half-crystal will then be much 
smaller than a (or 6). (This assumption will be shown later to be reasonable.) 
In these circumstances each half-crystal may be treated as an elastic continuum. 
Moreover, for the sake of simplicity the continua will be assumed to be elastically 
isotropic and to have identical elastic moduli. 

We assume that the « component of the force (i) on any particular atom 
depends only on its tangential separation from the nearest atom in the other slab, 
and that the z component is similarly a function of the x component of relative 
displacement only. The tangential forces acting on corresponding points of 
the opposed surfaces are therefore equal and opposite, in which case their 
2 component of relative displacement is zero and the z component of (1) 
vanishes. ‘This assumption neglects the change in the equilibrium value of their 
2 component of separation where the atoms are out of alignment in the two planes. 

It is convenient to express the « component of the force (i) as a function of U, 
and we note that it has the period c. We assume the function to be sinusoidal. 
This is equivalent to neglecting all but the first term of a Fourier series. Some 
aspects of this approximation have been discussed by Mackenzie (1950). If U 
varies slowly from atom to atom we may regard S, and S, as continuous and express 
the force in terms of a shear stress acting on Sp, namely, 


PEN = (ee cae yISit (ZT) 66 2 ace ee (1) 
where the constant jy is of the order of magnitude of the shear modulus p. If 


d=c and the sine law is a sufficiently good approximation, ip is equal to yw, in 
accordance with the treatment of Peierls and Nabarro. 


§3. THE GOVERNING EQUATION AND ITS SOLUTION 

The problem is now that of setting up the equation governing the equilibrium 
of the forces (i) and (ii). This problem, which according to $2 involves plane 
strain, may be solved by expressing the stresses in terms of a stress function x 
in the usual form (Timoshenko 1934) : 

aX ae eae a) 
Daa yp? Pos 2? Pea = axoz” Nae 

Furthermore, in the case of plane strain, the strain components can be expressed 
in terms of the stress components by an equation 


2Hlxe = (1 re a) Rae ber eee ee (3) 
where o is Poisson’s ratio and e,,,, is the x component of strain. 
x must further satisfy the equation V'y=0. The general expression for a 
particular solution of this equation may be written in the form (Timoshenko 1934) 


NetB paCjer™ (ine), oe (4) 
where m, B and C are constants. If we take into account that the stresses must 
vanish at x = = © and that p,,, vanishes at x = +$p on S, and Sp while p,, vanishes 
everywhere on S, and S, in accordance with the assumptions in § 2, we may choose 
the following proper form of x: 


y=6Cexp(F2anl/p)cos(2amx/p), wv eee (5) 
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where n is a positive integer and where ¢ is measured from S, for A and Sz for B. 
The upper sign refers to B and the lower one to A. This sign convention will 
be used throughout. (Where the signs are equal one will be omitted.) 

From (2) and (5) we now find 


pee = (- a a 2) exp ( + ) C cos oe ) 


as Qa 
Pes Th = exp ( z a) Ccos NX I a a Aa (6) 
Ci 7 7 alah 
ot oe a ease 


On the planes S, and S, equations (6) reduce to 


Pox = + (47n/p)C cos (27nx/p), 
P| oe cede a re ee 8 eg ct (7) 
and Pex =(270/p)C sin (27nx/p). 


If we substitute (7) into (3) and carry out the integration, choosing the integration 
constant such that wu is zero at x= +4), we obtain 


u=F{(1—o)[p}Csin(2anx/p), ne we ae (8) 


This is seen to be in accordance with the fact that the displacements must be 
equal and opposite, as follows from the assumptions of equal and opposite forces 
and identical elasticity. It is convenient to put 


Ca —{uclrl —«)4 i a eee (9) 


where the constant A, depends on. Finally we may express the most general 
form of (8) in terms of the Fourier series 


Cue eas 27x 
u= + Teer > ean eteneus (10) 
On using this we find 
Epes iets Scalia (11) 
P Ty=1 0 Pp 
where U is taken to be zero at x= —}p. The first two terms in (10) represent 


the linear variation which exists when the lattices A and B are rectangular, 
i.e. before the force (i) is allowed to act. If we substitute (9) into (7) the general 
expression for the shear stress is seen to be 


< 2nnx 
me (Une a EE SA si ee ae 12) 
Pp r( ) PU 9) na a p ( ) 


The problem is now that of finding U from (1), (11) and (12). We put 
B =ape/pu(l1—c), X =2ax/p, b=Ujle ~~ sae (13) 


Pre. tee” 


ee 


-!- 


_- 
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in which case the differentiation of (11) and the combination of (1) and (12) 
respectively yield 


2 ) 
Pa be cosmk =n (FE ) 


n=1 aX W Qn 
‘ PaO sat ces ns: (14) 
and eae A, sinnX = — ag sin Qin, | 


If we represent the function conjugate to f(X) by Conjugate {f(X)} and make use 
of the fact that a function and its conjugate are respectively the real and imaginary 
parts of an analytic function of a complex variable, then it is seen that 


n 
-10 eA Oe 10 
w/a 
Figure 2. 


Curve A. Relative displacement of an atom in Sp with respect to the nearest atom in Sa (U/c=y).. 

Curve B. Displacement of boundary perpendicular to itself (w/c). The curves are calculated for 
the case in which each dislocation is 20 atoms from the next (p=20a). co is taken to be 0:3° 
and py to be equal to up. 


Conjugate {cos mX}= sinmX and Conjugate {sinmX}=—cosmX. Byusing this. 
we may finally write the governing equations (14) in the form 


Lar (4 — x) =Conjugate {B-'sin2mp}. «2 (15) 


This can be expressed as an integral equation (Hardy and Rogosinski 1944) : 


CS) Maa ey ein Qrl(X = 1) cot itd 
dn (oe — ze) ze lk {sin 2my(X +#)— sin 2nyi(X — 1) cot H dt. 
pieced (16) 


A method of iteration yielded the solution of this equation as a power series in 1/.. 
The sum-function 


%=4+4(1/m)arctan {[((1+B +f ]tangX}  ...... (17) 
of this series is defined for all values of 8, whereas the series converges only for’ 


a limited range of the variable. A graph of (x) is shown in Figure 2.. 
Differentiation of (17) gives 


di jdX =(1j2n){((1 +B) BteosX}2. ws. (18) 
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This solution can be verified by substituting (17) and (18) into (15). We find, si 
as we should, that B-!sin2m and 2 (dp/dX —1/27) are respectively the real } 
and imaginary parts of an analytic function of a complex variable. This function : 
is 


=271 =[(T ep yh esl a eee (19) 


The Fourier coefficients A,, can now be obtained by expanding (19) into a } 
Fourier series and comparing the coefficients of its imaginary part with those 


of 2x (db/dX —1/2m) in (14). We find 
A, ={(1 4B) BSA Gay). ee (20) ( 


§4. STRESSES AND DISPLACEMENTS 


If we use the notation 


X =270)| Pow ands Zi 276 ge eee eee (21) 


and substitute (9) and (20) into (5) and (6), we may deduce 


__ pe Ss A” ene X 

XS pearpeeeye a er” Cosi xX. 
ead = aa FnZ M 
Per = — SO CU ANZ F 2Z)e*”” cos nX, . 
pl =6) n=1 7 
bite aec anata (22) i 
Pz= eee. SZ nAre®™ cos nX | 
e p(l—c) n=1 t 

and Pex = — Aa An(1 = nZ)et™ sin nX. 


PU ein j 


These Fourier forms of the stresses are useful in calculating the displacements 
and the strain energies in each half-crystal. 

The expressions for the sums of these series are also useful. The sum for x 
can be written down from the sum of the complex series of which y is the real 
part. ‘The other sums concerned are easily treated by differentiation of this 
series. ‘The expression for the stresses may, of course, also be determined from x 
by differentiation according to (2). We find 


x = {ucp/472(1 —c)}Z In R, 
Pow =— Q{Z[(1 + A*eF”) cos X —2Ae™“] F 2R® (cos X —AeF?)}, 
Pa=QZ{(1+ A?et”) cos X—2Ae*}, 
and Pex = —Qsin X{R?F Z(1 — A®eF*4)}, 


where QO =ycAe*7/(1—c)pR4 and R? =1 + A2et?4 — 2 Aer4 cos X. 
The displacements u and w can be determined by using (3) and the analogous 
expression for e.,. The integration of (3), for example, gives 


2p =(1—0) | Pye dx —o [p..dx +const. 
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The integrals in this and the analogous equation for w are most easily obtained 
by using (22). If we sum the integrated series we find 


Rate Aet4sin X cZe*4sin X 
Se) ae l—Ae**cosX) 4n(1—c)R? 
(1 —2c)c Rag\eecZ ae? (cosX Aer) 
d Se at Le Lk 
ee yen ag ae 


i a(1 = a) 
where the integration constants are chosen such that both u and w vanish at the 
dislocation line. 


The graph of 


w(x)/¢ = 


(1—2c) , [1+A?—2A cos (2nx/p) 
Sina) (— 4? 


shown in Figure 2 is of particular interest as it describes the shape of the boundary. 

An interesting feature of the present treatment is that it provides another 
method of deriving the stresses and displacements associated with a single 
Peierls-Nabarro dislocation. This is done by expanding the radicals, and the 
exponential and sinusoidal functions contained in (17), (23) and (24) as Taylor 


series and taking the limits of the resulting expressions as p tends to infinity. 
We find 


=} +(1/n) arctan (x/£5), 
w= oa In (F) = pee) igen ieee te (25) 


Pow =H2I(E + bo) — E1(E + Lo)? — x], 
Por =O5X(6 + Lo)” — x7} 
and Pew = Geir” —20(E + Lo}; J 
where g=pa/2r(1—o)r*, 7? =x? +(€ + Co)?, Sp =ap/492(1 —o) and a=b =c. 


Leibfried and Liicke (1949), who have also treated the Peierls—Nabarro case 
in which y=, discuss in detail the relation between its solution and the pure 
elastic solution. 


§5. THE INTERFACIAL ENERGY 


The energy associated with the boundary between the two half-crystals 
consists of two parts, namely, the strain energy in each half-crystal and the 
potential energy of attraction due to the sinusoidal force at the interface. The 
mean total energy per unit area of the interface is called the interfacial energy 
and denoted by r. 

The total strain energy residing in a slab extending a distance p in the 
x direction, unit distance in the y direction and to infinity in the positive and negative 
2 directions is given by 

7 00 
0 


Be in aie ep?) — 200 ap + ope dx dl 
E, =4| i ra o)(Prx Pz) OP yxPzz Da 
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where, as we are dealing with plane strain, the integrand is the strain energy 
per unit volume in each half-crystal (Timoshenko 1934). The integral is easily 


evaluated by using the series (24). We find 


ope AM Ce ie open on ee (26) 
PS Gelman 


The energy of attraction per unit area of the interface follows from (1), 


U Siac (are =e) 
| Pia(0, x) dU = 7 5(1— cos —— ). 


“0 
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Figure 3. Interfacial energy as a function of misfit (c/p) between the two surfaces in contact (misfit 


in one direction only). Energy in units of pc/4m? per unit area, where c is the lattice 
parameter, p the rigidity modulus of the bulk material and p the distance between dislocations. 


Curve A. The elastic energy of the two half-crystals (E¢/p). 
Curve B. Energy resulting from sinusoidal force at interface (E,/p). 
Curve C. Total energy of interface {r=(E,+£,)/p}. 


The total energy of attraction connected with the interface of the slab considered 
above is therefore 


0 da 


3p 
(ee es (1 — cos mp) de. 
— ip c 


If we substitute (17) in this integral and carry out the integration, we find 


E,=(Hocp/4n2)(1+8-(1 489). eee (27) 
The interfacial energy may now be written as 
T =(poc/477){1 + 8 —(1 +8?)? — Bln [28(1 +62)? -267]}, ..... (28) 


the logarithmic term being £,/p and the remaining part F,/p. The graph of 
7(c/p) in Figure 3 shows that 7 increases initially very rapidly with the density 
of dislocations c/p. 
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§6. A CROSS GRID OF DISLOCATIONS 


We still have to show that the preceding treatment also covers the case in 
which there is in addition a difference of spacing in the y direction, as will usually 
be the case. There is then an additional sequence of dislocations parallel to the 
* axis and associated with it a force analogous to (i) in §2. We take a term of 
the form (1) to represent the force acting parallel to each axis. This is equivalent 
to neglecting all but the first two terms of a surface Fourier series, each being a 
function of one coordinate only. Even when the sinusoidal approximation is 
not a good one, the resulting error in the interfacial energy will be small as long 
as the distance between dislocations is large compared with the interatomic 
distance, for then only a small proportion of the atoms at the interface are in bad 
alignment, i.e. outside the regions where the sinusoidal law approximates to 
Hooke’s law. 

To the degree of approximation of the sinusoidal law we may take U and V 
to be of the forms U(x) and V(y), in which case u and v must also be of the forms 
u(x, x) and u(y, 2), and hence w of the form w =w,(x, 2) + w,(y, 2), where V and v 
are analogous to U and u but refer to the y direction. That w has this form can 
easily be shown by substitution in the differential equations for u, v, and w. To 
the degree of approximation used, the equations are therefore linear, which allows 
the superposition of solutions and the addition of the energies associated with 
the two sequences of dislocations. 


II. Bounpary DUE TO ‘TWIST 
§7. THE MODEL 


The second problem is that of finding the stresses and energies associated 
with a twist boundary, i.e. the plane boundary between the halves of a crystal 
which are rigidly rotated with respect to each other about a common axis 
perpendicular to the boundary. 

The model used in treating this problem may be described as follows. Let 
a perfect crystal be cut along a plane S which is midway between the adjacent 
atomic planes S, and S, of the half-crystals A and B respectively. Let the 
meshes of S, and S, have square symmetry, the side of an elementary square 
being a. Cartesian coordinate axes are chosen in such a way that the plane 
z =0 is coincident with S and that the x and y axes are parallel to atomic rows 
in S, (and hence in Sp). For convenience we call these rows the x and the 
y rows respectively. Let B be twisted rigidly through an angle 3@ in one sense 
and A through 46 in the opposite sense. A force, analogous to force (1) of §2, 
now comes into play and is balanced by a force analogous to (11). If we do not 
allow force (i) to act, the meshes of S, and Sz will remain square as indicated 
by the dots and crosses in Figure 4(a), which represent the atomic positions in 
S, and S,, respectively. Let S be the plane of the paper, S, being below and S, 
above. If we now allow the force (i) to act, it causes a distortion which is 
particularly marked in the meshes of S, and S,;, but also extends some distance 
into A and B. The new atomic positions in S, and Sx are as indicated by the 
circles and plus signs respectively. We may now define u, v and w as in §2. 
We further define the x component of relative displacement of two atoms in 5S; 
and S, which were in alignment when 6 was zero to be the algebraic difference 
of their x components of displacement from the aligned position and denote it 
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by U. We may similarly define V as the y component of relative displacement. 
The distortion gives rise to a state of affairs analogous to that existing at the 
boundary with a difference of spacing in the two directions as in I. The 
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Figure 4(a). Diagram of the atomic configurations on either side of a twist boundary. The dots, 
circles, crosses and plus signs represent atomic positions. The dots are the positions on 
one side of the boundary and the crosses the positions on the other side before relaxation. 
During relaxation the dots and crosses become the nearest circles and plus signs respectively. 
The two perpendicular sequences of dislocations are most easily distinguished by viewing 
the figure at grazing angle. 


By 


Figure 4(b). Diagram showing only the sequence of dislocations which is parallel to the y axis. 
The broken and full lines represent the atomic rows on either side of the twist boundary, 
before and after relaxation respectively. i 


dislocations are, however, now of the screw type introduced by Burgers (1939). 
We note that the distance between dislocations is 


p=tacosec 46, Sen een 0) 
that V (or U) increases in every interval p along the « axis (or y axis) by 
b=asec 36, cereale SU) 
and that, as in I, the state of strain of A and B is periodic in x and y with period p. 
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It may be shown, by an argument similar to that of § 6, that the equations which 
we use to solve this problem are linear to the degree of approximation used and 
therefore allow superposition of the solutions corresponding to the two sequences 
of dislocations. It is therefore only necessary to treat one sequence, say the 
sequence parallel to the y axis, and for this purpose we fix our attention on the 
atomic rows which were originally parallel to the x axis, the corresponding picture 
being that of Figure 4(b). As we are considering only one sequence of dislocations, 
we take uw and U to be zero and V to be a function of « only. Further, v must 
vanish at points far from the interface, i.e. at = + 00. Since we are concerned 
with screw dislocations, we may also take w to be zero. This approximation, 
like that of §2, whereby the x component of relative displacement vanishes, also: 
neglects the change in the equilibrium value of the x component of separation 
where the atoms in S, and Sg are out of alignment and consequently the 
corresponding term in the interfacial energy. The z component of the force (i) 
is therefore assumed to vanish. We treat A and B as elastic isotropic continua, 
the conditions for the validity of this being the same as before. We must also. 
assume that the overall amplitude of the potential energy of an atom in Sz (or S,): 
with respect to the atoms in S, (or Sg) is independent of the angle of twist. In 
the continuum approximation, this implies that the overall amplitude of the surface 
density of potential energy at a point in S, (orS,) with respect to S, (or Sj) is 
independent of 6. If we further adopt the assumptions of §2 relating to the 
force (i), we may express (i) in terms of the shear stress, 


Pa Uea) = 2a )sim(2 Vib), 9 3 Se ea 3 (31) 
where Da iig COSA 0 ey Me ha eee (32) 


Equation (32) is the condition for the overall amplitude 72b/27? of the potential 
energy density, as obtained from (31), to be independent of 6. 


§8. THE GOVERNING EQUATION AND ITS SOLUTION 


The next step in the solution of this problem is to set up the equation governing 
the equilibrium of the forces (i) and (ii). "To do this we require expressions for 
the displacements and stresses which we obtain by solving the differential equations 
for the displacements in the elastic continua. In our problem u and w are zero,. 
and so these equations (Timoshenko 1934) reduce to Laplace’s equation in 
two dimensions, namely 02v/dx2+0?v/dz*=0. The appropriate solution must 
vanish at z= + co and on Sy, (and S,) at x= +4, and must have equal and 
opposite values at opposite points on S, and Sx. 

The most general. solution satisfying these conditions may be expressed. 
as a Fourier series : 


y & CU OSIM Wet hay seers. (33) 


where X, Z (and £) and the sign convention are the same as before. We find by 
differentiation 


dv pb = FnZ as 
Payee pee A,e*"sinnX, 
ov pee "or 
and Poy = 5, = + rome A eae" COS Gee SV Fete. (34) 
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“On S, and Sz wv and p,,, reduce to 


ae a 3 salah) 4 eee (35 a) 
2a n=1 1 
pope 5 
-and Pay = i 2 Asin 2Xi0 1, eee (35) 
n=1 
Using (35a), we may construct 
Vath oy is ea (36) 


Qn TT n=1 AW 


“the zero of V being taken at x= —4p and the first two terms having the same 
significance as in (11). 

The problem is now that of finding V from equations (31), (356) and (36). 
These equations are seen to be of the same forms as (1), (12) and (11) respectively, 
and the solution is, therefore, 


V =4b +(b/7) arctan {[(1+8-*)?+fB-4]tandX}.  ...... (37) 

The parameter 8 and the Fourier coefficient A, are similarly seen to be 
P=rub/pp = 27pm ‘tan 40sec 30 | 
cand A,,={(1 +6%)!—Byr=A” (say). | 


We may likewise determine the sums of the series in equations (33) and (34). 
We find 


Oats Me arctan ey: 
~ Qa 1— Ae*ZcosX}’ 
Pey = —(ubA/pR?)e*7sin X sete cares (39) 
ane Pny = + (ubA|pR2)e*4 {cosX—Ae*}, | 
where R®=1 4 A2%e* 4 2 Act cos X 


The expressions for the stresses and displacements associated with a single 
screw dislocation may be deduced by the method described in §3. They are 
found to be 

V =a +(a/7) arctan (x/Cy), | 


uv =(a/27) arctan {x/(€ + Co)}, | een (40) 
Day = —(wa2r)xr” | 
and Pray =(ual2r)(E + fo)/7°, J 
where P= ee (CA Ca)e | ands 6y sania. 


§9. THE INTERFACIAL ENERGY 
The interfacial energy 7 again consists of two parts as in $5. These may be 
calculated in the same way as before except that the strain energy density is now 
given by Vy =(p%,, +p2,)/4u. We must take the second sequence of dislocations 
into account. ‘This is a simple matter as the energies associated with the two 
sequences are additive to the degree of approximation used. We find 


7 =(pa/2n?){1 +B — (1 +6?) — Bln [26(1 + 8*)#-267]}.... (41) 
_A graph of 7(@) is shown in Figure 5, curve C. 
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Expression (41) reveals a shortcoming of the treatment, namely, that the inter- 
facial energy does not come back to zero when @ becomes 7/2, as it should on 
account of the square symmetry of the meshes in S, and S,. Instead it tends to 
a value j9@/27?, which will be the value of the potential energy density when the 
meshes in Sz and S, remain undistorted. At @=7/4, + is somewhat less than 
Mo@/2n°. A value of this order of magnitude is to be expected, for the average 
frequency with which the force (i) changes sign is now approximately the 
reciprocal of the interatomic distance, and the departure of the meshes from their 
undistorted forms will therefore be small. We may therefore take the formula (41) 


to be a good approximation when 0 is less than 7/4, and to be particularly good 
when @ is small. 


@ (degrees) 


Figure 5. Interfacial energy as a function of the angle 9 through which the two faces in contact 
are twisted with respect to each other. Energy in units of pa/27” per unit area, where pis 
the rigidity modulus and a the lattice parameter. (The curves are calculated for wy)=4.) 


Curve A. Elastic energy in the half-crystals (Ee/p). rns 
Curve B. Energy resulting from sinusoidal force at misfitting interface (Ea/p). 
Curve C. ‘Total energy of interface {r=(Ee+Ea)/p}: (41). 

Curve D. ‘Total energy according to the symmetrical formula (42). 

Curve E. Total energy according to the approximate formula (43). 


The fact that + does not return to zero at 7/2 is due to the approximation 
involved in expressing the force (i) in the form (31). ‘This approximation implies 
that we treat the atomic rows in S, and Sp, as though they were continuous 
(see Figure 4(b)) and the force between them as though it acted only between 
a point in a row of S, and the nearest point in a row of S,, and then only if the two 
rows were parallel when 6 was zero. Furthermore, the interaction may be des- 
cribed as a tendency to keep these points as close as possible to the zx plane 
(or zy plane) through the nearest point of complete alignment. The approximation 
is agood one when @ is small. The interaction between rows which were initially 
perpendicular, though negligible when @ is small, becomes greater than the 
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interaction between the rows which were initially parallel as soon as 6 exceeds 7/4- 
As a first approximation we take the interaction energy between a row in Sp and 
the two sets of rows in S, to be additive and construct the symmetrical formula 


7,=7(0)+7(4r—O0)—7(47) ees (42) 


for the interfacial energy. This gives a zero value at both @=0 and @=7/2. 
The graph of 7,() in Figure 5 (curve D) is seen to differ only slightly from the 
graph of 7(@) (curve C). 

The equivalent assumptions to those used by Shockley and Read (1949) 
in deriving equation (55)—see § 13 below—will lead in the case of a twist boundary 
to the formula 


72 (ual2n)\O(C—1n8),, 22) ee (43) 


Equation (43), like (55), is intended to be a good approximation when 6 is small; 
we may therefore evaluate C by comparison with (41). It is found to be 
C= In(euo/u27). The corresponding curve E is shown in Figure 5, and is seen 
to depart seriously from our curve for values of @ exceeding 5°. 


§10. AN APPLICATION TO ORIENTED OVERGROWTH 


Frank and the present author (1949) have developed a theory which explains 
the formation of completely oriented overgrowths on crystalline substrates 
when the difference in lattice spacing is smaller than a certain critical value. 
According to this theory, the first molecular layer that deposits deforms homo- 
geneously in such a way that its atoms are in the potential troughs of the substrate 
field. The second monolayer deforms in a similar manner, and so on, to form a 
completely oriented overgrowth. 

However, when the overgrowth attains a certain thickness it becomes unstable 
because of the increased strainenergy. Thisenergy may be released by the breaking 
away of the deposited film from the substrate. If there is no change of structure 
by recrystallization at this stage, the new state of affairs being that described in I, 
then, provided the overgrowth does not swing round, the initial orientation will 
be maintained. ‘The considerations of this paper may be used to show that if 
we assume the thick overgrowth has the same symmetry properties as its first 
molecular layer with respect to rotations about the normal to the contact plane, 
then it at least cannot swing round freely, since this is certainly accompanied 
by an increase in the energy. That this is so may be deduced as follows: the 
maximum value which the interfacial energy can attain in a rotation is probably 
not less than the maximum interfacial energy 7,,, ina pure twist boundary (Figure 5). 
‘The maximum increase in energy during a rotation is therefore probably not 
less than 7,,—7,, where 7, is the interfacial energy when the deposit has broken 
away from the substrate. Since the maximum value of + in I is less than in II, 
T,—71 1s always positive. The magnitude of 7, depends on the density of dis- 
locations, i.e. on the degree of misfit. Since the degree of misfit with which we 
are concerned is small—less than the critical value of 9% to 14°% below which 


the initial layer is oriented—7z, is small, and the increase in energy 7,,—7, 18 
therefore large. 
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III. Bounpbary pur To TILT 
See LE DaViO IDs L, 


The third problem is that of finding the stresses and energies associated with 
a symmetrical tilt boundary, i.e. the boundary between the halves of a crystal 
which are symmetrically rotated with respect to each other about a common 
axis in the boundary. 

Consider the case in which the undisturbed crystal lattice is of simple 
rectangular structure, the lattice parameters being a,, a, and a3, and in which 
the axis of tilt is parallel to a crystal axis (corresponding parameter a, say), the 
planes (a,, a) of the two portions making angles 14 with the boundary. We 
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Figure 6. Diagram showing three possible configurations of tilt boundaries. 
g : i 


restrict ourselves to the case a,=a,(=a), a, being unrestricted; the extension 
to the case in which a, and a, are different is quite simple. Let the boundary 
be the xy plane of a Cartesian coordinate system, the y axis being the axis of 
rotation. We denote the half-crystals by A and B. Asin I and II, two forces,. 
which we may call (i) and (1i), act on A and Dees bite force (11) may be defined 
as in §2, but the force (i) now acts essentially perpendicular to the boundary. 
The effect of the forces is again to reproducea configuration which may be described. 


42-2 
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as a sequence of dislocations, the dislocation lines being along the edges of the 
atomic planes terminating in the boundary. ‘Three possible configurations are 
shown in Figures 6 (a), (4) and (c). 

The form 6 (a) is the one shown in Burgers’ diagram (1939); 6 (5) and 6 (c) are 
alternative forms suggested by F. C. Frank (private communication). ‘The form 
(b) contains dislocations in their two possible symmetrical positions and (a) and (c) 
each contains dislocations in only one of the possible positions. Since by using 
(b) we obtain a precisely planar boundary, this should be a configuration of lower 
energy than (a) or (c) if the energies of the dislocations in their two possible 
positions are equal, as appears from the calculations of Nabarro (1947). The 
latter result is presumably a consequence of the approximations used. If the 
energies in these two positions differ, the configuration of lowest energy may 
well depend on 0, being that of Figure 6 (a) or (c) when @ is small, and (b) when 
6 is large. 

We define a conventional normal displacement W as follows. ‘The deformed 
half-crystals A and B, lying below and above the surface S respectively, may be 
considered to define systems of curvilinear coordinates in the two half-spaces. 
Standard mathematical methods of interpolation could be used, firstly to define 
continuous surfacés passing through all the atom centres on each principal 
crystallographic ‘plane’, and then to define all intervening surfaces of the 
curvilinear coordinate system. In these coordinate systems let Z, =0,1,2.... 
and Zi =0, 1,2... . denote those surfaces which, at a distance from the boundary 
S become atomic planes making angles +40 with S. Now, in each interval 
2p of x, where p is the distance between dislocations, the separation of continuous 
lattice surfaces across the boundary increases by 2b, while Z, +Z} decreases 

_by 2. The parameters p and 6 are given by the formulae (29) and (30) 
respectively. Let surfaces S, and Sp be described in the two half-crystals, 
parallel to the boundary surface S and each at a distance 4 from it. Let 
Z\(x) and Z}(«) be the coordinates of points in these two surfaces, then 
W = — b{Z} (x) +Z}(x)}. 

Between each successive pair of dislocations there is a value of x, which may 
be taken to be half-way between them and for which W/b is an integer. The 
distance between atoms in the two half-crystals is there the same as it would be 
in the undeformed crystal, and the normal force (i) should therefore be zero at 
such points: it changes sign at these points, and should vary linearly with W 
in their neighbourhood. _ It changes sign again at each dislocation, corresponding 
to a half-integral value of W/b. We assume the normal force to be proportional 
to sin (27 W/b), this being the simplest function which has the desired properties ; 
our knowledge of the inter-atomic forces is insufficient to justify use of any more 
elaborate function. We accordingly replace the half=crystals above S, and below 
S, by Hookian isotropic elastic continua which have suffered small displacements 
at these surfaces. "The accompanying tangential displacements (and corresponding 
forces), excluding those which are a consequence of the normal force, are 
disregarded. The error thus introduced is only serious when 6 is large: an 
approximate way of compensating for it is indicated below ($13). We also 
adopt the assumption that the overall amplitude of the surface potential energy 
density is independent of 8. 

__ As before, we express the sinusoidal force (i) in terms of a stress p,,(0, x). 
We fix the amplitude of the sine-function such that this function approximates 
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to Hooke’s law in the linear region half-way between dislocations. We further 
make use of the fact that p,., and p,, are equal at S where Pp: Vanishes according 
to our assumptions and of the equation for e,, which is analogous to (3). We find 


p.A0, x) ={a/n(1—20)}sin(2nW/b), Ss... . . . (44) 


where is given by equation (32). Certain aspects of the sinusoidal approximation 
are discussed in § 14. 


§12. THE GOVERNING EQUATION AND ITS SOLUTION 


‘The next steps in the solution of this problem are very much the same as those 
inI. ‘The proper stress function must now be chosen from (4) in accordance with 
the assumptions of §11, i.e. in such a way that the stresses vanish at z= + «, 
that p.,, vanishes at S, that w vanishes at S where x= +4, and that w has equal 
and opposite values at opposite points on S, and S,. The most general form 
of the stress function can be expressed as a Fourier series. If we use a convenient 
form for the Fourier coefficients, the mth term of the series is of the form 

a mis a (ltnZ)e*"@sinnX,  —...... (45) 
where X, Z (and ¢) and the sign convention are the same as before. By 
differentiation according to (2) we find 


ub 


ee = — A, (1 FZ )e*" sin nX, 

pb ay 
w= A, (LenZjet™sinnX + ...... (46) 
P= — hs A,(1 02) 


and Pia A,nZe'"4 cosnX. 


ub 
p(1—9) J 
If we now substitute from (45) into the equation analogous to (3), integrate with 
respect to ¢, let £ tend to zero and choose the integration constant such that w(x): 
vanishes at x= +4), we find 


We (DA If 2anysinnXs 9 ee ciel (47) 
By using this we may now construct 


Fi ee SO sins Ngee ae Ae o (48) 


27 TT n=1 1 


where the different terms have the same significance as before. From (45) 
we may also write down the Fourier representation of (44), namely 


Pa nia Seance: Wee are (49) 


1—c) n=1 


We must solve equations (44), (48) and (49). These are of the same forms as. 
- before and the solution is, therefore, 


W = 4b +(b/m) arctan {{(1+f-2)'+B-]tandX}. 2... (50) 
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The parameter f and the Fourier coefficient 4, are similarly seen to be 
B =zp(1 —2c) tan 46 sec $0/9(1 - 2) 
and A, ={(1 +B%)'—B}"=4" (say). | 


The nth terms of the Fourier series representing the stresses are given in (46). 
‘The displacements and the sums of these series can be obtained as described in § 4. 
We find 


Pow = —QOsin X{R?F [1 — A?eF™]}, . 
Pov = —QOsin X{R* + [1 — A?e**7]}, 
Pow =OZ{[1-+ A2e**4] cos X —2Ae*}, 
b Ae*t“sin X TAet ean X 
w= mae arctan ==>. us oe 
and We aber In (=a) + Tome [cos X— Ae*?]} 
PeRee (52) 


where Q=yubAe*?/p(1—c)R* and R?=1 +A?et?4 -2Ae*% cosX and where 
u and w are taken to be zero at the origin. 


The stresses and displacements associated with a single dislocation may be 
deduced from (52) by the method described in $4. We find 


Pre = SBE ene 246 Laane | 
Pez = = guint 2c(C + fy)}, 
Pow =9E{(E £ bo)? — 2°}, 


wg Lg 208 AAV SUES C0) ON ee 53 
"ar laied *\G) lar | me 
a x x 
w= se {arctan (; z 7) Un 6 are 

ae W = ha +(a/m) arctan (x/C), 


where g=pa/2r(1—o)r*, r2=x?+(C+0))? and f)=a(1—2c)py/4u9(1—<). 


§13. THE INTERFACIAL ENERGY 


The interfacial energy consists of two parts and is calculated as in §5. We find 
7 = {ig@/2n*(1 —20)}{1 +B — (1 +62)! — B In [28(1 +62)! -287]}. ...... (54) 


The formula (54), like (41), gives a non-vanishing value of 7 at O=7/2. We 
may also correct this error by constructing a symmetrical formula of the form (42) 
for the interfacial energy of the tilt boundary. 

; Shockley and Read (1949) obtained an estimate of the interfacial energy of 
a tilt boundary by assuming that the energy associated with an element of volume 
which is nearer to a dislocation than the distance between dislocations does not 
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vary with the tilt and that the strain energy outside these regions is negligible 


They find 
(Tr={pa/4r(1—o)}(B-In6), = ...... (55) 


where the quantity B, which they do not evaluate, depends on the energy in 
the region surrounding the dislocation line, where Hooke’s law becomes 
invalid. (55) is meant to be a good approximation when 6 is small, and may 
therefore be evaluated from (54) by comparison with (55). We find 
B=In {e(1—o)uo/2u(1—2c)}. The graph of (55), shown in Figure 7 (curve E), 
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Figure 7. Interfacial energy as a function of the angle of tilt 9. Energy in units of pa/2n*(1—20) 
per unit area, where a is the lattice parameter, » the rigidity modulus and o Poisson’s ratio. 
(The curves are calculated for c=0-3 and py=p.) 


Curve A. Elastic energy in the two half-crystals (Ee/p). ed ; 

Curve B. Energy resulting from sinusoidal force at misfitting interface. 
Curve C. Total energy of interface {r=(Ee+F)/p}; formula (54). 

Curve D. Symmetrical expression for total energy; formula (42). 

Curve E. Total energy according to the formula (55) of Shockley and Read. 


together with graphs of (54) (curve C) and (42) (curve D), is seen to be comparable 
with (54) and (42) over a wider range of @ than was found for the corresponding 
curves in II. 

For the values xy = 4 and o = 0:3, used in plotting the graphs, B is equal to 0-415, 
Shockley and Read found that when they take B equal to 0-5 their curve is in 
good agreement with the experimental observations of Dunn and Lionetti (1949). 
According to our formula B is equal to 0-5 when py =p and o =0-322. 


§14. DISCUSSION 

It is appropriate to consider now what are the advantages and disadvantages 

of the present treatment, particularly as compared with the calculations of Shockley 
and Read and of Burgers, and what is the general significance of the results. 

To begin with, the resolution of the boundary into a sequence of dislocations 

is a natural consequence of the present analysis. The lattice fault of I is a state 

of bad alignment between nearest atoms on either side of the boundary and is 

almost wholly contained in the distance in which U varies from ¢/4 to 3c/4.__ This 
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interval, which is given by (2p/7) arctan {(1+f-?)!'—B 1}, may be termed the 
width of a dislocation, and analogous definitions may be used for dislocation 
width in the cases II and III. The width, thus defined, as well as the energy 
density at all points, varies with the distance between dislocations. The treatment 
of Shockley and Read does not allow for this variation. When the distance 
between dislocations becomes large, the dislocation widths approximate to 
ap|tug(1 —c), ap/9 and a(1 — 2c)44/}492(1 —c) for dislocations of the types I, IL and 
IJI respectively. 

The width of a dislocation of the type III may be regarded as a measure of the 
rapidity with which the sinusoidal approximation allows for the atomic ‘planes’ 
which terminate in the boundary. In the case of a single dislocation, for example, 
the separation of atomic ‘planes’ on either side of S increases in a distance Cy 
by 4a, or in a distance a by a(2/7) arctan (a/2¢,), as can be seen from (53). The 
latter increase in separation and the distance ¢) respectively approximate to 5a/6 
and 0-35a when py=pu and c=0-3. The corresponding allowances are of the 
order which is to be expected and may therefore be regarded as partly justifying 
the use of the sinusoidal approximation. 

It is also to be counted an advantage of this treatment that singularities are. 
avoided, and the departure from Hooke’s law in regions of high strain is allowed 
for, even if only in an approximate manner. 

The question should be considered whether in the regions treated as continua 
there are strains too large for the validity of Hooke’s law. In the strain associated 
with the relative displacement of the two atoms in Sz on either side of a single 
dislocation may be taken as the maximum strain occurring in this case. From (25) 
it is seen that this is {u(4a)—u(—4a)}/a =arctan(a/2l9), which, for uy=p and 

.¢=0-3, is approximately 1/5. Analogous definitions may be used for the 
corresponding quantities in II and III which are then found to be approximately 
1/4 and 1/2-5 respectively. The strains of I and II probably are, and those 
of III certainly are, beyond the limit at which Hooke’s law becomes invalid. 

One necessary condition to permit our treating the half-crystals as elastic 
continua is, as given in §2, that the relative displacement of neighbouring atoms. 
in each half-crystal must change only slowly from atom to atom. The difference 
between the relative displacement of the two atoms considered in the previous. 
paragraph and the relative displacement between one of these and its opposite 
neighbour is approximately the maximum variation that exists in case I. From 
(25) this difference is seen to be u(3a/2) — 2u($a) +u( — 4a), which is approximately 
—a/9. ‘The corresponding differences in II and III are about —a/7 and —a/2-6. 
The continuum approximation is therefore reasonable everywhere in the half- 
crystals of I and II, but is certainly bad for III inthe neighbourhood of a dislocation. 

The assumption that we may treat the surfaces of the half-crystals as continuous. 
may be regarded as a separate assumption. It allows us to express the force (1) 
in terms of stresses and requires U, V and W to vary slowly from atom to atom, 
1.€. 1t requires the width of a dislocation to be large. The dislocation widths, 
which are respectively 1-4a, a and 0-35a for uy =u and o=0-3, are rather small; 
this is the main objection to the treatment and other treatments in which the 
Peierls assumptions are used. 

__ A further objection is the uncertainty as to whether the sinusoidal law of force 
is good enough in the region where it does not approximate to Hooke’s law, 
i.e. in the small region around a dislocation where, for example, sin 27% cannot. 
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be replaced by either 2m or —27(1—7). If the approximation is bad it will 
introduce appreciable errors into the interfacial energy and in particular into 
that part associated with the potential energy of the atoms which lie in this region 
of high potential energy. ‘The width of a dislocation may be regarded as a relative 
measure of the extent of this region. The error will be large when the dislocations 
are close together since the fraction of the total number of atoms at the boundary 
which are in the regions of high potential energy becomes large. A consideration 
of the width alone would predict that the smallest error is obtained in case III. 
However, it is very likely that the sinusoidal approximation is the worst for this 
case. 

We finally conclude that a single mathematical treatment enables us to calculate 
the stresses and energies associated with three different types of boundaries 
and that these results are in all probability good approximations when the 
distance p between dislocations is large compared with the interatomic distance a. 
The errors, depending particularly on the sinusoidal approximation, may become 
appreciable when p is comparable with @ and are certainly greatest in III. The 
outstanding feature of the interfacial energy is its rapid initial increase with the 
density of dislocations. It is unfortunate that, except for the work of Dunn and 
Lionetti (1949), no experimental observations with which the present results 
could be compared have been made. In view of the importance of boundaries 
of misfit in questions concerning oriented overgrowths or the physics of polycrystals 
further experimental measurements are much to be desired. 
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The Theory of the Propagation of First and Second Sound 
in Helium II.—Energy Theorems and Irreversible Processes 
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ABSTRACT. ‘The first part of this paper contains a discussion of the conditions deter- 
mining whether sound waves in helium II are propagated isothermally or adiabatically. In 
the remainder of this part an expression is found for the energy flux from a closed region, and 
the impedance concept is applied to waves of second sound. 

The second part of the paper is devoted to the study of the effects of irreversible processes 
on the propagation of waves of first and second sound. ‘The irreversible processes discussed 
are those due to viscosity and thermal conduction, both in an extended fluid and in a narrow 
tube. 


§1) IN DRODUCLION 

HE work given here is a continuation of that of an earlier paper of similar 
| title (Dingle 1948, referred to as I). In the present paper, the emphasis 
is on a unified theory of the propagation of first and second sound in helium. 
The interactions between the two types of sound are included in the theory. 
The first part of the paper contains a discussion of the thermodynamical nature 
-of the fluctuations due to sound in helium IJ, and some considerations on the. 
flow of energy; these will be required in the second part of the paper. ‘The 
proof of the energy theorem of $3 is based on some results of a paper on the 

hydrodynamics of helium II (Dingle 1949). 

The irreversible processes discussed in the second part are those due to 
well-established effects—viscosity and thermal conduction. In efforts to find 
expressions for the thermal convection, other irreversible processes have been 
postulated. Meyer and Band (1947, 1948) have assumed the existence of a time 
of relaxation inversely proportional to the square of the heat flow, and Gorter 
and Mellink (1949) have assumed the existence of a mutual force of friction 
between the normal fluid and the superfluid proportional to the cube of their 
relative velocity. In the present paper we shall not discuss these possible causes 
of irreversibility, since in the absence of theoretical justification they must be 
regarded as empirical formulae describing the additive result of many irreversible 


‘effects. In any case both these postulated effects will vanish for sufficiently 
small amplitudes. 


PART I—GENERAL CONSIDERATIONS AND ENERGY 
THEOREMS 
§2. THE THERMODYNAMICAL NATURE OF THE SOUND 
ERUC IUPAC @iNiS 
In normal fluids the fluctuations in density and-pressure due to the passage 
-of sound waves are adiabatic, because of the small thermal conductivity of the 
medium. However, helium II transports heat very easily, and we must inquire 
whether the fluctuations may still be considered as adiabatic.* The situation is by 


* This problem was discussed by G 0] x i i 
2 Se ee y Groenewold (1939) before the nature of heat flow in helium II 
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no means the same as that which would exist in the hypothetical case of a liquid 
with a very high heat conductivity of the familiar type (i.e. a heat transfer not 
involving movement of the fluid), but otherwise normal. For in helium II the 
transfer of heat is itself due to a wave motion—second sound. Quantitatively, 
the problem may be investigated as follows. Neglecting the corrections due to 
viscous effects, we obtain the following equation for the sound velocities (cf. (A3) of 
the Appendix, also Landau 1941) : 


ee ana) 

The two solutions are given by 
2u* =tgg? +Uy 3” + [(thy5 — Wag”)? +4aq"(Uis? —Mar)}*, se eee (2.2) 
MCC m eto. tp sep. i c==(Opiop)e. and . Up? =(OP/Op) qs. 1a’ seas (2.3) 


‘The difference between the adiabatic and isothermal compressibilities is small, 
so that u,,2~u,,?._ Expanding in terms of their difference, we obtain for the 


velocity of first sound, given by the solution corresponding to the positive sign 
of (2.2), 


ere) jee 6. 4 


eee 2 
Ugg (Uy s" — Uy 7") . 
20S TE ty? + 


Uy 5°(Uy 5" om U7") (2 4) 
Uys" — Ugg” Uys” — Ugg” 


Uy? = Uy5" + 

Thus first sound in helium II is propagated neither quite isothermally nor 
adiabatically, the precise nature of the fluctuations depending on the ratio of 
the velocity of first sound to that of second sound. In the temperature region 
so far investigated experimentally, (T>1°K.), (te)max~M1/13, so that ug9? <uy,¢7. 
In this case, we see from the first of the alternative formulae (2.4) that u, =1,, 
and thus first sound will be propagated adiabatically, as in other fluids. 

If, on the other hand, there exists a temperature region in which the velocity 
‘of second sound is considerably greater than that of first sound, as is unlikely, 
the second relation of (2.4) shows that the first sound fluctuations would be 
isothermal. 

The velocity of second sound is given by the solution corresponding to the 
negative sign in (2.2), 1.e. by 


2 2 ee ee 
Uys —Uyp Ujyr — Uso 
te ae 4 lon? 5 | cea Cr (2.5) 
Uis — 420 Us" — Uso 
The deviations of the actual velocities from those calculated on the assumption 


of isothermal or adiabatic propagation may be calculated thermodynamically, 
the ratio of u,, to u,, being given by 


& 2 
Re he homer oh este (2.6) 
Uy 7" Gy Ga 
where « is the coefficient of thermal expansion. Thus ignoring terms of the 
second order, 
2 2 2 4 shi 2 
Tipe ace aH es Uyp Ol 
ie 20 “1S =e + ——* =. Zl, 
estas [+ ] se [1 typ — Ugg?” a oe 


2 
Uy 5° — Ug)” Cy 


ue Ta? 
ioe etlone E a, eae oR FP a” PB aurtoxtarrie) (or o1-0) (2.8) 


and 
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Hence in the temperature region so far investigated, and probably at all 
temperatures, the velocity of first sound is greater than that calculated from either 
the adiabatic or the isothermal compressibility, and the velocity of second sound 
is less than that calculated from the usual formula u.2=T\S%p,/Cp,. It was 
shown by Stokes (cf. Rayleigh 1896) that in ordinary liquids sound waves would 
suffer extreme attenuation if propagated under conditions which were not very 
nearly adiabatic, or very nearly isothermal. It is interesting to note that this 
conclusion does not apply to the propagation of sound in helium II; in fact we 
have shown above that there is no attenuation in spite of the fact that the waves 
are propagated neither isothermally nor adiabatically. This difference in 
behaviour exists because the transfer of heat in helium II is primarily due not 
to conduction but to a convection current, which in this liquid is itself propagated 
as a wave motion. 

The actual magnitudes of the corrections to the velocities are negligible. 
For instance, at 1:5°K. we may take «~—3 x10-8deg+, C, ~ 0-28 cal.gm™, 
u,~2:5 x 10cm. sec~1, giving Tx2u,2/C, ~7 x 10-4, showing that the difference 
between the two compressibilities is minute. 


§3. A. THEOREM ON DHE hEOWVORSENERGY. 


In the theory of electricity there is a well-known theorem that the energy 
flowing per second from a region enclosed by a surface A is given by the surface 


integral | (c/47)ExH .dA: the integrand is the Poynting vector. In this section 
A 


we find an analogue to this vector representing the energy flux due to waves of 
first and second sound. 

As was shown in a recent paper (Dingle 1949)+ simple expressions may be 
found for the derivatives of the potential energy of unit mass of helium. With 
W as the potential energy per unit mass, these may be written, (cf. (6) and (7), 
of loc. cit.), 


(@Widc\s=p' lp? a) a (3.1), 


(QW/eS),=T', | ae ee (3.2) 


and 


where p’ and 7” are the fluctuations in pressure and temperature respectively. 
Next we introduce the coordinates x and y defined by the relations 


[ROTO OS (3 .3) 
and 


p Sy =PnpsS *(Xp a Xs), OF OG D1 G (3.4) 


where x, and x, are the displacements of the normal and superfluid respectively, 
S the total entropy per unit mass and S* the difference in entropies of the two 
fluids. The coordinate x gives the position of the centre of gravity of the fluids 
and y is a coordinate proportional to their relative displacement. In this system, 
the conservation equations may be written 


Opt + pV.x=0) | oi ere (3.5) 


and 
OS/OLSV.y =0a) |, ae (3.6), 


} For the remainder of this section ‘loc. cit.’ refers to this paper. 
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Thus the rate of change of the potential energy per unit volume is 


aa (0%) 208 
ae Op /spot’ op at 
by (3.1), (3.2), (3.5) and (3.6). 

The rate of change of the density of kinetic energy is 


= SV. w pT Vy ake cs: (3.7) 


7) ; : 
Ran reyiene 3(PnPn? +p.U5”) = vg (ORs » Un + Pgvs - Us) 
ot 


v | 
= = s(ROle pte: 7 SEVIS on): idee: (3.8) 


(by (17) and (18) of loc. cit.), 
=x.Vp+pSy.VT 


by (3.3) and (3.4) of the present paper. Hence the rate of change of energy 
in a region of volume 7 enclosed by a surface A is given by ; 


-- (total energy) = | lev .x+pST'V.y) +(x.Vp+pSy.VT)] dr 


=| V.(p'a +pST'y) dr = | (p'x +pST'y) dA. 
A 


A 


Thus the analogue of the Poynting vector of electromagnetism is the vector 
(p'x +pST'y). The physical meaning of this result is clear: p’x is the rate at 
which ‘ mechanical’ work is done in moving the centre of gravity of the two fluids 
against a pressure excess, whilst pST’y is the rate at which ‘thermal’ work is 
done in moving entropy (relative to the centre of gravity) against a temperature 
eXcess. 


$4. THE CONCEPT OF IMPEDANCE APPLIED TO SECOND SOUND 


In I, it was shown that the rate of flow of energy per unit area is, for a plane 
wave of second sound, 

1pnp(Vnt)*t, _ 1 ppsS*(T)” _ 1 puaC( Ti) = 115708), soe (4.1) 

2 Ps 2 Pnttg 2 ie i. 


where V,,° and T,,° specify respectively the peak values of the velocity of the 
normal fluid and the temperature excess. ¢$°=pV,,°S is the peak value of the 
flux of entropy, and H® =pV,,°S/T is the peak value of the flux of heat. Hence 


Flux of energy of second sound =} (peak fluctuation in temperature) 
x (peak flux-ofentropy) §...... (4.2) 
=1 (peak fluctuation in temperature/actual temperature) 
~(pealsheat tix)" e272. (4.3) 
Comparing (4.2) and (4.3) with the well-known expression for the electrical 
energy flux, 1.e. 


Flux of electrical energy =4(peak voltage) x (peak current density) 
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we see that for second sound there are two convenient ways of defining the 
impedance; either 

Fluctuation in temperature 

EEO Eee eaten 


Z =Impedance = Actual temperature 


giving 
17° T\2 1 (HH? 
Pasay & [rn yf pices) Oey 
Flux of energy of second sound = IF H 12 ( 7) DEAE (4.4) 
or 
Z =Impedance = Fluctuation in temperature/flux of entropy 
giving 


Flux of energy of second sound =47T°d® =4(T°)?/Z=3Z(6°)?. ...... (4.5) 


These results may also be deduced from (3.10) of the preceding section. 
The energy flux vector (p’x +pST’y) may be divided into a ‘mechanical’ energy 
flux vector 


DP L=P (PaV ap g)/P= = a ee ee (4.6) 
and a ‘thermal’ energy flux vector 
pal y=pSl (VaaV..). a ee eee (427) 


In first sound the ‘thermal’ energy flux is negligible, and hence the flux of energy 
is p'V =4p°V°, if the vibrations are simple harmonic. In second sound the 
‘mechanical’ energy flux vector is negligible, and hence the flux of energy is. 
pST'(V,,—V,)=pST'V,, since j =p,V,+p.V.=0; thus for simple harmonic 
waves of second sound the energy flux is $oST°V,°, in agreement with (4.1). 

Of the two equivalent definitions of the impedance for second sound 
. propagation, the second brings out clearly the nature of the work done—i.e. work 
done in transporting entropy through differences in temperature. Corresponding 
to ordinary mechanical power given by the product of force and fluid velocity, 
in second sound we have the thermal power given by the product of the temperature 
fluctuation and the entropy flux. 


PART II.—IRREVERSIBLE PROCESSES 
§5. THE VISCOUS ATTENUATION OF PLANE WAVES OF SOUND 
IN AN EXTENDED MEDIUM 


We assume that normal fluid has viscosity 4, and that the superfluid is non- 
viscous*. ‘Then, as has been shown in a previous paper (Dingle 1949) the 
equations to be satisfied are 


DF Vp— (VV, +4VVV)=0, eee (5.1) 
ae Vj =O a a ee (5.2) 
ee + pSV Vo, | ea (5.3) 
- -svrs 3? a0; 225. ie eee (5.4) 


* The supposition of a non-zero value of th 1 1 i i i 
e superfluid viscosity would be incompatibl 
with the assumption that the superfluid has no entropy. : : 
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These are respectively (28), (1), (2) and (4) of the paper quoted above, with 
V,, written for x,, the velocity of the normal fluid, and V, written for the 


velocity of the superfluid. By definition we have also p =p, +p, and j =p,V,, +p,V.. 
Representing fluctuations by symbols with a dash, (5.1)-(5.4) reduce to 


Vila —ep)atpVe—p /u=0,~ panes (5.5) 
TESS Oe | oes be a or (5.6) 

O] PomO Wes af O81 | a ee eae (5.7) 
UVert Se ep pi, 8 cae tien (5.8) 


where w is the velocity of propagation and ¢ =4iw/3u2p, w being the pulsatance. 
Eliminating V,, and V, from these equations, we obtain 


iat) freer pouSs bt @. ae. (5.9) 
and WES (pa Py Sime EU Sp laip. Ai— wea ue (5.10) 


We shall obtain approximate solutions of (5.9) and (5.10) by neglecting the 
right-hand terms, these giving only the coupling between first and second sound. 
This procedure will be justified in the Appendix. To this approximation we 
may also neglect the coefficient of thermal expansion, so that it is legitimate to 
take p’ =p'dp/op and S’ =T’0S/0T =CT'/T, where C is the specific heat and 
T the actual temperature of the fluid. Thus 


EN €)—- OP CP Ola | ee ae (5. leh) 
2 
and us ( -«) CANN» ae AON (5.12) 
4f Ps 

Since eu? =4i1.w/3p, we obtain as our two solutions 
UP =OpiOpsiual(Sp- ae ve ee (S213) 

TS*p, . 4iuw ps 
a= = +4 a ie Wa he 5.14 
ag ea Cen - 3p Pn ( ) 


Thus the real parts of u, and uw, are altered only in the second order of 
approximation. The distance in which the amplitude falls to 1/e of its 
initial value is given by 1/I(w/u,)=3u,3p/2u? for first sound, and by 
1/I(w/uy) = (3u22p/2u")(p,/p.) for second sound. Hence the attenuation of 
first sound due to this cause will not vary appreciably with temperature, since 
u, is approximately constant, but that of second sound will alter considerably 
with temperature, the attenuation being very large near the A-point, where 
u,3/p,~0. 

A rough calculation of the magnitude of the distance in which the amplitude 
falls to 1/e of its initial value may be made for a specific temperature. For 
example, at 1:5°K. we may take u, =2°5 x 104 cm. sec"; u,=2 x 10° cm. sec; 
p=0-145 gm. cm, p,/p, =7-4 and p =10-° poise, obtaining for the characteristic 
distance a value 8-8 x 10!7/f2cm. for first sound propagation, and 6-1 x 101°/f? cm. 
for second sound: f is the frequency in cycles per second. At this temperature 
the attenuation is important only at extremely high frequencies. Near the 
A-point, however, the attenuation of second sound will be serious for lower 
frequencies also. 
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§6. THE VISCOUS ATTENUATION OF FIRST AND SECOND SOUND 
IN A THIN TUBE* 


We assume that the tube is sufficiently wide for the ‘shell’ of fluid near the 
walls, which is affected by the boundary condition V, =0 at the surface, to be 
considered as thin. Then the forces acting on the bulk liquid are those due to 
the pressure gradient and to the force transmitted by this shell of fluid near the 
boundary. The latter force may easily be calculated if we consider the walls 
_as effectively flat; this approximation, and also that of assuming the ‘shell’ to be 
thin, may be justified @ posteriori. If the sound wave is taken to be travelling 
along the x axis, we have 
pO V,,/0t =|e ad? V,,/dx?, 
since the viscous force acts only on the normal fluid. Taking d/dt =u, 
i d? V jax aa 1p, Vi, /p. 
The solution satisfying the boundary condition V, =0 at x =0 is 


V7, =V,%1—exp{—x(1+iwpy/2w)M vee (6.1) 
where V,° is the velocity of the normal fluid well away from the shell of liquid 
-dragging at the perimeter of the tube. _ If this shell is thin enough to be considered 
_as located at the plane x =0, the force transmitted by it is, per unit area, 


d : 
be ( 2) =pm(1+2)V,(wpn/2p)i. we ee (6.2) 
xX / x=0 


Thus the viscous force per unit volume is 


P[w 3 . 
a (3H) a +)r, 
where P/A is the ratio of perimeter to cross section, and V,, has been written 
for V,°. This force opposes the acceleration. Writing 


ce a (5) 'e- ee, 2 a eee (6.3) 


‘the equation of motion of the bulk liquid reduces to (5.5) with e’ replacing « of §5. 
The elimination of the variables proceeds as before until the equivalents of (5.11) 


and (5.12) are obtained. However, unlike c, «’ does not involve u, and the solutions 
are now given by 


u?=(1=<«)=cpleps = a eee (6.4) 
-and Us* = (1 Se’ palo.) 7 S7ps/ Co, eee oe (6.5) 
‘Thus the real part of the velocity of first sound is altered by a factor of approximately 
a Lol (5) 
2p A \ 2 
.and that of second sound by a factort 
aS pee Gal 
2pnp A \2w ) * 


oaseal steady flow of helium II through thin tubes has been discussed by F. London and Zilsel 


+t The results for second sound which were given in I are erroneous, as the increased coupling 
between the sounds due to the viscous effects was not taken into account. ‘The equation 

: 4 Mnyatats —€'pV y+ u*p’—p’=0, 
obtained by elimination of 7 between (5.5) and (5.6), appeared to show that the effects of viscosity 


“(contained in e’) are associated with density and pressure fluctuations, i.e. with first sound, and thus 
with second sound only through thermal expansion. 
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For first sound, the distance in which the amplitude falls to 1/e of its initial value 
is 1/I(w@/u,) =1,pr(2/p, uw)! for a circular tube of radius r. The corresponding 
distance for second sound is (ugpr/p,)(2p,/um)'. The attenuation of second 
sound due to this cause will be serious near the A-point, where wu, ~0, and at very 
low temperatures, where p,~0. With the values given in §5, the distance in 
which the amplitude falls to 1/e of its initial value is, at 1-5°K., 1-6 x 10%rf-! cm. 
for first sound, and 5-3 x 1047f-} cm. for second sound; r is the radius of the tube 
in centimetres and f the frequency in cycles per second. 


§7. THE COUPLING COEFFICIENTS BETWEEN FIRST AND 
SECOND SOUND 
We have seen in §5 that the coupling between the sounds is caused by at least 
two different factors: the existence of thermal expansion (which means that the 
temperature fluctuations of second sound will be accompanied by slight pressure 
and density changes), and the viscosity uw of the normal fluid. In this section 
we calculate the magnitude of this coupling by a method due to Lifshitz (1944). 
Introducing parameters a, b and c defined by the relations V,, =aV,, p’ =bV.,, 
T’ =cV,, we have from (5.8) 
Die SC oe ee oe yee oe ne (71) 
and by (5.5) 
DU peep EPC = Upurea (7-2) 


Finally by (A1) of the Appendix, 


1G), w00-9 1-43) ] ve[o-o G8), -*268),]- 


Brean (7.3) 


Since (dp/dT),, and (0S/ép),= —(eV/0T), are proportional to the coefficient 
of thermal expansion, the products of these quantities with « may be neglected. 
Noting also that u,-? =(0p/0p)7(1 —«), we have 


b(u2/u,2—1)—pure(ateC/ST)=0. see, (7.4) 


Thus by (7.1) 
pe uwp(a +eC/ST) 


In a similar way, we obtain the following values for a, 6 and c: 


O,plyeus ean (14 So,Tu,? e o,fu,? 


Aenten alae Tg Ba BT Ries ye 
v Pn Spps(y2 — 4g”) pC(u,? — uy”) C(u,? — uy?) 
A Apert (7.6) 
Be ete 2P\ PR Uggs poe Oop ae a 
ae ae Pu Pn Spn (u,” Ree U9”) : : S(u? 2. Uy”) (7.7) 
U2 (4, _ % Uy thes 
es ae S (u,? — uy?) )’ J 


where €, =4piw/3puy?, €2 =4ptw/3puy?, op =% + Ce, TS, anda, =a +Ce/TS. 
At 1-5°K., «, and «, may be calculated using the values of §5, giving 

<,=9:2 x 10-ifande, =1:5 x 102%. With C =0-28 cal.gm-t and S =0-06cal.gm7! 

the additions to be made to the thermal expansion coefficient «~—3 x 10-3 deg"! 
PROC. PHYS. SOC. LXIII, 6—A 43 
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to take account of viscosity are Ce,/TS =2:8 x 10-if and Ce,/TS =4-7 x 10- Mf 
for first and second sound respectively. 

The additional coupling due to viscosity is negligible except for second sound 
waves of very high frequency, or at temperatures close to the A-point, where ¢«, 


becomes very large since u,~0. 
When the liquid is contained in a narrow tube, the additional coupling due 


to viscosity may be calculated merely by replacing «, and e, by e’, where 


(a 2 f pnt\* on 
€ == (5%) (@—1). 


At 1:5°K., «’=135(i—1)r1f-?_ and Ce’/T.S =42(i—1)r-1f-. Thus when the 
liquid is contained in a tube, the coupling between the sounds may be much 
greater than that calculated from thermal expansion alone, especially if the tube 
is narrow and the frequency low. 


§8. THE -RATIO OF THE AMPLITUDES OF THE ELRST AND ESECOND 
SOUND PRODUCED BY DIFFERENT PROCESSES* 


(a) A surface fluctuating in temperature. 

Let A, and A, be the amplitudes of the velocities of the superfluid in first 
and second sound respectively, these being produced by a source consisting of 
a plate undergoing fluctuations in temperature of amplitude 7’. The boundary 
conditions to be satisfied are: (i) The sum of the temperature fluctuations due 
to the two sounds must equal the temperature fluctuation TJ’. Since 

T=cV,, G Ay AC Ag = eee (8.1) 
. (ii) There must be no net flow of fluid at the surface; i.e. j=p,V,+ pV, =0- 
Since V,,=aV, this gives 
(ps ** @p,) Ay + (es Gppn) Ao = 0 5 See (8.2) 
From (8.2) we obtain the relation 


T;' ere A, =, Co Pe ores Co (ur =u")? ~ € 


TE, C1 pst+@2py, Toyo, uy>ug To yogt uy” 
At 1:5°K., taking o,~o,~a~—3x10-° deg}, T,'/T,’ ~1-38 x 104. 


(6) A surface vibrating in its own plane. 


Let V be the velocity of the plane. Then, since only the normal fluid drags, 
our boundary conditions are 


(i) 4,4, +a,A4,=V and (ii) A, +A,=0. 
Hence the relations 


To) Cp An 1G _ UgC(uy?—uy") Citta 


(Pe C8, Ci Soild, ao) Su 
At 1:5°K., taking o,~a, T,'/T,’ ~0-83 x 102. 


These numerical values of T,'/T,' are minimum estimates, as the increased 
coupling due to viscosity has not been taken into account. This additional 
coupling decreases the efficiency of these methods of second sound production. 


* The method of this section is also due to Lifshitz (1944). 
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§9. THE RATIO OF THERMAL TO MECHANICAL ENERGY IN SOUND 


If we consider a volume enclosed by a spherical surface of unit area, the rate 
of loss of energy is given by the product of the velocity of the sound concerned 
and the energy density within the sphere, i.e. by u(p,V,2+p,V2). Hence by (4.6), 
the proportion of ‘mechanical’ energy emerging from the sphere is, in the 
notation of §7, 


P' (PWV n + Ps V,) 5 b(pya + ps) 
UP(PrYVin> +pV2) — up(pna? +p)” 


For second sound this is equal to 


2 2, 2 2 
Oo Pn Uy"Us (9 2) 
Ss : oe Si Ue Pi Ph roertivORee 2) 
8 al 2 


The proportion of ‘thermal’ energy is, by (4.7), 


pST(Vn—V,) _ p,Se(a—1) 
U( py Ve ay O8 Vo) U(py,a” =F Ps) ; 


For first sound this is equal to 


a2 (See y et (0.4) 
Speen," 

When o,~o, (there is exact equality if the only coupling taken into account 
is that due to thermal expansion), the proportion of ‘thermal’ energy in first 
sound is almost identical with the proportion of ‘mechanical’ energy in second 
sound. At 1-5°k. this proportion is 3 x 10~® if only thermal expansion coupling 
is taken into account. 


10. DHE AT LENUATLON- OF SOUND DUE TO THERMAL CONDUCTION 
INV AN EXGLENDED MEDIUM 
We shall assume that in addition to the heat flow caused by convection there 
is a heat flow due to conduction which is directly proportional to the temperature 
gradient, with a factor of proportionality K per unit cross section per unit time. 
Considering unit volume, the rate of loss of entropy is —0(pS)/dt. This 
must be equated to the sum of the loss pSV.V,, due to convection (assuming that 
only the normal fluid carries entropy), and the loss due to conduction 


T-aQ/at = — KTVT", 


where T is the mean temperature of the liquid. Hence 


O(pS)/0f pSViV,—KIOAVT’ =O ee. (10.1) 

The equations describing the hydrodynamics of helium II (cf. §5) are now 
DAV et Vie =p i =O) i acta (10.2) 

Up ORV ee DeV —Oe  n a tilly (10.3) 
p'S+pS’—pSV,/u+iKkel'/Tw=0,  ...... (10.4) 

UV Sle pip =U.e5 Sere Gash, (10.5) 


* The elementary theory of this effect for second sound alone was given in I. §5. 


43-2 
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By the methods of the Appendix, we obtain 


oe efi (Rah ] [ae ee - 22 | (2) =U 
‘ees Ee Gs pce a: (2 S : Cy Pn pCy Op ely 


eee (10.6) 
Writing 
yg? =p.TS?|p,Cy,  4=iKelpC,, —4s?=(p/Op)s» tax” = (OD/Op) os 
Akasa (10.7) 
we have 
5 2(u29” — 7)(t1s* — Ua") 
2u? = Ugg’ — ne Uys" + Es Ugg? +79 + Sai cea 
Rera Re (10.8) 
The velocity of first sound is obtained by taking the positive sign in (10.8), thus 
Cee y Se. gee te (10.9) 
uy, Us 2(uys? —u,? +7)Uys" 


‘Since u12 — 47? =u,7*Ta?/C,, the distance in which the amplitude of first sound 
falls to 1/e of its initial value is 1/I(w/u,)=2pCyu,/Kw*Ta?. 

The velocity of second sound is obtained by taking the negative sign in (10.8), 
thus 
(w3” = 1) (15° = Uy7") (10.10) 


Oe Sat eet 
Us — Uo 7 (u, 2 = ue Zi 7) 95 2 8 A! ake ore: 


in which the third term is negligible in comparison with the second. ‘The 
distance in which the amplitude falls to 1/e of its initial value is therefore given 
by 2pC,ug3/Kw?. 

At 1:5°K. we may insert the following values in the formulae: 
p=0:145 gm.cm-, C,=0-28 cal.gm™, K=6x10-°cal.cm-*deg + (the value for 
helium I, where the conduction is not masked by convection), « = —3 x 10-3, 
u, =2°5 x 10*cm.sec}, wy=2 x 10? em:sec-t: Then we obtain for the character— 
istic distance the value 3 x 10!*f-2cm. for first sound, and 2-73 x 10"f-2cm. for 
‘second sound, where f is the frequency in cycles per second. 

The attenuation of first sound due to this cause is therefore negligible except 
at very high frequencies, but that of second sound will be serious near the A-point, 
where u,~0, and probably also near the absolute zero, where C,~0 with u, 
probably remaining finite. 


$11. THE ATTENUATION OF SECOND SOUND IN A THIN TUBE DUE 
TO THE THERMAL CONDUCTION OF THE WALLS 

The loss due to the thermal conduction of the walls of the tube may easily be 
estimated with the help of the impedance concept of §4. Two impedances are 
involved: Z),, that of the liquid, which will give the power travelling down the 
tube, and Z,,,, the impedance of the tube walls, which will give the power flowing 
out sideways, the power lost by the liquid. 

(1) The impedance of the liquid. From the work of § 4, Ayiq = T\iq/TH, where 
ff is the heat flux. Neglecting pressure and density fluctuations, Hepsi. 


Theory of Propagation of First and Second Sound in Helium IL 649 


whence H=p,pTSV,/p,=p,pTS?T" |p,tty (by (3:8) = pCT a, sinceu,?= 7S2p./Cp, - 
(This calculation is really unnecessary, since the result is physically obvious, 
pCT’ being the excess heat per unit volume and uy being the velocity with which 
It 1s transported (cf. Osborne 1948).) Thus 


PAs: 1 oe ROT Ne i cee ae ae are (11.1) 
The flux of energy travelling down the tube is therefore 
UT ZqcdpTOupA(TeV, wae (11.2) 


MaxX 
where 4 is the cross-sectional area of the tube. 


(ii) The impedance of the walls of the tube. The equation of heat conduction 
CB fe Kei 27 
“adi = em, V I. oe: eb 8) 'e! 6-6 (13) 
has the solution 
i a dha exp[ —a(1 +t) LPC [2K 5o}4]. a ele (4) 


Here we have neglected the curvature of the walls of the tube. The heat 
flux H is given by the relation H = —K,,(0T"/@z),-9=BT’,9x(1 +i), where 
B — (F@P. 1 Coo Keo)? Hence 

Zl EL = AB (ete ae. oS aca shee (11.5) 
Taking P as the perimeter of the tube, the energy flux into the walls is, per unit 
length, 


sol 


3(Tinax)”P/Zoo1 = 3BT(1 +2)( Trax)? P. 


max 


The flux of energy lost to the walls is given by the real part of this expression, 


eee Pio) bk re Meee Sak a PRS en ae Be (11.6) 
(11) The attenuation. By (11.2) and (11.6), 
Loss of energy per unit length ge ofl Ne (11.7)? 
Energv travelling downtube  pCu,A\Tiax/ © 977” ‘ 
Here 77...) Ty, is the ratio of the maximum temperature fluctuation in the liquid 


in the neighbourhood of the walls of the tube to that in the bulk liquid. Peshkov 
(1948) has assumed that these temperatures are thesame. ‘This simple assumption 
is probably false, as it neglects the effects of viscous forces, which will cause the 
amplitude of second sound vibrations to die out near the walls—the flow of 
heat in second sound is proportional to V,, the velocity of the normal fluid, and 
this must vanish at a surface if there is a true viscosity. It is possible, however, 
that ‘slip’ between the liquid and the walls may not be negligible; in this case: 
the value of (T¥,,./Tiax) Will lie between zero and unity. 


$12. THE ATTENUATION IN A THIN TUBE DUE TO THE COMBINED: 
EFFECT OF VISCOUS AND CONDUCTION LOSSES IN THE LIQUID 
In the last section we have shown that the temperature fluctuations of second 
sound die away near the walls of the tube, since V,, falls off in magnitude owing to 
viscous effects. This coupling between the thermal and viscous effects gives. 


* The distance in which the amplitude falls to 1/e of its initial value is given by twice the: 
reciprocal of (11.7). 
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rise to a loss additional to that calculated in § 6 on the basis of viscous losses only. 
The fundamental boundary condition is that V,, =0 at the surface, and the thermal 
losses result from the fact that the excess temperature 7” is related to V,. From 
(5.3) we obtain the approximate relationship 


OL LSe0V, 
Bt Coxe 
By (6.1) 
V,, =V,,°[1 — exp { —x(1 +2)(wpy/2)*3]. 
Hence 
Sa, - : 
T’ = 7 V1 ~i)(wpy/2n)![1 — exp {—x(1 +2)(wpy (2m) +++ (12.1) 


‘Thus the heat flow per unit time in the boundary layer is, per unit length, 
—K;,(@T/0x) 9 =P VOT Sp, Kylee ree or (12:2) 


By I, §6, the fraction of this heat lost to the walls is 2Z,/(Z,+2Z,). (12.1) and (12.2) 
give Z,=T’/HT =1/(1 +2)TK((p,/2u)?, and Z, is given by (11.5). Hence the 
transmission coefficient is approximately (2/K,)(uK,C,p,/p,)*, and the rate of heat 
loss per unit volume of liquid is BTV,°, where B =(2P/A)(S/C)(K,C,p,p,/p)*- 


‘The equation of entropy conservation is now 
O(pS)/ot+pSV.V,+BV,=0  ...... (1223) 
and the hydrodynamical equations become 
V(0,—<'p) +pV,—p'|u =, 
up’ — p,»V,—p.V.=9, 
u(p' S +pS’) —V,(pS +7Bu/w) =0, 
uV,+ST’ —p’'/p =0. 


Eliminating V, and V,, we obtain 


"= (Cpidp\l—<\-) ek eee (12.4) 
and 
df= (TSC \(pia 2) ee (12.5) 
where 
oP ase Pe PO ee ea , 
a = P A a (5 (2 = \+ wp A : G (K€ 5.7.) CO) Gieo so (12.6) 


For second sound, the distance in which the amplitude is reduced to 1/e of 
its initial value is thus ; 


A 4 2s (KsC.pspn\? | 
1/I(w/uUs) = puy (5) {p. (=) is a (See) i eet (12.7) 


whilst the corresponding distance for first sound is, to this degree of approximation, 
unaffected by the thermal conductivity. 

The attenuation due to this cause is very important for second sound propagation. 
With the values given in §§5 and 10, and taking K,~10~4 cal.cm*deg™? and 
C,~2 x 10° cal.gm™ for glass, (219/C,)(K.C.p,p,/4)! ~3 at 1-5°K.; this is to be 
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compared with p,(j.w/2p,,)! ~5 x 10-3f! at the same temperature. Thus at all but 
the highest frequencies the conduction effect associated with that due to viscosity 
is far more important than the direct effect of the viscous force. This result is 
hardly surprising when we recall that second sound is primarily a thermal effect. 


S133 DHE REELECTION OR {SOUND 


In this section we shall examine the problem of the reflection and transmission 
of a wave of sound incident normally of a plate of width d. Throughout, we 
shall use the concept of impedance introduced in §4. The results obtained may 
be applied to such problems as that of determining the disturbance of first and 
second sound propagation due to a thermometer or to the reflection of the waves at 
a metal plate forming one end of a resonance chamber. 

The boundary conditions at the two surfaces are given by the continuity of 
the ‘current’ J and the ‘potential’ Z7. Taking Z, to be the impedance of the 
plate and Z, that of the surrounding liquid, we obtain for the ratio of reflected to 
incident current 


irs (Zee Zeya) ok 2 eRe (13.1) 
I, (Ze2+Z2\1 —m2) 422,2(1 +m) 


and for the transmission coefficient 


I, 4mZ ZL. (13.2) 


of e) o> ees 


I, (Z2+Z2)(1 — m2) +2Z,Z,(1 +m?’ 
where m is the factor (in general complex) by which the amplitude of a single wave 
is reduced on passing across the plate. 


(1) Application to reflection of first sound. We take the velocity of motion as 
the ‘current’, and the fluctuation in pressure as the product of impedance and 
current. The impedance of the liquid is Z)=p/V =pott) where py and uy are 
respectively the density of the liquid and the velocity of first sound init. Similarly 
for the plate, Z,=p.u,. Since the current of first sound undergoes only a change 
of phase in passing through the plate, m =exp (—iw d/u,). 


(ii) Application to reflection of second sound. Here, as in §4, we take the flow of 
heat as the ‘current’, and the fluctuation in temperature divided by the mean 
temperature as the ‘potential’. By (11.1), the impedance of the liquid is 
Z,=1/pTCu., and by (11.5) that of the solid is Z,=1/87T(1+z2), where 
B=(4@p.:C oo Kooi)?» The heat current in the solid is attenuated; by (11.4), 


m= Cx d(1 +1)($HPso1Cso1 aa) | 
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APP EN Des 
THE JUSTIFICATION OF THE SOLUTIONS GIVEN IN §5 AND $6 


In §5, we obtained our solutions by ignoring some coupling terms in (5. 9) 
and (5.10). We shall now demonstrate that this procedure was justifiable. A 
similar justification may be given for the work of § 6. 

The values of u are given by the condition that (5 9) and (5.10) are consistent. 
We are at liberty to take any pair of the variables p’, p’, S’ and T” as independent. 


Taking 
Le dp / .) / 
= (sn), 7 ‘ (35 He 
; Op no (OS: 
= (5p)? + (55) 2% 
we obtain from (5.9) 


a «pas cpu? OS] 1, 
| was 1 Se les a) = o7 |t=0 pate (Al) 


and from (5.10), 
oS weSe@ oS uweS 0 
offen Lier ~ p (Gir e Pp lal 
EiGwieeesboees b> - 


These equations will be consistent if the determinant formed by the coefficients 
of p’ and T” is identically zero. Putting V =eu? =4ipw/3p, neglecting terms in V?, 
and then making some thermodynamical transformations, we obtain 


HTS? op a ESUe io 
4 _ 4,2 Ps a p Ps op 
fi : fe Cy us ee “s | z (Es Pn ar 


Vo TS? (ap\ (aT ap 
eee E + (2). ps(5-) - > (35) |-2 Sane ity (A3) 


The last two terms in (A3) and the difference between the adiabatic and isothermal 
compressibilities are proportional to the coefficient of thermal expansion, and 
may be neglected in determining the attenuation due to viscosity. Hence 


2u* = (uy°)? +(u99)? + Vp/pn + [(ur°)? — (2°)? +V (on —ps)/pa], +++ ++ (A4) 


where u,°=(0p/dp)? and u,°=(TS%p,/Cp,)?. The two solutions are those given 
by (5.13) and (5.14). 


and 
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ABSTRACT. The angular distribution of radiation from a synchrotron has been measured 
with an ionization chamber for a number of different target materials and thicknesses, and 
the results have been compared with a simple theory in which the radiation from a parallel 
beam of electrons striking a thin plate has been calculated. Detailed comparison with theory 
is not possible, because the precise way in which the electrons hit the target cannot be 
determined, nevertheless qualitative explanations have been found for most of the effects 
observed. ‘The most important of these is the relative insensitivity of forward intensity and 
beam width to target size and material. 


ol UNER OD UC DilOiN 


N this paper the results of a preliminary study of the effect of target size and 
material on the radiation from a synchrotron are described. Several 
interesting effects have been observed, though further work would be needed 

to obtain a full understanding of the problem. Unfortunately the way in which 
the electrons hit the target is very complicated, so that a detailed comparison 
with theory is not possible. Nevertheless a simple idealized theory has proved. 
to be a useful guide. 

Polar diagram measurements taken with absorbers sensitive only to the higher 
energy quanta have been reported by other workers (for example, Allen- Williams. 
and Appleyard 1949 a,b, Lees and Metcalfe 1950); their results are in fair 
agreement with those given here. 


§2. MEASURING TECHNIQUE 


The machine used in these experiments is described by Fry e¢ al. (1948). 
It was normally run at an energy of 10 Mev., with an internal target, and electron 
spiral pitch at the target of 16x 10-*cm. ‘The radiation intensity was measured. 
with a cubical ionization chamber of side 9 cm., the front wall of which could be- 
varied in thickness. The chamber could be moved about on a table in front of 
the machine, and the polar diagram measurements were taken with it 1 metre 
from the target. The chamber subtends an angle of 5° at the target, so that for 
narrow diagrams a check of the shape near the maximum was made at a distance 
of 24 metres. 

The chamber simply measures ionization produced, it gives no indication at 
all of the quality of the radiation falling on its front wall. Some indication of 
this can be obtained by varying the wall thickness and examining the absorption 
curve produced, though such information is not usually precise enough to be of 
much value. It is easy to distinguish the hard gamma rays (1 to 20 Mev. say), 
from the softer radiation, but it is difficult to be much more precise than this. 
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§3. THEORETICAL DISCUSSION OF TARGET RADIATION 


The electrons in a synchrotron give up their energy to the target partly by 
ionizing the atoms of which it is composed, and partly by making radiative 
collisions (bremsstrahlung). In addition they also undergo considerable elastic 
nuclear scattering. It is this elastic scattering which plays the predominant 
part in determining the polar diagram of the radiation. 

It is instructive to consider what happens when a parallel beam of electrons 
falls on a thin plate of uniform thickness, although this is far simpler than the 
conditions which exist in a synchrotron, where the target geometry is much more 
-complicated. 

If the target thickness is expressed in ‘radiation lengths’ (defined as 


A|(4aNZ?r,? In 1832-1), 


where 7)=Cclassical electron radius =2-82 x 10-"%cm., N=Avogadro’s number 
= 6-02 x 1023 and «=fine structure constant =1/137) the angular distribution 
of elastically scattered electrons (for small angles of scatter) after passing through 
.a thickness t¢ is, as shown in the Appendix, 


2 


NEY kerry 


exp (707 44.02), ee eeenveeees Cl 
where E is the electron energy (assumed large compared with the rest energy) 
expressed in Mev. ‘This expression is obtained by taking a multiple convolution 
-of the formula for single scattering, and is therefore inaccurate when ft is very small 
-or 0 very large. 

Now if the target is thin compared with one radiation length, the energy loss 
in it will be small, and particles will be equally likely to have radiative collisions 
at all values of ¢. If then F(t, @) is integrated over ¢, and | F(t, #) dt is combined 
with the radiation polar diagram for a single collision, the result will be the overall 
polar diagram of the radiation from the target. Assuming intuitively that the 
radiation polar diagram for a single collision is a gaussian curve of r.m.s. angle 
\/zu/E (where w is the rest energy of the electron) the resulting polar diagram, 
normalized to be unity at 0=0, may be written 


P(E6) =[ — Eif —(E6)?/440t +2u2)} + Eif — (E0)?/2u2}]/In (1 +220¢/y2). 


From this we see that the polar diagram width is inversely proportional to the 
incident electron energy, and is the same for all materials of equal thicknesses 
expressed in radiation lengths.* Polar diagrams for different values of ¢ are 
shown in Figure 1. : 

The variation of forward intensity with target thickness may be readily found. 
It is convenient in this connection to define a quantity used in aerial theory called 
the ‘gain’ of the angular distribution; this is given by 


G=4n | | P(E0)2n sin 0 db 


*This treatment is due to Schiff (1946), though in the analysis given “here a simpler form is 
taken for the polar diagram of a single radiative collision. This simplifies the analysis and is 
justified by the fact that no accurate information is readily available. Recent calculations by 
‘Stearns (1949) show that the r.m:s. angle of radiation from a single collision is a slowly varying 
function of E/, the photon energy and the nuclear charge, and also that the angular distribution 
as not gaussian. The theory presented here must therefore be treated with reserve. 
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and would be equal to unity for an isotropic distribution. Now when 6 is small 
so that sin@~6, it may be shown (see Appendix) that 


G late 220 rat\ Se eye kes, (3) 


Kk 
~ 1100 
Now if ot is the proportion of the electron energy converted to radiation in a 
thickness ¢ of the target, the fractional forward radiation intensity per unit solid 
angle is given by 
pCO Ai Mie ov 4) eee (4) 
‘The quantity Gt/E? is plotted against target thickness t in Figure 2. The curve 


Intensity 


-200 -100 0 100 200 


MeV. x Degrees 


Figure 1. Theoretical polar diagrams for various thicknesses of target; horizontal scale 
in Mey. X degrees and vertical scale in arbitrary units of intensity. 
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Figure 2. Theoretical curve of variation of forward intensity of radiation with target thickness. 


Note: the horizontal and vertical scales should be multiplied by 2. 


rises steeply when ¢ is small, but at larger values of ¢ where the polar diagram 
width is controlled mainly by the scattering process, the forward intensity is 
relatively little affected by the target thickness. 
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The energy distribution of the radiation from a thin target will clearly be that 
from a single collision. For a thick target this distribution will be degraded, 
because some electrons will lose appreciable energy by ionization before radiating, 
and others will make more than one collision. The effect of the degradation from 
both these causes will be most marked in the radiation at large angles to the 
forward direction. 


§4. EXPERIMENTAL RESULTS 
(i) Preliminary Experiments on the Nature of the Radiation 


Measurements were first made at 10 Mev. of the variation of ionization chamber 
reading with wall thickness. These showed that under some conditions there 
was a large amount of very soft radiation, which could be entirely stopped by about 
1:5 gm/cm2 of matter, superimposed on the expected hard radiation. Further 
investigation showed that this soft radiation had a broader polar diagram than the 
hard component and appeared to come from about 10cm. outside the target. 
It was concluded therefore that the soft radiation consisted of electrons which 
had passed through the target, and been deflected sufficiently by it to escape from 
the magnetic field altogether, finally passing through the glass donut wall. The 
path of these electrons would, of course, be curved while they were under the 
influence of the magnetic field, causing them to appear to come from a point 
outside the target. 

Increasing the ionization chamber wall about 1-5 gm/cm? did not affect the 
shape of the polar diagram, so that for all subsequent experiments this thickness. 
was used. 


(11) Polar Diagram Measurements for Different Targets 


In Figure 3 a number of polar diagrams with different targets all measured at 
an electron energy of 10 Mev. are shown. It will be seen that these are somewhat 
similar in shape to the theoretical diagrams of Figure1. The apparent thicknesses. 
of the targets deduced from the widths to half intensity are of the right order of 
magnitude, though the radiation at large angles is greater than might be 
expected. 

The electron orbit contracts towards the target with a spiral pitch of about 
1-6 x 10-°cm.*, so that if the target is very thick, only its outer edge will be 
struck. The polar diagram will therefore be narrower than it would be if the 
whole target were being used. Allen-Williams and Appleyard (1949 a, b) have 
pointed out that this effect would account for the unsymmetrical shape of the 
diagram which is noticeable with the thick targets; electrons deflected inwards 
pass through more of the target than those deflected outwards, and thus more 
energy is radiated on one side than the other. 

A slight broadening of the diagrams from thin targets might be expected, 
since as shown in § 4(v), it is possible for electrons to pass more than once through 
the target. Such electrons, when they arrive atthe target for the second time, 
would in general be moving at an angle to their original direction, so that effectively 


the beam striking the target would be divergent. This would broaden the polar 
diagram. 


* The effective spiral pitch will however be much greater than this, since the electrons will be 
executing radial oscillations. 
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‘These considerations for thick and thin targets explain the fact that there is 
less variation in beam width from the various targets than might be expected 
from a direct application of the theory. ‘The excessive radiation at large angles 
might be partly due to multiple radiation collisions in the target, but is probably 
also due to electrons which, after passing through the target, make radiative 
collisions in the doughnut walls. 


(i) Effect of Changing the Electron Energy 


Diagrams from a 0-03 cm. diameter cylindrical tungsten target were measured 
at energies of 10 and 5 Mey., and the beam widths to half maximum amplitude 
were found to be 18° and 36°. @ thus remained constant as demanded by the 
theory. 


(iv) Relative Forward Intensities for Different Target Materials 


The relative forward intensities from targets of different materials were 
measured with the aid of movable targets. The results for cylindrical targets of 
0-085 cm. diameter made of tungsten, copper and aluminium, measured over 
the energy range 3-20 Mey. are shown in Figure 4. Although the thicknesses of 


Degrees 


Figure 3. Polar diagrams of various targets measured at an electron energy of 10 Mev. 


Curve A B € D E 
ateri Gold Aluminium Tungsten Tungsten ‘Tungsten 
ae Foil Wire Wire Wire Wire 
Thickness (cm.) 0-0013 0-085 0-013 0-05 0-2 
Thickness (radiation lengths) 00042 0-009 0-039 0-15 | 0-6 
Beam width to half maximum iBed. qe 18 20°5 24 
Apparent thickness (radiation lengths) 0-02 0-008 0-05 0-18 — 


the targets measured in radiation lengths are very different, the forward intensities 
from them are comparable. This is what would be expected from Figure Ux, 
though the value for aluminium is too high. This is probably due to the multiple 
passage effect described in the next section. No accurate comparative measure- 
ments were made with tungsten targets of different diameters, though it was found 
that the forward intensity from a target of diameter 0-085 cm. was only about 15%, 
higher than that from a target of diameter 0-013 cm. 
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(v) Experiments to Determine whether Multiple Passage through the Target Occurs 


It is suggested in § 4(ii) that electrons may pass through the target more than 
once. An electron might get scattered through a small angle with little loss of 
energy, and be refocused on the target (possibly after several revolutions) by the 
magnetic guide field. The scattered electrons leaving the target travelling at 
too wide angles or with too little energy would hit the donut walls or resonator, 
but those within about 3° which had not lost much energy would eventually pass: 
through the target again. Since the period of oscillation bears no relation to 
the period of rotation (it will in fact be a function of the oscillation amplitude 
since the field index m varies with radius), the target might not be hit again for 
several revolutions, depending on its size. If now a larger target is placed at a 
slightly smaller radius than the original one, there will be a greater chance of 


Intensity 


=30) 0 30 
Degrees 


Intensity Ratio 


Figure 5. Curves obtained from the ‘ multiple 
passage ’ experiment using an aluminium. 
target. 

o i shee ° - A. Orbit spiral 2107? cm., outer 
; target hit. 


Figure 4. Ratio of forward intensity from 
0-085 cm. diameter targets of aluminium and 
copper to that from an equal tungsten target. 
The target thicknesses in radiation lengths 
are: tungsten, 0:27: “copper, 0°06: 
aluminium, 0-009. 


B. Orbit spiral 1:6 10-5 cm., both 

targets hit. 

The relative intensities of the two curves 
are uncertain, but they are plotted for 
convenience with the same maximum. 
value. 


electrons hitting the large target on subsequent revolutions, unless the scattering 
angle is very small and the oscillation amplitude is less than the radial distance 
between the targets. 

Accordingly two targets, each consisting of a vertical wire, were placed in the 
synchrotron at a very small radial separation, but an azimuthal separation of 15°; 
it was thus possible to see which target was being struck by measuring the output 
polar diagrams. An aluminium target showed considerable multiple passage when 
the electron spiral pitch was 1-6 x 10->cm., but very little when the spiral pitch 
was increased to 2 x 10-3cm.* with the aid of pulsed contractor coils above and 
below the orbit. The diagrams obtained in this experiment are shown in Figure 5 
and some results are tabulated below. 'Thicknesses are in radiation lengths. 


Outer target Thickness Inner target Spiral pitch Results 
0-0125 cm. W 0-038 0:05 cm. W 1:6 10-%cm. No multiple passage 
0-085 cm. Al 0-009 0-05 cm. W 1:6 10-5cm. Multiple passage occurs 
0-085 cm. Al 0-009 0-05 cm. W 2x10 cm. Very little multiple passage 


* Changing the spiral pitch by this amount appeared to make no difference to the output or 
polar diagram of the thicker targets. 
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‘The radiation from the aluminium target is in a narrower polar diagram when there: 
is little multiple passage, as might be expected. 

A second similar experiment with a gold foil target of thickness 0-0013 cm. 
(0-0043 radiation length) likewise showed multiple passage effects, both at spiral 
pitches of 1-6 x 10->cm. and 2 x 10-3cm. 


95. CONCLUSION 

Experiments have shown that it is not possible to exert much control over the 
intensity or angular distribution of synchrotron target radiation by varying the 
target material or geometry. Measurements on a series of targets show a vari- 
ation by a factor of about two in both beam width and intensity. These results 
are in good general agreement with a simple theory. The theory cannot however 
be applied directly, because multiple passage occurs in thin targets, making them 
effectively thicker, whereas only the outer edge is effective in thick targets, which 
therefore appear thinner. 
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APPENDIX 


DERIVATION OF EQUATIONS (1) AND (2) 

Rossi and Greisen (1941, §23) show that the distribution of electrons after 

multiple scattering through a thin plate is approximately gaussian. For high 
energy electrons the distribution function may be written 


F(t, 8) = exp(—E?67/8802), eae ae (1) 


Fee 
880rt 
where E is the energy of the electrons in Mev., and ¢ is the thickness of the plate 
in radiation lengths. ‘The function has been normalized so that its integral over 
all angles is unity. 

Now if the target is so thin that the energy loss of electrons passing through 
it is small, and the scattering angle is not too large, there will be an equal amount 
of radiation from a thickness 5¢ of the target at all values of t. The angular 
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distribution of the radiation from a layer of thickness 6¢ at a depth ¢ in the target 
will be a convolution of the expression in equation (1) and the expression for the 
radiation polar diagram for a single radiative collision. (The convolution of two 
two-dimensional gaussion distributions of r.m.s. angles 4/20, and 1/26, 
(written as exp [ — 02/20,2]* exp [ — 6?/20,”]) may be shown to be 
677054 

~when 6, and 0, are small.) The polar diagram of a radiative collision is assumed 
‘to be of the form 


2 exp [—6°/2(0, + 42°)] 


2 
at Spee 
| (6) One exp (= S072)... aed ees secre (5) 
where jz is the rest energy of the electron. 
If now odt is the proportion of the electron energy converted to radiation in 
-a thickness 5¢ of the target, the radiation per unit solid angle from such a layer will 
be given by 


2 E2 
a SR 2p2 * _ 17 292),,2 
dR=odt 440at exp ( — E?6?/440t) Dat exp (— $£76?/17) 
Eodt 
ae aN ste ee 16272 2 
Pa 220 + 2) exp [ — $07E?/(220¢ + p?)], 
whence by proceeding to the limit and integrating from 0 to ¢ 
ok? : : : 
IR 4400 [ — Ei{ — $276?/(220¢ + ?)} + Ei{ — $270? /.?}]. 


Equation (2) may be obtained by normalizing this expression to be unity when 
6 =0. 

The gain may be found by putting 6=0 in equation (6); this gives 
ok? 
4407 


Now the fraction of energy radiated per electron is ot, so that Ry_)= Got/4z, 
whence 


Rico= In (1 +2202/p2). 


2 


E 
os 110¢ 


In(1 +220¢/u2). 
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The function — Ei( —y) is defined as | y-te-" dy and 1s tabulated in Jahnke 


y 
and Emde, Tables of Functions, p.6. When @=0 equation (2) is indeterminate 
and the asymptotic relation — Ei(—y,) + Ei(—y.) =In(y,/y,) must be used. 
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CETTE RS= tO THE EDITOR 


Measurement of Polar Diagram of Synchrotron Gamma 
Radiation 


In support of work by Lawson (1950), described in this issue, it was considered 
desirable to measure the polar diagram of the radiation from a machine similar to his by 
using a detector sensitive to high energy quanta only. For this purpose copper was irradiated 
to give the **Cu positron activity with half-life 10-5 minutes as suggested by Allen-Williams 
and Appleyard (1949 a, b). Baldwin and Klaiber (1948) have measured the cross section for 
the activation, and copper may be considered to be effectively a resonance detector at 22 Mev. 
A further purpose of the experiment was to use detectors smaller than those used by Lawson 
or Allen-Williams and Appleyard and to measure the diagram in two planes at right angles. 

The machine used in the experiment was described by Fry et al. (1948). The electrons 
were allowed to shrink on to a tungsten target by turning off the radio frequency just before 
the peak of the magnetic field cycle, and struck the target with an energy of 24 Mev. The 
‘copper samples were about 1 cm. square and 1 mm. thick, weighing 1 gm., and the activity 
after irradiation was measured with a standard G.M.4 end window counter. Detectors in 
various positions were given a standard irradiation as measured by an ionization chamber. 
The results are plotted as initial induced activity against position of sample in the Figure, 
‘curve (a) showing the horizontal polar diagram and (6) the vertical diagram, the vertical 
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Polar diagrams of radiation from 24 Mev. synchrotron in region 11-24 Mev. 
(a) Horizontal polar diagram. (b) Vertical polar diagram. Obtained with tungsten target 
consisting of a 0°005 in. foil bent into a cylinder of 0-36 cm. diameter. Orbital 
spiral pitch of the order of 10~° cm. 


diagram being measured at the peak of the horizontal one and vice versa. ‘Two sets of 
results were taken at different distances from the x-ray source, the ordinates being corrected 
by assuming an inverse square law. ‘The agreement may be taken as a measure of the 
accuracy of the experiment. It will be seen that the horizontal curve is not symmetrical ; 
it shows the broadening due to deflection of the electrons towards the thicker parts of the 
target, as mentioned by Allen-Williams and Appleyard, and by Lawson. The vertical polar 
diagram is, however, quite symmetrical. 


ic E Research Establishment, ID), Ih, Wisisey. 
oo pais Harwell. L. H. METCALFE. 


17th November 1949. 


ALLEN-WILLIAMS, D. J., and APPLEYARD, R. K., 1949 a, Brit. . Rad., 22, 106; 1949 b, Proc. Camb. 


Phil. Soc., 45 (Part 2), 305. 
7 _C., and Kratser, G. S., 1948, Phys. Rev., 73, 1156. 
Cee Ciuzor. J. W., Goward, F. K., and Darn, fe 1948, Nature, Lond., 161, 504. 


Lawson, J. D., 1950, Proc. Phys. Soc. A, 63, 653. 
PROC. PHYS. SOC- LXIII, 6—A A4 
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Identification of Photo-Disintegration Stars in 
Nuclear Emulsions 


Nuclear emulsions, irradiated by a continuous spectrum of gamma-rays extending to 
energies between 23 and 25 mev., have been found to contain charged particle stars arising 
from the photo-disintegration of boron, carbon and oxygen (Goward, Titterton and Wilkins. 
1949, 1950, Goward, Telegdi and Wilkins 1950). In this work the criterion of “momentum 
balance’ was used to develop accuracy of measurement, and the present letter states the 
limits of resolution and identification so far achieved. This information is essential to a 
discussion of a new reaction in nitrogen, given in a companion letter henceforward referred 
to as B (Wilkins and Goward 1950). It may also be of assistance in similar nuclear emulsion 
work where the measurement accuracy is not explicitly determined. 

If the identities of the particles of a photo-disintegration star are known, the particle 
energies F'.(r=1, 2...) in Mev., and the momenta, p,, in units of (Mev x mass number)", 
may be determined. The momentum, p,, of the gamma-ray is (Ey+,)/43-1, where 
Ey=E, and E, is the threshold of the reaction, so that, for conservation of momentum 
A=(2p,)—p,=0. Thus | A| is a convenient measure of the overall accuracy. 

When the reliability of an experimenter has been established by statistics on attested 
reactions, his ability to distinguish other reactions may be inferred. Such a test is necessary 
since wide variations have been found in the statistics of individuals; these variations should 
disappear if accuracy statistics are used to cultivate a measurement technique. Some 
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Figure 1. Statistics of | A] for boron loaded EI plates. 
(a) 228 carbon star measurements by the authors. Superposed is a curve of the form 


N=(K| Al)? exp —(K| A))?. 
(6) 51 oxygen star measurements. 
(c) Measurements on 31 stars attributable to °B(y, D) or “B(y, T) reactions. 
(d) Carbon star measurements by three other experimenters, considered representative. 


statistics for measurements by the authors on carbon and oxygen stars in boron-loaded 
Ilford E1 emulsions are shown in Figures 1 (a) and 1 (8). |A| was found to be almost 
independent of Ey or E,, unless a particle had a rangé less than about 2 microns. An 
approximate theoretical curve is fitted to Figure 1 (a) by assuming that the measurements 
give three rectangular components of A independently and with equal accuracy. No 
re-measurements were taken after | A] had been evaluated, so that the statistics give a realistic 
estimate of the accuracy expected in unknown reactions. Figure 1 (d) shows statistics 
which are considered typical of uncultivated measurements by otherwise experienced 
individuals. All the authors’ measurements were made, using a magnification of 1400, on 
unselected stars; emulsion shrinkage factors were measured directly by a micrometer. 

The allocation of stars to appropriate energetically possible reactions is conveniently 
accomplished by first considering that all stars consist of alpha-particle tracks only and 
calculating corresponding | A| values. The results, for all three-particle stars in ie Ei 
emulsions of Figure 1 are shown in Figure 2. Almost complete resolution of the carbon 
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stars is apparent, with |A|<1-0. For all stars with | A] >0-8, appropriate tracks are recon- 
sidered as ?H, *H, *Li, etc., and the possibility of overlooking a fourth track is also investi- 
gated (see B); the appearance of the tracks serves to confirm the | A] analysis. Figure 1 (c): 
shows, for example, how a group tentatively attributable to the reactions UB(y, T) and 
B(y, D) is resolved in this way. The “B(y, pn) reaction gives an exceptional group of 
stars in which |A | +0, owing to the unobserved neutron present. Stars attributed to one or 
other of the boron reactions are indicated by full shading in Figure 2, but the extent to which 
they may be distinguished from each other, in relation to the results of the present letter, 
will be described in a later publication. 


Number of Stars 


Figure 2. |A| values for 3-particle stars, treating all tracks as alpha particles. 


Unshaded rectangles indicate 228 carbon stars, fully shaded rectangles indicate 47 boron 
stars, and half-shaded rectangles indicate the 8 events discussed in B. 


Figures 1 (6) and 2 represent all three- and four-particle stars observed, including those 
discussed in B, but excepting several stars readily identified (from appearance, energy and 
momentum balance considerations) as radioactive contamination stars or collisions (between, 
for example, protons and protons). Only 3% of the events, shown in Figure 1 (a) as the 
tail of the carbon star statistics lying outside the fitted curve, can be considered as of doubtful 
identity. 


Atomic Energy Research Establishment, F. K. Gowarp. 
Harwell, Didcot, Berks. Je ine WitING 
3rd February 1950. 


Gowarb, F. K., Tetecpi, V. L., and WILKINS, J. J., 1950, Proc. Phys. Soc. A, 63, 402. 

Gowarp, F. K., TirrEertTon, E. W., and WILkins, J. J., 1949, Proc. Phys. Soc. A, 62, 460; 1950, 
Ibid., 635 172. 

Wikxins, J. J., and Gowarp, F. K., 1950, Proc. Phys. Soc. A, 63, 663. 


Identification of Nitrogen Photo-Disintegration Stars 
in Nuclear Emulsions 


Of the 283 three-particle stars found in the emulsions discussed in the companion letter 
(Goward and Wilkins 1950), six have been attributed to the reaction 


UN hy — *Li1--2*He—1lo-l wey, = ane (1) 


The following discussion of the certainty of this assignation provides an excellent example 
of the value and limitations of the identification methods described. eS 
For the present purposes, a sufficiently accurate range—energy relation for “Li nuclei in 
nuclear emulsions was derived directly from that for alpha-particles (by adapting to ®Li 
nuclei the electron capture and loss data for alpha-particles reviewed by Knipp and Teller 
(1941)). The momentum of a °Li nucleus of a few Mev. energy is nearly twice that of an 
alpha-particle of the same range. An important consequence follows, namely that the main 
difficulty in identifying stars from reaction (1) is to avoid confusion with the four-particle 


reaction 
16) --Ay —> 44#He—14-5mev., § jj‘  ....6. (2) 
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rather than (as might be expected) with the three-particle reaction 
2C+ hy + 34He—7°3 Mev. eras tees (3) 


As shown in Figure 1, the latter differentiation is quite conclusive with the present resolution. 

A reaction (2) star, however, may be mistaken for a three-particle star when the fourth 
track is (a) nearly coincident with another track or (b) too short to be detected. In either 
case, if the vector difference in momentum between the undetected particle and any other is 
small, satisfactory momentum balance when interpreted as a reaction (1) event is probable. 
Small vector differences in momentum have, in fact, proved to be common, since a ground 
state ®Be nucleus is a frequent intermediate product in the formation of oxygen stars. 
If E* is the energy released when an *Be nucleus disintegrates into two alpha-patticles of 
energies H, and Ez, then 


F*=4( EE Es) = (EE) cos by ee eee (4) 


where ¢ is the angle of separation of the alpha-particle tracks. Here &* is only of the order 
of 0-1 Mev., and the resulting tracks will therefore be little separated. Two coincident 
tracks (6=0) of slightly different lengths may be obtained, and have in fact been detected on 
close examination ot appropriate stars. (Further details of the oxygen disintegrations must 
be left for a future publication.) 
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disintegrations considered as of stars with energy of 
(a) 4N disintegrations. incident quanta. 
(b) 12C disintegrations. (a) eight suspected ‘4N 
(c) 18O disintegrations. disintegrations. 


(6) 51 1°O disintegrations. 


With the above possibilities in mind, eight stars remained with no detectable fourth 
track, good momentum balance as nitrogen events and poor balance as carbon events, as 
shown in Figures 1 (a) and 1 (6). The energy and momentum of a possible undetected 
fourth track may be calculated from equation (4), putting 60 and E*=0-1 Mev. The 
momentum unbalance data obtained on this assumption is shown in Figure 1 (c), and 
comparison with Figure 1 (a) gives no conclusive result. “Moreover, the identity of two of 
the stars, indicated by shading in the Figures, at once becomes doubtful when it is observed 
that their possible fourth tracks would be practically undetectable (less than one micron 
long). The remaining six stars, if from reaction (2), should show detectable fourth tracks 
1:5 to 3 microns long (alpha-particle energies of 1:4 Mev. (4:7 microns) and 0:6 Mev. 
(2-4 microns) correspond to a ®Li energy of 2:6 Mev.). Coincident tracks of about this 
length have proved to be distinguishable, as mentioned above, and on this evidence the six 
events may be assigned to reaction (1). 

A final test is to plot the total energy of the particles of each star, Ey, plus 16-1 Mev. 
i.e. the incident gamma-quantum energy £,, as in Figure 2 (a). The same two events fe 
above (which have low EF, values) are again suspect, but the identity of the remaining six is 
confirmed. ‘This is best illustrated by comparison with Figure 2 (6), which shows E,, (or 
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Ey + 14°5 Mev.) for a greater number of oxygen stars. Consistency with the known peak 
energy of the gamma-ray spectrum, and with the relative threshold energies of the two 
reactions is satisfactory. On the basis of the six selected stars, reaction (1) may be assumed 
to occur with a cross section of the oider of 10-28 cm? at around 25 Mev. 

A possible alternative reaction in nitrogen is 


MN +hy > #H+34He—-17-6 Mev. you 3 (Gy) 


There is no definite evidence for this reaction, although it would be possible to consider 


some of the shorter alpha-particle tracks in oxygen stars as deuterons without significantly 
altering the momentum balance. 


Atomic Energy Research Establishment, J. J. Witkrns. 
Harwell, Didcot, Berks. F. K. Gowarp. 
6th February 1950. 


Gowarp, F. K., and Wirxns, J. J., 1950, Proc. Phys. Soc. A, 63, 662. 
Knipp, J., and Tetrer, E., 1941, Phys. Rev., 59, 659. 


Alpha-Particles from 7’Al+D 


Groups of alpha-particles from the reaction ?’Al(d, «)?®>Mg have been observed, using 
a differential ionization chamber with electron collection in argon. A thin target of 
aluminium, about 0-5 mm. air equivalent evaporated on brass, was bombarded with 
deuterons of 0-93 Mev. and the alpha-particles emitted were observed at 90° to the incident 
beam within a cone of semi-angle 10°. An air absorption cell of fixed length and variable 
pressure was mounted between target and chamber, so that particle ranges from 2:0 to 
6:0 cm. could be measured. 


The variation of counting rate with absorber pressure is shown in the Figure. Seven 


Counts perza.of Beam (arbitrary units) 


0 10 20 30 40 50 60 70 80 
Absorption Cell Pressure (cm.Hg) 


Alpha-particle groups emitted in the bombardment of aluminium by deuterons 
of 0:93 Mev. energy and from a ThC source. 


groups of particles can be distinguished. ‘Those lettered B, C, D, E and F can be attributed 
to the 27Al(d, x) reaction. Group A is due, mainly at least, to the *C(d, «) process in a 
thin layer of carbon contamination. The detector was biased so as not to pecan pulses 
from single protons or deuterons, but group G can be ascribed to protons from the O(d, p) 
reaction, sufficiently intense to give multiple pulses. Similarly, the rise of counting rate 
at pressures below group A is caused by scattered deuterons. (By changing the bias on 
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the detector, it was always possible to distinguish alpha-particles from protons or deuterons.) 
Groups A and G, and also the rise of counting rate below A, were observed when a blank 
target was bombarded under the same conditions. 

In order to convert the pressure measurements to mean ranges, the alpha-particle groups 
from ThC and from “C(d, x) were used. Group H (see Figure) was obtained by mounting 
on the target a thin source of thorium-active deposit; the position of the carbon group was 
established accurately during the blank run. The ranges of these groups were taken as 
4:70 and 2-91 cm. respectively, the latter being based on recent measurements by Buechner 
et al. (1949), and the Cornell 1938 range-energy relation (Holloway and Livingston 1938). 

The ranges, alpha-energies and QO values for the Al(d, «) groups are tabulated below. 


Group Absorber pressure Mean range  a-energy 


(cm.) (cm.) (Mev.) (Mev.) 
B 23°8 3:34 4-83 4°75 
€ 29°8 3-66 5°14 5-08 
D 41-0 4-24 5-64 5:68 
E 47-7 4-60 595 6-04 
F 58:8 5:18 6:43 6°62 


An error of about —-0-:05 Mev. is associated with any of the energies or Q values given in 
the Table. 

It seems certain that group F corresponds to the transition to the ground state of 
25Mg; the QO value of 6°62 Mev. may be compared with the figure 6°46 Mev. given by 
Pollard et al. (1949), who observed four groups in all at 3 Mev. deuteron energy; in the 
present experiment the lowest of these would be at the position of the carbon group. 

The QO values given in the Table indicate excited states in *°Mg at 0:58, 0:94, 1:54 and 
1-87 Mev. above the ground state. The first two of these agree well with the results of 
Allan e¢ al. (1949). The energy release in the ground-state transition can be taken with the 
Q value 5-03-+0-05 mev. for the ?4Mg(d, p) reaction (Allan and Wilkinson 1948), to give a 
O of 1:-59+0-07 mev. for the reaction ?’Al(p, ~). A recent study of this reaction (Freeman, 
private communication) has yielded a value of 1:58-+-0-01 Mev., with which our result is thus 

_-consistent. 

One of us (P.B.T.) is indebted to the Royal Commissioners for the Exhibition of 1851 

for the award of an overseas scholarship. 


Cavendish Laboratory, A. P. FRENCH. 
Cambridge. Pee B weuRbAGY. 


24th February 1950. 


ALLAN, H. R., and WILKINSON, C. A., 1948, Proc. Roy. Soc. A, 194, 131. 

ALLAN, H.R., Winxkinson, C. A., Burcuam, W. E., and Curtinc, C. D., 1949, Nature, Lond., 163, 
210. 

BUECHNER, W. W., Strait, E. N., SPERDUTO, S., and Mam, R., 1949, Phys. Rev., 76, 1543. 

Ho.ioway, M. G., and Livinecston, M. S., 1938, Phys. Rev., 54, 18. 

Potrarb, E. C., Sattor, V. L., and Wyry, L. D., 1949, Phys. Rev., 75, 725. 


Alpha-Particles from °F Bombarded with Deuterons 


It is well known that deuteron bombardment of !8F causes the emission of several 
alpha-particle groups (Burcham and Smith 1938). These have been ascribed to the 
reaction “F(d, «)!"O, but the process F(d,n), followed by break-up of the resulting 
*°Ne* into 16O and an alpha-particle, is also energetically possible. We have made a study 
of the alpha-particles to find which, if any, of the groups were due to this second process, 
and whether they could be related to known levels in 2¥Ne (Bonner 1940, Powell 1943). 

The method was simply to observe the dependence of particle energy on deuteron 
bombarding energy. In the process 19F(d, n) nearly all of the bombarding energy is carried 
aviay by the neutron so that the variation of alpha-particle energy is small. In the reaction 

F/d, x) on the other hand, most of the kinetic energy of the colliding nuclei is given to the 
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alpha-particle. In our experiment alpha-particles from a thin CaF, target (0:5 mm. air 
equivalent) were observed at 90° to the incident beam, and bombarding energies of 0-63 
and 0:93 Mev. respectively were used. The expected shifts of mean alpha-particle energy 
under these conditions are +-0-22 mev. for the direct transition and —0-005 mey. for the 
two-stage process. For alpha-particle energies from ~3-5 Mev. to ~7 MeV., measurements 
were made with a differential ionization chamber and absorption cell; the region between 
2 and 3 Mev. was studied with a magnetic analyser and scintillation counter (Burcham and 
Freeman 1949). 


The results are shown in the Figure in the form of numbers—energy curves. No attempt 
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Energy distributions of alpha-particles emitted by fluorine bombarded with deuterons 
of 0:63 and 0:93 Mey. 


has been made to link the magnetic analyser measurements with the rest. In addition 
measurements on the lowest peak of the curve obtained using the differential chamber were 
repeated, using a higher bias setting of the discriminator (full circles) to ensure that the 
position of the maximum was not affected by scattered deuterons. Seven peaks can be 
distinguished. ‘That labelled A is due to an intense group of protons and was therefore 
not studied in detail. Group D arises from carbon contamination on the target, as was 
verified by its position and the increase of its intensity with time. ‘The remainder can be 
attributed to fluorine. It will be seen that the positions of peaks B, E, F and G change 
noticeably with bombarding energy, whereas peak C is broadened (as one would expect), but 
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only slightly shitted; From the diagram the following table can be constructed; £}, indicates 
bombarding energy. 


Group a-energy (MeV.) Difference 
E\y=0-63 Mev. E»=0-93 Mev. (Mev.) 
B 2-70 2-93 0:23 
G 3-90 3-96 0-06 
E 4-71 4-91 0-20 
F Sesh 5:54 0:23 
G bes 6:16 0-23 


One can conclude that groups B, E, F and G are due to the reaction F(d, «)?’O, the 
corresponding Q values being 2°81, 5:26, 6:03, 6:79 Mev. The last three of these are 
consistent with values obtained by Burcham and Smith (1938), who did not study lower 
energies. Assuming that group C arises from the two-stage process described above, 
one obtains a value of 9-7 Mev. for the corresponding level in ?°Ne. This is not in good 
agreement with results obtained on neutron spectra by Bonner (1940) and Powell (1943), 
but it is difficult to find another mechanism which could be responsible for this group. 
The shift of the peak of group C, though slight, appears to be quite genuine and may point 
to the presence of a relatively weak unresolved 1°F(d, «)-group in this region. No group 
was found in the present experiment corresponding to the 7:8 Mey. level in ?°Ne reported by 
Powell, nor has this level been identified in experiments on scattering of alpha-particles by 
oxygen (Brubaker 1938, Devons 1939, Ferguson and Walker 1940), so it seems likely that 
alpha-emission from this level is forbidden. 

Thanks are due to Miss J. M. Freeman for making available the magnetic analyser used 
in these experiments. This work has been made possible for two of us (P.M. and P.B.T.) 
by assistance from the British Council and the Royal Commissioners for the Exhibition of 
1851 respectively. 


Cavendish Laboratory, A. P. FRENCH. 
Cambridge. P. MEYER. 
24th February 1950. P. B. Treacy. 
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The Energy-Release in some (pa) Reactions in Light Nuclei 


Q-values for the reactions ?7Al(px)?4Mg and ?°Na(pa)2°Ne have been determined, as 
already described (Freeman and Baxter 1948), by measuring the differences between the 
ranges of the alpha-particle groups obtained in these reactions and that from the reaction 
“F(px)O*; the range of alpha-particles from the latter was found absolutely by comparison 
with a thorium-C source. The corresponding energies were taken from the 1938 Cornell 
range-energy curve for alpha-particles (Holloway and Livingston 1938), and the Q-values 
then computed. Recent determinations, based on direct energy measurements, of the 
energy-release in the reaction !*F(pa)!*O* (Burcham and Freeman 1949 a, Chao, Tollestrup, 
Fowler and Lauritsen 1949, Strait, Van Patter and Buechner 1950) have shown that the 
result obtained from the range measurement is too low, possibly because the range—energy 
curve is not accurate in this region. The range comparisons mentioned above can now be 
used to recalculate the alpha-particle energies, and hence the Q-values, for the aluminium and 
sodium reactions, because the ranges of the alpha-particle groups emitted in the three 
reactions are close together and only the slope of the range—energy curve over a limited 
region is required. This slope appears, from our range measurements and from other work, 
to be fairly reliable in the 2 Mev. region. ‘The most accurate Q-value available for the 
“F(pa)'*O* reaction is 1-977-+0-008 Mev. (Chao et al. 1949), and with this the energy- 


releases in *7Al(pa)*4Mg and ?*Na(pa)?°Ne become 1:57-+0:03 Mev. and 2:37+0-045 Mev. 
respectively. 
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As a check on these results direct energy measurements of the alpha-particle groups have 
now been made with the aid of a magnetic analyser, calibrated with particles of known energy 
as previously described (Burcham and Freeman 1949 b), and thin targets of Al and NaCl. 
Corrections were made for target thickness and carbon contamination, and the energy of the 
alpha-particles emitted in the reaction ??Al(pax)24Mg at a proton bombarding energy of 
940 kev. and an angle of 135° was then found from several measurements to have the mean 
value 1-991 +0-010 mev., giving a Q-value of 1°585-+0-015 Mev., in agreement with the 
previous figure, and consistent with the value 1:59-+0-07 mev. (French and Treacy 1950) 
deduced from the reactions ??Al(da)?*Mg and 24Mg(dp)2®Mg. For the reaction 28Na(pa)2°Ne 
the energy of alpha-particles emitted at 135° with a proton bombarding energy of 590 kev. 
was measured as 2:28+0-03 Mev., giving a Q-value of 2:34-0:04 Mev.; the probable error 
in this case was larger because of greater uncertainty in the analyser calibration in this region. 
The mean Q-value from the range and energy measurements is 2°35-+0-04 mev. 

Adding the new Q-value for !°F(pa)!®O* and the energy of the gamma-radiation from 
the excited O* gives the energy-release in the reaction !°F(pa)!®O as 8:10-L0-03 mev., 
according to Chao et al. (1949). The energy-release in this reaction has been determined 
here by comparing the range of the alpha-particle group with that from the reaction 
*Li(px)*He and using for the latter the Q-value 17-28+0-03 mev. found by Smith (1939). 
Thin targets of calcium fluoride and oxidized lithium were mounted on a rotatable target 
head so that either target could be irradiated with a proton beam of 840 kev. energy, and the 
alpha-particles emitted at 90° passed through a mica window of 3-9 cm. stopping power into 
a shallow air-filled ionization chamber whose distance from the window could be varied. 
The difference in the extrapolated ranges of the two groups was found to be 2:828-+-0-007cm. 
of standard air. Corrections were applied for straggling and variation in stopping power of 
the mica window; the range-energy relation, which is well-established in the required 
region, was used and a value of 7-075 -+0-030 mev. was deduced for the energy of the alpha- 
particles from the fluorine reaction. This result gives a Q-value of 8:06-0:04 Mev., which 
is consistent with the value quoted above and in good agreement with the recently published 
value (8-068 Mev.) of Strait et al. (1950). 

In the course of this work similar observations were made on the reactions *N(pa)?C 
and O(pa)!®N. Thick targets of aluminium nitride and copper oxide were used, and a 
mica window of 2 mm. stopping power. In each case the range of the alpha-particle group 
was compared with that from a thorium-C source (extrapolated range assumed to be 
4-796 cm.) inserted in place of the target button. The Q-values thereby deduced were 
4-96+0-05 Mev. and 3:97-0:05 mev. for the nitrogen and oxygen reactions respectively,. 
in good agreement with the values of Burcham and Smith (1939). 

The Q-value results presented are summarized below : 


27Al(pa)?4Mg : Q=1-'585-40-015 Mev. 

23Na(pa)29Ne: Q=2:35 +0-04 Mev. 
ME (pa)O: O=8:06 +0-04 mev. 
O(pa)/N: QO=3-97 -£0-05 mev. 
UN(pa)#C : O=4:96 +0:05 mev. 


Part of this work was carried out during the tenure of a Studentship awarded by the- 
Commonwealth Scientific and Industrial Research Organization of Australia. 


Cavendish Laboratory, Joan M. FREEMAN. 
Cambridge. 
15th March 1950. 


Burcuam, W. E., and Freeman, J. M., 1949a, Phys. Rev., 75, 1756; 1949b, Phil. Mag., 40, 807. 

Burcuam, W. E., and Smiru, C. L., 1939, Nature, Lond., 143, 795. 

Cuao, C. Y., ToLiestrup, A. V., Fowier, W. A., and Lauritsen, C. C., 1949, Bull. Amer. Phys. 
~ Soe, 24; Wo: 8, 11: 

FREEMAN, J. M., and Baxter, A. S., 1948, Nature, Lond., 162, 696. 

Frencu, A. P., and Treacy, P. B., 1950, Proc. Phys. Soc. A, 63, 665. 

Ho.tioway, M. G., and sly Sree 1938, Phys. Rev., 54, 18. 

SmituH, N. M., 1939, Phys. Rev., 56, : ; 

Srrair, E. N.; Vane ParrER, D. Ate ana Burcuner, W. W., 1950, Bull. Amer. Phys. Soc., 25,. 

INIOR Bin Bio 


670 Letters to the Editor 


The L-M and K-M Discrepancies in the Rare Earths 


It has long been known that for many of the heavier elements the measured valucs of 
v/R for the Myy and My absorption edges differ from the values calculated from the relations 


My= Ly— Lp, 
My =Lyy— Lae 
My=Lm— Le, 


by amounts which appear to be greater than can be accounted for by experimental error. 
This ‘ discrepancy ’ can only arise from differences in the end levels of the transitions from 
Ly, mr and Myy,y induced by absorption. Where the difference A, given by 


A= Mineas oe Meale> 


is negative and small it-seems possible that it may be accounted tor by the operation of 
selection rules, the end levels being in the neighbourhood of the conduction bands, but in 
cases where it is positive (and large) Phelps (1934), McGrath (1939) and Cauchois (1942) 
have shown that it is necessary to assume transitions into lattice levels beyond the conduction 
levels. In reviewing the data available in 1934 Phelps (1934) showed that A is positive from 
83Bi to 7Ta (the latter being the element of lowest atomic number for which data have 
hitherto been available) and predicted that it would very likely vanish for the rare earth 
elements below 74Lu. In a later review Cauchois (1942) showed that A is very likely 
negative for °?U, °1Pa and °°Th, positive from **Bi to 7°T a (reaching a maximum of more than 
40 ev. round about 74W) and expressed the opinion that it would vanish probably in the 
neighbourhood of 7°Yb. 

During the past two years we have been carrying out an investigation into the fine 
structure of the M;y and My absorption edges of a number of the rare earth elements in the 
form of their oxides and have incidentally measured the wavelengths of the edges themselves. 
By combining these with such L data as are available we have estimated A for the elements 
59Pr, 8°Nd, ®Sm, ®8Eu, ®*Gd and **Dy. Of these only ®*Sm had previously been investi- 
_ gated by Rule (1945), and our results for this element are in close agreement with his. The 
values for A are given in Table 1 in v/R and in electron volts. It will be seen that in all 
cases A is negative and with one exception exceeds 5 ev. The wavelengths of the L edges 
and lines have been taken from the tables of Cauchois and Hulubei (1947), with the exception 
of those for °*Dy where we have taken the more recent values of Barrére (1947). 


Table 1 
Myy My 
A(R) A (ev.) A(R) A (ev.) 
59Py 270-52 = 750. = 0-49 == 6:65 
59Nd —0-66 —8-9, —0-45 — 615 
625m —0:71 —9-6; —0-56 — 7:6, 
8Eu —0:53 —7:26 —0°52 — 7:0, 
4Gd —0-73 —9-9, —0:75 —10-2, 
A Di —0:49 ST OEO, —0:27 — 3-6, 


Obtained as they are by combining our own measurements with those of many different 
“observers it is difficult to estimate the probable error, but we consider that it is unlikely to be 
much more than +5 ev. It would appear likely therefore that from 5°Pr to ®*Dy the dis- 
crepancy though possibly small is negative. In many cases the M edges observed by us 
are accompanied by a ‘ white line ’ on the long wavelength side, due probably to transitions 
to unoccupied levels in the Nyy, yt electron groups, and a correction to the measured value 
of the edge should be made on this account, as described by Sandstrém (1935); since, 
however, the profiles of the lines are not sufficiently well brought out on our photographs we 
have used the uncorrected values in our calculations of A.” The effect of the correction 
would be to increase the wavelength by one or two parts in 10,000, which would have the 
effect of making the A’s slightly more negative. To these results we may add the very 
doubtful case of Yb; here the Myy edge was badly obscured by that of the K-edge of the 
aluminium foil on which the ytterbium was deposited, while instead of the My edge we have 
been able to get only what appears to be a narrow and sharp absorption line, in agreement 
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with an earlier observation of Lindberg (1931). Assuming the long wave edge of this line 


to be the My edge we find the corresponding A to be —9-92 ev. 
For three of the elements studied by us, namely Pr, Nd, and Sm, estimates of the Myy 


and My levels based on K emission data have been given by Rouault (1939). .The discre- 
It will be seen 


pancies between these values and our measured values are given in Table 2. 
that they are all negative and differ from the L-M discrepancies by amounts probably of 


the order of magnitude of the experimental error, three of them being greater than and three 
less than the L—M discrepancy. It would certainly seem unwise at this stage to attempt to 
attach significance to the fact that individual values do show differences. 


Table 2 
M Mv 
AQ/R) —— Aev.) A(R) A ev.) 
=eh Be —1:05 —14-2, —0:96 — 13-0; 
60Nd —0°45 — 61, —0-36 — 4-8, 
82S mm —0-93 —12:6, —0:56 — 7:6, 


The accompanying diagram shows a plot of the L—My discrepancy against Z taken from 


‘Cauchois’ table (1942) with our own results added. A similar plot for the Myy discrepancy 
With such uncertain and limited data any generalization is necessarily 


resembles this closely. 
unsafe, but the diagram does suggest that the ‘ normal’ condition is for the end level in 
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Myy,y absorption to be a tew electron volts below that for Ly,; absorption but that for 
elements 72 to 83 exceptional circumstances mask the negative discrepancy and convert it 
In view of the difficulties likely to be encountered in filling the 


into a large positive one. 
gap between °°Th and ®*Bi, the elements 72Hf and 74Lu promise to be of particular interest 


E. A. STEWARDSON. 


and it is hoped shortly to investigate these substances. 
H. F. Zanpy.* 


Physics Department, 
Leicester University College. 
22nd February 1950. 


* Now at the University of Teheran. 
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Inhomogeneity of Deformation in Metal Single Crystals 


The presence of strain hardening, and the occurrence of asterism in x-ray Laue Pes 
graphs shows that plastic deformation of metal crystals by slip is not merely a simple 
geometrical translation of parts of the crystal along the slip plane. The process is accom- 
panied by local changes of orientation which do not, as the phenomenon of asterism itself 
indicates, amount to a random fragmentation. In some cases asterism has been shown to 
be due to a macroscopic bending of the slip lamellae, but no such evidence has been obtained 
for cubic metals. "The question arises whether there is a connection between asterism and 
the phenomenon of deformation bands (Barrett 1939) or kinking (Orowan 1942) which leads 
to local changes of orientation. An obvious method of studying this question is the X-ray 
technique developed by Berg (1934) and later used for metallurgical investigations by 
Barrett (1945). In what follows the results of an investigation of local orientation changes 
in deformed aluminium crystals by means of a modification of the Berg method will be given. 

The principle of the method is shown in Figure 1. An x-ray beam from a line focus falls 
on a crystal which is rotated till a Bragg reflection occurs. The reflection represents a point 
for point image of the crystal surface, and changes in orientation, e.g. local rotations or 
bending, are revealed as striations in the image. The most satisfactory results are obtained 
on large crystals with flat faces which have been electropolished. The image is detected on a 
fluorescent screen which is then replaced by a Kodak Maximum Resolution Plate placed 
close to the crystal, but shielded from the direct beam. The resolution of the method is 
such that the image can be usefully enlarged 50-100 times, and so microscopic features of the 
deformation can be observed. 


LINE SOURCE 
RAYS 


= 


PHOTOGRAPHIC PLATE 


CRYSTAL 


Figure 1. 


Large aluminium crystals (99:5°% pure) lightly deformed in tension (1-10% elongation) 
were examined in this way, and in all cases fine networks of bands were observed in the 
images. Figure 2 shows an image enlarged to 50 diameters from a crystal after 4°, elon- 
gation. ‘The plane on which the bands form can be determined by examination of two 
surfaces at right angles. The fact that this plane is usually a (110) shows that the bands are 
not slip bands. ‘They are practically invisible in the optical microscope up to elongations of 
about 5% but can be detected after heavier deformation as very nairow deformation or 
kink bands in which the lattice is sharply bent. In the crystals examined they were usually 
about 0-05 mm. apart while their width was about 0:005 mm. Careful adjustment of the 
single crystal enables also the slip lines to be detected (Figure 3) and thus their relation to the 
deformation bands can be studied. Fine deformation bands have been observed in all 
aluminium single crystals so far examined, and are present after as little as 1% elongation 
in tension. ‘They become broader and more clearly visible in the microscope in regions 
where the deformation is less homogeneous, e.g. at the shoulders of the specimen and in the 


vicinity of small included grains, while in polycrystalline specimens, the familiar macroscopic 
deformation bands are often observed. 
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It seems that there is a connection between the occurrence of fine deformation bands and 
the existence ot asterisms in x-ray Laue photographs. The author has found that asterisms 
occur in aluminium single crystals after as little as 1°/ extension, and that in all cases the 
asterism possesses a fine structure which is not due to absorption edges of the crystal or of the 
film. In the continuous smear of the asterism, intensity maxima can be detected; the 
number increases as the size of the beam is made larger. It has been shown that this occurs 
immediately after the deformation and is not due to a recovery process such as polygoni- 
zation, which is only observed in deformed aluminium crystals on annealing at 300° c. or 
higher (Cahn 1949). In the light of the above experiments, this phenomenon can be explained 
if it is assumed that, as a result of the deformation, the crystal breaks up into small blocks or 
layers varying slightly in orientation, which produce the dark spots in the Laue asterisms. 
The blocks are joined by narrow bent regions which give rise to the continuous background 
of the asterisms. It seems clear that these regions provide the visible curvatures as distinct 


igure 2. Aluminium crystal—4% elongation in tension. Figure 3. Aluminium crystal—4% elongation in tension. 
X-ray photograph (x 37-5). (Cu radiation). X-ray photograph (x 37:5). (Cu radiation). 


from the hypothetical ‘local curvatures’ postulated by Yamaguchi, Taylor and Burgers 
as the cause of x-ray asterisms. It is perhaps relevant to mention that the author has found 
that with cadmium single crystals of suitable orientation extended by 100% in tension, the 
Laue spots can remain practically free from asterism. Cadmium crystals deformed in 
tension when examined by the above x-ray method show no deformation bands. 

It is well established that deformation bands or kinking plays an important part in the 
deformation of polycrystals. The above results show that this mode of deformation is 
significant also in the tensile deformation of single crystals, and thus even in these circum- 
stances, the plastic deformation is often far from homogeneous. 

The author wishes to thank Dr. E. Orowan, for helpful advice and criticism, and acknow- 
ledges an I.C.1. Fellowship during the tenure of which this work was begun, and a Fellowship 
from the Armourers and Braziers’ Company which is enabling its continuation. 


Cavendish Laboratory, R. W. K. HoNEYCOMBE. 
Cambridge. 
28th February 1950. 
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Singlet System B of ZrO 


In a previous communication (Afaf 1949) the author drew attention to the discovery of 
a 1X1 system, since designated as A, in the spectrum of ZrO. This system is in the 
ultra-violet region, with the (0, 0) band head at A 3682. The vibrational constants in the 
lower aX and the upper b!= states were evaluated from the data of nineteen band heads. 
The constants in the b!X state were we’ = 843-75 cm and xe’we’ =3-15 cm, and those in 
ald state were found to be 978-07 and 5:04 cm7! respectively. Further details are given in 
a later article (Afaf 1950). It has been shown that the a'Z state lies somewhere about 
1,300 cm-! above the x°II, ground state. The b* state should therefore be about 
28,451 cm~! above the x*II, state. Further efforts have now been made to assign the 
remaining unclassified bands. The most prominent heads of these are at AA 4460-4 (2), 
4736-9 (12), 5185-0 (7), 5553-1 (10), 5610-1 (10), 5629-5, (14). 

The author has now established the identity of another singlet system for ZrO, of which 
the A 5185 band head is the (0,0) band. This system has been designated as system B. 
In all only nine band heads were observed for this system. ‘The record of the various 
band heads as measured by the author in the region of the bands of system B, along with the 
analysis offered, is presented below. It will be noted that a few more bands have been 
recorded for triplet system «a of ZrO, extending the system to the longer wavelength. 


Noi I Vee Classification 
4736-91 (6) 21105:0* Several band heads of system « (AA 4730-4900), given 
4769-41 (1) 20961 -0* elsewhere (Lowater 1932), are omitted here. 
4798-91 (0) 20832 -3* B (3, 1) 
4863-09 (3) BYOB 
4895-31 (1) 20422-0 
4969-82 (2) 20115-9 133 (605 (0) 
4996-28 (2) 20009 °3 B (2, 1) 
5022-95 (1) 19902-9 BiG; 2) 
5052-95 (1) 19785-6 « (0, 2) Ra 
5073-94 (1) 19703-1 « (0, 2) Rp 
5077-60 (1) 196889 « (0, 2) Re 
5102-28 (0) 195936? « (1,3) Rp? 
5103-77 (0) 19587-9* « (1, 3) Re 
5128-80 (0) 19492°3 « (2, 4) Rp 
5131-67 (0) 19481-4 (2, 4) Re 
5155-39 (0) 19391°8 « (3, 5) Rp 
5156-28 (1) 19388 -4* 

5158-64 (0) 19379-6 a (3,5) Re 
5185-04 (6) 19280-9* B (0, 0) 
5205-65 (0) 19208°3 

SPAWaIy (1) 19180-5* BiGlead) 
5240-17 (0) 19078 -1 y B (2, 2) 
5436-87 (3) 18387-8* B (3, 2) Xa 
5438-23 (2) 18383 -2 

5439-47 (3) 18379-0* y (5, 2) Ra? 
5448-88 (1) 18347°3 B (0, 1) 
5456:56 (6) LSS2ie 5% B (1, 0) Re 
5461-55 (4) 18304-7* B (A, 0) Xe 
5465-38 (2) 18291 -9* B (4, 3) Xp 
5478-39 (0) 18248°5 Biiee2) 
5485-74 (7) 18224-0* B (2, 1) Re 


The bands marked with an asterisk (*) were also recorded by Miss Lowater (1932). 


The Deslandres ‘Table of system B is given below. 


e 
(@) UA 2 uv / 

aN A’1/ 1 A’3/2 2 

A 0 19280-9 (933-6) 18347°3 
412 (835-0) (833-2) 834-1 
1 20115-9 (935-4) 19180-5 (932-0) 18248-5 : } 

A’3/2 (828-8) (829-6) (829-2) 
2 200093 (931-2) 19078-1 

A’5/2 (823-0) (824-8) (823-9) 
3 20832°3 (929-4) 19902-9 


(934-5) (930-9) 
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The vibrational differences A’(v’+ 4) in the upper state of system B, viz. 834+1, 829-2 
and 823-9 cm, agree very closely with those in the upper state b!X of system A, viz. 
836-6, 830-1 and 823-9 cm~!. This shows that the upper state of system B is bE. Since 
the band heads of system B are definitely single and no trace of Q branches can be found, it 
follows that the lower state of system B should also be 1, and is henceforth termed cD. 
The state c'% (v”=0) should lie about 9,170 (28,451-19,281) cm~! above the x®I1, ground 
state. 

System A and system B involve respectively the following transitions ; b!% — a!X and 
b*X +c'X. The possible intercombination system, viz. c!& — a1, will lie in the infra-red 
region near about 13,000 a., and should be observed with some difficulty because of the 
somewhat weak intensity in the bands. 

A further point of interest is the fact that, while at a temperature of 2400 a. enough 
molecules might be maintained in the a! state (1,300 cm~ above x*I,) to cause absorption 
to take place in the Late Type Stars in the bands of system A, such an absorption is extremely 
unlikely to take place in the bands cf system B on eccount of the greater excitation of the 
lower level (9,170 cm™ above x*II,). Absence of absorption in the spectra of Late Type 
Stars in bands of system B would be consistent with the present analysis. 

The formula relating band heads of system B is 


vi” =19330-5+ {839-2 (v’ +4) —2-56 (v’ +4)?} — (938-1 (v” +3) —1-80 (v7 +)". 
The investigation is being continued. 


‘The author expresses his gratitude to Assistant Professor R. W. B. Pearse for valuable 
guidance throughout this work. 


Department of Physics, M. AFar. 
Imperial College of Science and Technology, 
London, S.W.7. 
27th March 1950. 


Arar, M., 1949, Nature, Lond., 164, 752; 1950, Proc. Phys. Soc. A, 63, in the press. 
LowatTer, F., 1932, Proc. Phys. Soc., 44, 51. 


REVIEWS OF BOOKS 


Theory of Atomic Nucleus and Nuclear Energy-Sources, by G. GAamMow and 
C. L. CrircHFIELD. Pp. vit+344. 3rd Edition of Structure of Atomic 
Nucleus and Nuclear Transformations. International Series of Monographs 
in Physics. (Oxford: Clarendon Press, 1949.) 30s. 


This book is presented as the third edition of Gamow’s pioneer treatise of 1931 on the 
structure of the atomic nucleus. Quantum mutatus ab illo! It was then possible in a 
slender volume to give an exhaustive account of our knowledge of nuclear properties. 
Now, with its size doubled, the book brings only a selection of the most important topics 
of nuclear physics in a much condensed form. The senior author, moreover, felt the need 
of calling a collaborator to the rescue, and, as in the second edition, the treatment has more 
the sedate character of a textbook than had that survey of a strange and novel world which 
the first edition presented. However, the spirit of adventure, which gave this early attempt 
its special charm, has not entirely disappeared : the cosmogonical chapter, which forms 
one of the most valuable additions of the present issue, is as Gamovian as one could wish. 

The difficult task of selecting and arranging the subject-matter, and of reducing the 
account of the various aspects presented to their most essential features, has on the whole 
been successfully accomplished by the authors ; they have produced a readable, well- 
balanced and authoritative survey of the many-sided subject. ‘There are, however, some 
weaknesses in the composition of the book and in the presentation of some questions. 
For instance, chapters III and VIII, which deal with nuclear forces and with nuclear 
collisions, are rather sketchy and confused, and the way in which the all-important concept 
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of the compound nucleus in nuclear reactions is introduced, is very unsatisfactory. When 
discussing the proton and neutron as elementary nuclear constituents, the mistaken 
impression is conveyed to the unwary reader that the conception of these particles as distinct 
species is somehow essentially different from their description as two states of the * nucleon ’ 
species. Apart from such minor criticisms, one can but praise the clever selection of the 
illustrative examples from the experimental material, and the reliability of the information 
given. In fact, it is so difficult to catch the authors dozing that I can only point to a single 
wrong statement, viz. that the symmetrical meson theory would not be consistent with the 
saturation properties of the nuclear forces. 

On the other hand, I must violently object to the first pages of the book, purporting 
to be an ‘ historical’ introduction. It is hardly conceivable that scientists (our authors 
are not the only sinners in this respect) loath to put down a single technical statement 
without carefully checking its accuracy, can allow themselves to be so lax when dealing 
with historical matters. The result of this attitude in the present case has been to produce 
a startling hotch-potch of Democritean atoms and Aristotelian elements, and an unfortunate 
confusion between the French chemist Proust and the British physician Prout. 

Fiom another point of view, it is a matter for regret that the authors should uphold the 
illogical nomenclature electron—positron, instead of adopting the correct usage of the word 
electron for any sign of the charge, and the words negaton—positon to distinguish between 
these signs. This terminology has been lately recommended by the Cosmic-ray Commission 
of the International Union of Physics, and it is a bad example indeed to show so little 
regard for the efforts, however modest, made by scientists to put their own house in order. 

L. ROSENFELD. 


Introduction to Atomic Physics, by S. ToLaNnsky. Pp. x1+371. 3rd Edition. 
(London: Longmans, Green and Co., 1949.) 15s. 


The third revised and enlarged edition of this popular textbook for undergraduate 
students is welcome. ‘I‘wo new chapters have been added, one on the neutron and the 
second on nuclear fission and atomic energy. There are new sections dealing with the 
' radioactive estimation of the age of the earth, K-capture and isotope separation. The 
chapter on cosmic rays has been revised, with sections on the photographic plate technique 
and mesons. 

In a book of manageable size and such wide scope, dealing equally with the electronic 
and nuclear properties of the atom, each section must necessarily be sketchy and brief. 
‘The whole book integrates the present knowledge and outlook and gives a perspective of this 
field which makes it an excellent introductory text. W. I. B. SMITH. 


Higher Physics, by E. NIGHTINGALE. Pp. xvi+808. (London: G. Bell and 
Sons Ltd., 1948.) 27s. 6d. Also in separate parts: Pt. I. Mechanics and 
Properties of Matter, 7s. 6d.; Pt. I. Heat, 7s. 6d.; Pt. III. Light and Sound, 
10s); Pty LVs Electricrivel2s6d: 


The clarity of style which has characterized Mr. Nightingale’s many school texts in the 
past is in evidence in his latest effort. This book is of intermediate standard and is primarily 
intended for sixth form scholars. 

The wide field of subject-matter which is covered within the confines of this single 
volume has necessitated considerable ‘ condensing ’ in certain sections, but on the whole this 
process has been carried out quite efficiently. However, it is inevitable that some sub- 
divisions become unduly compressed when room is found for such topics as the theory of 
flight, the quantum theory of spectra, etc. 

An undesirable effect of the lag, so prevalent in the post-war era, between the preparation 
of an MS. and its publication is exemplified in the present instance by the author’s failure 
to record the recent official acceptance (October 1946) of the absolute system of electrical 
units as the practical system, to the exclusion of the International Scale. Also no mention 
is made, rather surprisingly, of the M.K.s. system of units. 
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‘The author makes frequent use of the methods of dimensions in the various sections of 
the book and the use of graphs in exhibiting experimental results is also stressed in the text. 

The errors noted are remarkably few, but one which should be noted is the mention, 
on page 633, of the use of an iron former for the suspended coil of a galvanometer. 

Line diagrams which are capable of reproduction by the student are freely employed 
thoughout the text, and another feature which commends the book as one to be treasured by 
the scholarship candidate is the large selection of examination questions at the end of each 
chapter. R. W. B.S. 


A Textbook on Heat, by J. H. Awprry. Pp. x-+302. 1st Edition. (London, 
New York and Toronto: Longmans, Green and Co., 1949.) 15s. 


This is a book which has been written because the author, who is a specialist in the 
subject, has a conviction that heat is one of the fundamental parts of physics and therefore 
all physicists will find something of value in a fresh study of its principles. 

The first seven chapters give an excellent account of the definition and measurement of 
temperature and of modern calorimetric practice which is based upon first-hand knowledge. 
The author then turns to the more theoretical topics, introducing the chapters on thermo- 
dynamics with a clear explanation of the meaning of partial differentiation which ensures 
that the mathematics which follows does not obscure the physics of the problems. Un- 
fortunately in some of the later chapters this clarity is lost. In Chapter XVI, which deals 
with the conduction of heat, the author includes a large number of special mathematical 
methods of solving certain problems of heat flow. In some cases the account is too brief, 
for instance, half a page scarcely gives much idea of the method of the Fourier transform. 
In others there is a tendency to lead the reader half-way through a lengthy piece of analysis 
and there to leave him. 

Throughout the book the mathematics is seldom set out in such a way as to be easily read, 
being too closely interwoven with the text. There is also a shortage of diagrams, and the 
choice of these at times appears a little strange; in the first three chapters there is only one, 
and that is of amaximum and minimum thermometer, while further on, on p. 84, a graph taken 
from some published paper illustrates, not very clearly, a small point in the text. These are 
probably matters of economy which may be remedied in a second edition. 

On the whole this is a book of an original character which will make a useful supplement 
to the existing textbooks. It would, in particular, make refreshing reading for students at 


the end of a degree course, but there is something in it for readers at all stages. 
E. M. WILLIAMSON. 


Die Methoden zur angenaherten Losung von Eigenwertproblemen in der Elasto- 
kinetik, by K. HoHENeMsER. Pp. 89. 2nd Edition. (Berlin: Springer- 
Verlag., 1943; New York: Chelsea Publishing Company, 1949.) No 
price. 


In the classical theory of elasticity, the eigenvalue problem arises when a function y 
(representing, for example, displacement or stress or strain) is determined by a linear, 
homogeneous differential equation involving variables x,, x2, ..-. X,, the coefficients of y 
and its derivatives being given as functions of the variables x, and of certain parameters 
d,, Ay, -.+-A,. In addition, a number of boundary conditions for y and its derivatives are 
given, and the eigenvalue problem consists in finding the values of the parameters A, which 
satisfy both the differential equation and the boundary conditions. In practice, it 1s con- 
venient to subdivide problems of this type into two categories : (a) statical problems in 
which the variables do not include time, e.g. Euler’s treatment of the buckling of a thin rod, 
and (6) dynamical or vibration problems, in which time enters, e.g. Bernoulli’s solution of 
the transverse vibration of a stretched string. In the former case, the parameters A, are 
usually associated with critical loads or stresses at the onset of instability, and, in'the latter, 

‘with the frequencies of free vibrations of the system. 

The aim of Professor Hohencmser’s monograph is to survey the modern mathematical 

techniques available for the solution of the dynamical problems, and the scope of the work 
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may be judged from the specific cases which are treated as examples of the general methods : 
torsional and transverse vibrations of bars of various forms of cross section, the critical 
speeds of rotation of shafts, the vibrations of simple types of frame-structures, and the 
transverse vibrations of plates. 

The first part of the book is devoted to a review of the general methods : the integral 
equations of the vibrations of elastic bodies, including induction (or influence or Green’s) 
functions; iteration processes and the methods of extreme values for the calculation of 
eigenvalues; differential equations methods; replacement of continuous elastic system by 
a system with concentrated inertia and restoring forces; the solution of frequency 
determinants. This section is illustrated by treating the specific examples mentioned 
already. 

The book is exceedingly well written and the treatment throughout is very general in 
character; it has the usual advantages and disadvantages of this approach. On the one 
hand, the general framework of the subject is clearly shown and it forms an admirable 
introduction for readers who are mathematically inclined; on the other hand, it is apt to 
prove disappointing to a practising Physicist or Engineer who requires information relating 
to the solution of a particular problem. A very valuable feature of the book is the extensive 
bibliography of the subject, although it is rather surprising to find that Rayleigh’s Principle, 
by G. Temple and W. G. Bickley (Oxford University Press 1933), is not mentioned; the 
value of the bibliography would have been increased had an author index been included. 

R. M.D. 
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ABSTRACTS FOR SECTION B 


The Behaviour of the Cathode Spot on an Undisturbed Liquid Surface of Low Work 
Function, by K. D. Froome. 


ABSTRACT. The cathode spot of transient arcs using the liquid sodium—potassium alloy 
for cathode has been studied by means of the Kerr cell apparatus capable of taking a 
sequence of photographs. It is found that for a given current the spot takes, as in the case 
of mercury, the form of a line or broken line of total length roughly proportional to the 
current. ‘This is also true for a spot in a magnetic field. Again as for mercury, if the 
current rises slowly the length of line increases proportionally to the current. If the rate of 
growth of current is greater than a value lying approximately between 10® and 2 1068 
amp/sec., then fresh spots form and spread out radially from the point of formation into thin 
semicircular or circular lines moving with a maximum radial velocity of about 104 cm/scc. 
The apparent current density of emission lies between 2:5 and 5x 10% amp/cm2, and 
consequently is higher than for mercury. 


The Calculation of Magnetic Lens Fields by Relaxation Methods, by Mary B. 
HEsse. 


ABSTRACT. ‘The field and lens constants are calculated for two typical magnetic lenses 
as used in electron microscopes, using the relaxation method developed by Southwell for the 
solution of potential problems. The results are compared with those based on simple 
analytic approximations to the field distribution in the work of Glaser and Ramberg, and 
are shown to agree closely as regards those characteristics on which the discussion of the 
performance of the lens is based. 


An Electrostatic Fluxmeter of Short Response-time for use in Studies of Transient 
Field-changes, by D. J. Maan and B. F. J. SCHONLAND. 


ABSTRACT. ‘The fluxmeter consists of a conducting system of small capacity which is 
alternately exposed to and screened from the electric field by the movement of a rapidly 
rotating metal shield, thus producing an alternating E.M.F. of period 0°83 millisecond. ‘The 
output, after amplification, is displayed on the screen of a cathode-ray oscillograph and the 
record carries an automatic indication of the sense of the field every 7:5 milliseconds. 

At the maximum sensitivity employed the device gives a deflection of one centimetre in a 
field of 20 volts/metre, with a background noise-level of 0:15 cm. (3 volts/metre). 

The response-time of the instrument is one half-cycle (0°42 millisecond). 

The fluxmeter has been developed to give the rapid response and high sensitivity 
required for studies of the electric fields of thunderstorms in the intervals between the 
separate strokes of a lightning discharge. 


Wave Propagation in a Slipping Stream of Electrons : Small Amplitude Theory, 
by G. G. Macrarvane and H. G. Hay. 


ABSTRACT. A theoretical study is presented of the TTM-waves that can travel along a 
slipping stream of electrons. A slipping stream is defined as one in which the electrons 
move in parallel paths with velocity which varies with distance transverse to the motion. 
It is found that amplifying waves can travel along a slipping stream for all frequencies. It 
is also found that the slipping-stream tube has the remarkable property of combining the 
characteristics of a two-beam tube and a travelling-wave tube. This is due to the occurrence 
in the stream under suitable conditions of resonance layers, which act as highly reactive 
impedance sheets and can guide 'TM-waves of slow phase velocity in the same way as the 
spiral or corrugated surface in a travelling-wave tube. 

The slipping-stream tube behaves differently according to the fractional velocity range 
of the electrons. If the velocity varies linearly from Vy to V_4 across the stream and 
a=(V_g—V_@|Vat V_«)<0-42 it behaves primarily as a two-beam tube with plasma 
resonance frequency w /2, and therefore having a cut-off frequency W,=W/aV/2 and 
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maximum gain 2*1w»/V5 decibels per unit length, where Wo 1s the plasma resonance frequency 
of the slipping stream; in addition it has a low gain for w >a). However, if «>0-42 it 
behaves primarily like a travelling-wave tube and the rate of gain is about 0-53w/V—a 
decibels per unit length at all frequencies above the plasma resonance frequency. ‘This 
rate of gain is low compared with the gain possible with a travelling-wave tube in which the 
electron velocity is V_, but it is achieved without the need for any external slow-wave 
waveguide. 

Wave propagation along a slipping stream inside a waveguide, which in the absence of 
the electrons can guide a TM-wave of slow phase velocity, is also discussed. It is shown 
that this slipping stream travelling-wave tube has very much the same characteristics as a 
travelling-wave tube with a uniform electron beam. Maximum gain occurs when the phase 
velocity Uv» of the wave in the empty guide is about equal to the velocity of the electrons 
nearest to the reactive impedance sheet of the guide—in the example considered this is the 
fastest electron velocity V,. A parametric set of curves is given relating the complex 
propagation constant to the frequency for different ratios U9/V q. 


Theory of the Production of an Ionized Layer in a Non-Isothermal Atmosphere 
Neglecting the Earth’s Curvature, and its Application to Experimental Results, 
by J. A. GLEDHILL AND M. E. SZENDREI. 


ABSTRACT. A new theory of ionospheric layer formation is developed, in which the 
temperature is assumed to vary linearly with height. The equations are compared at each 
step with those obtained by Chapman in his theory of layer formation in an isothermal 
atmosphere. The equations for the maximum of electron density and its height are also 
given. The effect of the parameters on the shape of the layer is shown in graphical form. 

The equations are somewhat complex in form, but an ingenious graphical method has 
been devised suitable for the application of the theory to results given in the form of a 
Booker and Seaton parabolic distribution of electronic density with height. By application 
of the theory to mean hourly values for four summer months in South Africa, values are 
obtained for the temperature gradient, the t2mperature at 200 km., and its variation over the 
middle part of the day. The results obtained are in accordance with previous estimates and 
offer numerical confirmation of the theory that the atmosphere expands bodily upwards 
during the middle part of a summer day. 


A New Method of Measurement of the Variation with Wavelength of the Refractive 
Index and Absolute Stress Optical Coefficients of Amorphous Solids, by 
J. E. H. Braypon. 


ABSTRACT. White light is passed through a Michelson interferometer in one of whose 
arms a thin lamina of the substance in the shape of a tension member is supported by a 
loading frame. The emergent beam, examined by a spectrometer, shows a continuous 
spectrum crossed by interference fringes. For the unstressed lamina the position of the 
fringes gives the variation of refractive index with wavelength. For the stressed lamina the 
behaviour of the fringes gives the stress optical coefficients. No modification of the inter- 
ferometer, and only three sets of observations on the fringes, are required. For the calcu- 
lation of absolute values of the refractive index the values of two wavelengths must be 
known. : 

Full details of the method are given, and photographs of tests on Columbia Resin 39 are 
included. ‘These show that CR39 becomes positively birefringent under tension. 


Vibrations of Free Elliptical Plates, by Mary D. WALLER. 


ABSTRACT. The systematic observations which have been made on the normal 
vibrations of free elliptical plates confirm Chladni’s results. "The nodal systems may be 
divided into the same four classes of symmetry as those of rectangular plates. 

The correspondence which exists between the vibrations of elliptical, rectangular and 
circular plates is traced. Combined modes of vibration do not occur on elliptical plates as 
they occasionally do on rectangular plates. ; | 
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